PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

Similar documents
Chapter 4. Fourier Series

Castiel, Supernatural, Season 6, Episode 18

Today in Physics 217: separation of variables IV

UNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014

Question 1: The magnetic case

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

Solutions to Final Exam Review Problems

Null Spaces, Column. Transformations. Remarks. Definition ( ) ( ) Null A is the set of all x's in R that go to 0 in R 3

Most text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t

Løsningsførslag i 4M

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

MATH2007* Partial Answers to Review Exercises Fall 2004

A) is empty. B) is a finite set. C) can be a countably infinite set. D) can be an uncountable set.

The Non-homogeneous Diffusion Equation

Zeros of Polynomials

A second look at separation of variables

PAPER : IIT-JAM 2010

SECTION 2.6 THE SECOND ALTERNATIVE

Orthogonal Functions

Assignment Number 3 Solutions

(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

CHAPTER 5. Theory and Solution Using Matrix Techniques

Mth 95 Notes Module 1 Spring Section 4.1- Solving Systems of Linear Equations in Two Variables by Graphing, Substitution, and Elimination

MATH 10550, EXAM 3 SOLUTIONS

REVIEW OF CALCULUS Herman J. Bierens Pennsylvania State University (January 28, 2004) x 2., or x 1. x j. ' ' n i'1 x i well.,y 2

TEMASEK JUNIOR COLLEGE, SINGAPORE JC One Promotion Examination 2014 Higher 2

Calculus 2 Test File Spring Test #1

Quantum Mechanics I. 21 April, x=0. , α = A + B = C. ik 1 A ik 1 B = αc.

Kinetics of Complex Reactions

PROBLEMS AND SOLUTIONS 2

Lesson 10: Limits and Continuity

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

Objective Mathematics

Chapter 9: Numerical Differentiation

Ray Optics Theory and Mode Theory. Dr. Mohammad Faisal Dept. of EEE, BUET

3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials

6.3 Testing Series With Positive Terms

Numerical Methods in Fourier Series Applications

A.1 Algebra Review: Polynomials/Rationals. Definitions:

Appendix: The Laplace Transform

ECE 901 Lecture 4: Estimation of Lipschitz smooth functions

Calculus 2 - D. Yuen Final Exam Review (Version 11/22/2017. Please report any possible typos.)

FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING

Chapter 10: Power Series

COMP 2804 Solutions Assignment 1

It is often useful to approximate complicated functions using simpler ones. We consider the task of approximating a function by a polynomial.

( ) = p and P( i = b) = q.

7.) Consider the region bounded by y = x 2, y = x - 1, x = -1 and x = 1. Find the volume of the solid produced by revolving the region around x = 3.

Chapter 10 Partial Differential Equations and Fourier Series

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Sequences and Series of Functions

IIT JAM Mathematical Statistics (MS) 2006 SECTION A

Properties and Tests of Zeros of Polynomial Functions

Linearly Independent Sets, Bases. Review. Remarks. A set of vectors,,, in a vector space is said to be linearly independent if the vector equation

Numerical Methods for Ordinary Differential Equations

Honors Calculus Homework 13 Solutions, due 12/8/5

Markscheme May 2015 Calculus Higher level Paper 3

MATH 6101 Fall 2008 Newton and Differential Equations

Engineering Analysis ( & ) Lec(7) CH 2 Higher Order Linear ODEs

Solution: APPM 1360 Final Spring 2013

INMO-2018 problems and solutions

September 2012 C1 Note. C1 Notes (Edexcel) Copyright - For AS, A2 notes and IGCSE / GCSE worksheets 1

Mathematics Extension 1

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

Phys 6303 Final Exam Solutions December 19, 2012

Subject: Differential Equations & Mathematical Modeling-III

Unit 4: Polynomial and Rational Functions

CHAPTER 10 INFINITE SEQUENCES AND SERIES

GRAPHING LINEAR EQUATIONS. Linear Equations ( 3,1 ) _x-axis. Origin ( 0, 0 ) Slope = change in y change in x. Equation for l 1.

MATHEMATICS. 61. The differential equation representing the family of curves where c is a positive parameter, is of

Lecture 11. Solution of Nonlinear Equations - III

[ 11 ] z of degree 2 as both degree 2 each. The degree of a polynomial in n variables is the maximum of the degrees of its terms.

RADICAL EXPRESSION. If a and x are real numbers and n is a positive integer, then x is an. n th root theorems: Example 1 Simplify

Math 234 Test 1, Tuesday 27 September 2005, 4 pages, 30 points, 75 minutes.

Orthogonal Function Solution of Differential Equations

5.6 Binomial Multi-section Matching Transformer

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

Section 11.8: Power Series

ENGI 9420 Engineering Analysis Assignment 3 Solutions

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Chapter 2 The Solution of Numerical Algebraic and Transcendental Equations

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

A widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α

MAT1026 Calculus II Basic Convergence Tests for Series

SECTION 2 Electrostatics

MATH spring 2008 lecture 3 Answers to selected problems. 0 sin14 xdx = x dx. ; (iv) x +

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

Different kinds of Mathematical Induction

Flight and Orbital Mechanics. Exams

Class #25 Wednesday, April 19, 2018

Lecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting

Chapter 6 Infinite Series

Lecture #5: Begin Quantum Mechanics: Free Particle and Particle in a 1D Box

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii

Math 110 Assignment #6 Due: Monday, February 10

Math 113 Exam 3 Practice

Transcription:

Diola Bagayoko (0 PARTAL DFFERENTAL EQUATONS SEPARATON OF ARABLES. troductio As discussed i previous lectures, partial differetial equatios arise whe the depedet variale, i.e., the fuctio, varies with ore tha oe idepedet variale. such cases, partial ψ dψ derivatives ( - as opposed to total derivatives ( - appear i the differetial equatios. d Please refer to the lecture o Maxwell s E & M equatios ad the wave, Poisso, ad Laplace equatios we derived fro the for further illustratio. For exaple, the Laplace equatio for the scalar potetial of E & M is y, z or y, z y, z z The separatio of variales that follows is differet fro the oe we used for solvig Newto s type equatios, d x( t d = f ( t, = f ( t, or d = f ( t dt dt dt t is also differet fro the oe we used for solvig hoogeeous first order differetial equatios (that are liear, y P( x y y = P( x y = P( x dx y oth the aove exaples, we separated the depedet ( or ad the idepedet (t or x variales. We gave the ae of first separatio to this for of separatio of variales. The oe discussed elow cosists of separatig the idepedet variales ( x, y, z, or t as i the Laplace equatio aove. t is essetial to ote that the geeral separatio of idepedet variales is oly the first step i solvig partial differetial equatios. This separatio leads to ordiary differetial equatios that are solved (a y usig the first separatio followed y itegratio, ( y utilizig e pt or e px sustitutio ethod for liear, d order, ordiary differetial equatios with costat coefficiets. (c y eployig the Froeius series ethod for ordiary differetial equatios (O.D.E for which the coefficiets are ot costat, ad

Diola Bagayoko (0 (d y usig the solutios otaied (as i a.,., ad c. aove to satisfy oudary or iitial coditios. As show through the exaple elow, satisfyig oudary coditios ofte leads to ifiite series!. The Separatio of ariales Through a exaple. a. Preliiaries A plate of legth a ad width is coected i a circuit as show. Solve the Laplace equatio for this plate. ( at y ad at y =, it is ad at x ad x = a y Solutio a x Due to the geoetry of the plate, we select the (x, rectagular Cartesia coordiate syste, as show aove, to express ad solve the prole. Due to this choice, the scalar potetial is = (x, ad the Laplace equatio is everywhere iside the plate. At the edges of the plate, the oudary coditios (.c. are: at y ad at y = (x, 0 ad (x, = at x (0, = = at x = a (a, = Key: Boudary coditios (.c. for which the depedet variale (the fuctio is zero are called hoogeous.c. Exaples are (x, 0 ad (x,. Boudary coditios that are ot hoogeeous are said to e ihoogeeous or o-hoogeeous. Exaples aove are (0, = ad (a, =. (cidetally, see your E & M lessos to uderstad that the groudig of the upper ad lower edges of the plate eas that the potetial is zero at those edges!. Also ote -- icidetally -- that coditios at extreal values of space variales like (x, y, z or (r, θ, ϕ are referred to as oudary coditios. Coditios for t or t = t 0 are called iitial coditios. Now that we have selected a coordiate syste, idetified the oudary coditios, let us proceed with the separatio of variales.

Diola Bagayoko (0.. Separatio of variales: asic steps [.0] a Assue (x, = ( x ( [.], i.e., is assued to e a product of fuctios ad each of which depeds oly o oe variale. Note well that this assuptio is ot always correct. Hece, whe it is ot correct, the the separatio of variales is ot possile. Get ( x = [ ] ( ad ( = ( x[ ] c Sustitute the i the D.E. ad divide y ( x ( = ( * ( x *. The first ter depeds oly o x ad the secod o y oly. So, the separatio of variale is possile ad we ca cotiue with the search of solutio y this ethod. d The last lie aove is possile oly if = k ad = k. [.] deed, if they are ot costat, the, for exaple, will take a value as x varies- - such that 0, give that x does ot appear i (ad does ot chage. Give that each ter depeds oly o oe variale, we ca switch to total derivatives: d dx d Further, we ust have k k. Hece, oe is positive ad the other oe is egative. We deterie, give the.c., which oe should e positive or egative. 3

Diola Bagayoko (0 Key: Select the sigs of the costats i such a way that you otai oscillatory (cos, si, e ix, e it fuctios i the directio of the hoogeeous.c. ad i the directio of ost cyclic variales [i.e., ϕ, of the π cycle, i ( ρ,ϕ, z cylidrical coordiate syste. our case aove, we eed oscillatios i the directio of y (the.c. for the extreal values of y are hoogeeous. To have oscillatios i y we eed the solutio of d = k to e oscillatory: " k = or " 0 k = px Fro the results of the lecture of e sustitutio, k ust e positive, or k ust e ixy egative, for to e e or cos( y, si( y where is a costat. Ay egative uer k ca e writte as ius a square (all squares are positive or zero. So, set k = -k. d d Hece, we ow have = k ad = k dx These two equatios are d order, ordiary differetial equatio with costat coefficiets. (We have leared how to solve the with e py or e px sustitutio ethod. Hece, these equatios give = Acos( k B si( k ad = ce De, with =, as their geeral solutios. The geeral solutio of the Laplace equatio is the = [ Acos( k B si( k]*[ Ce De ] We ust look for the particular solutio that satisfies the oudary coditios (.c.. Key: Always satisfy the hoogeeous.c. first. So we wat (x,0=0 ad (x,=0 or 0 = [ A* B * 0][ Ce De ] A or = B si ky[ Ce De ]. B si k[ Ce De ] si k or π si k k = π or k =. Note well the apparitio of! πx πx This value of k leads to (x,= B si( [ Ce De ]. Key: With the appearace of iteger, the questio ecoes the followig: Which value(s of do we take to write πx πx (x,= B si( [ Ce De ]? Ad the aswer is - - all of the satisfy the hoogeeous oudary coditios! So, we have to take the all! Hece, (x, is a ifiite series ad πx πx = = B si( [ Ce De ]. [.3] Notice the idices, eve for the coefficiets B, C ad D. 4

Diola Bagayoko (0 Let B C = ad B D = (We are just reaig costats, the πx πx = = ( e e si( Let us ow cosider the ihoogeeous oudary coditios: ( At x = = = si(. Usig the properties of orthogoal fuctios (i.e., i this case si( si * ( 0 π w( ad itegrate fro 0 to : y with W(= over the iterval [0, ], we ultiply y πy πy si*( = si*( ( si( [.4] The itegral o the left had side of Equatio [.4] is πy cos( π 0 which is πy cos( π 0 = π π [ ( ] [ ] for for [.4.] The right had side of Equatio [.4] is ( N δ = ( [.4.]. See the chapter o Fourier Series to otai N = (N is π as π is the rage, so it is if is the rage is i Equatio [.4]. Coiig [.4.] ad [.4.] leads to 4 ( = for odd ad 0 for eve. [.5] π At x = a We have πa πa ( a, = = ( e e si(. = Usig the orthogoal properties of si( ad proceedig as doe aove for x, we get odd eve 5

Diola Bagayoko (0 = a πa 4 for odd π e e [.6] 0 for eve Equatios [.5] ad [.6] for a syste of algeraic equatios with two ukow, ad. Solvig this syste for ad as = ( e odd πx, which we call or relael as e si( ad, gives potetial [You have to solve for ad ad rewrite the aove fial solutio usig their values i ters of,, a,, ad π] 6