Strauss PDEs 2e: Section Exercise 4 Page 1 of 5. u tt = c 2 u xx ru t for 0 < x < l u = 0 at both ends u(x, 0) = φ(x) u t (x, 0) = ψ(x),

Similar documents
CHAPTER 4 FOURIER SERIES

Chapter 4. Fourier Series

Strauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.

Here are some solutions to the sample problems concerning series solution of differential equations with non-constant coefficients (Chapter 12).

Question 1: The magnetic case

Assignment Number 3 Solutions

April 1980 TR/96. Extrapolation techniques for first order hyperbolic partial differential equations. E.H. Twizell

Topics in Fourier Analysis-I 1

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

Chapter 10 Partial Differential Equations and Fourier Series

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Løsningsførslag i 4M

PROBLEMS AND SOLUTIONS 2

Self-Consistent Simulations of Beam and Plasma Systems Final Exam ( take-home )

Discrete Fourier Transform

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

Solutions to Final Exam Review Problems

Discrete Fourier Transform

Section 5.5. Infinite Series: The Ratio Test

Chem 442 Review for Exam 1. Hamiltonian operator in 1 dimension: ˆ d first term is kinetic energy, second term is potential energy

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

x x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,

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

19 Fourier Series and Practical Harmonic Analysis

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

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

MATH 6101 Fall 2008 Newton and Differential Equations

Name of the Student:

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

GENERATING FUNCTIONS

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

Zeros of Polynomials

2. Fourier Series, Fourier Integrals and Fourier Transforms

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

Time-Domain Representations of LTI Systems

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

ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS

ON THE EIGENFUNCTION EXPANSION METHOD FOR THE CALCULATION OF GREEN S FUNCTIONS

The Random Walk For Dummies

1. Hydrogen Atom: 3p State

*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully

A second look at separation of variables

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

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.

Lecture Notes for Math 251: ODE and PDE. Lecture 32: 10.2 Fourier Series

modes shapes can be obtained by imposing the non-trivial solution condition on the

Prof. Dr. I. Nasser Phys 571, T131 5-Nov-13 Green function.doc. 2 ρ(r)

Solutions to Homework 1

( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.

Math 113, Calculus II Winter 2007 Final Exam Solutions

Numerical Methods in Fourier Series Applications

Castiel, Supernatural, Season 6, Episode 18

Existence of Nonoscillatory Solution of High Order Linear Neutral Delay Difference Equation

8. Applications To Linear Differential Equations

The Riemann Zeta Function

Math 124B January 17, 2012

radians A function f ( x ) is called periodic if it is defined for all real x and if there is some positive number P such that:

11. FINITE FIELDS. Example 1: The following tables define addition and multiplication for a field of order 4.

Math 113 Exam 3 Practice

Chapter 8. Euler s Gamma function

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

C. Complex Numbers. x 6x + 2 = 0. This equation was known to have three real roots, given by simple combinations of the expressions

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

MATH2007* Partial Answers to Review Exercises Fall 2004

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

MAT1026 Calculus II Basic Convergence Tests for Series

Sequences and Series of Functions

: Transforms and Partial Differential Equations

Notes 18 Green s Functions

Solutions for Math 411 Assignment #2 1

Frequency Response of FIR Filters

Seunghee Ye Ma 8: Week 5 Oct 28

14 Separation of Variables Method

Complex Numbers Solutions

Expected Number of Level Crossings of Legendre Polynomials

Created by T. Madas SERIES. Created by T. Madas

4 Separation of Variables

Quantum Mechanics Lecture Notes 10 April 2007 Meg Noah

f(x)g(x) dx is an inner product on D.

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

CHAPTER 5. Theory and Solution Using Matrix Techniques

Fourier Series and their Applications

1 Boas, p. 606, problem

PAPER : IIT-JAM 2010

Radiative Transfer Models and their Adjoints. Paul van Delst

Physics 116A Solutions to Homework Set #9 Winter 2012

n 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent.

11.6 Absolute Convrg. (Ratio & Root Tests) & 11.7 Strategy for Testing Series

Some polynomials defined by generating functions and differential equations

Appendix: The Laplace Transform

Math 210A Homework 1

ENGI Series Page 5-01

Fourier Series and the Wave Equation

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

5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)

LECTURE 21. DISCUSSION OF MIDTERM EXAM. θ [0, 2π). f(θ) = π θ 2

Math 220B - Summer 2003 Homework 1 Solutions

An Extension of Panjer s Recursion

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.

where c is a scaling constant, 0, 0,. r c sinh cos csinh cos cos, csinh cos sin, ccosh sin U csinh sin sin, csinh sin cos,0

Transcription:

Strauss PDEs e: Sectio 4.1 - Exercise 4 Page 1 of 5 Exercise 4 Cosider waves i a resistat medium that satisfy the probem u tt = c u xx ru t for < x < u = at both eds ux, ) = φx) u t x, ) = ψx), where r is a costat, < r < πc/. Write dow the series expasio of the soutio. Soutio The PDE ad its boudary coditios, u, t) = u, t) =, are iear ad homogeeous, so the method of separatio of variabes ca be appied to sove it. Assume a product soutio of the form, ux, t) = Xx)T t), ad pug it ito the PDE ad the boudary coditios. Separate variabes ow. u tt = c u xx ru t XT = c T rxt u, t) = X)T t) = X) = u, t) = X)T t) = X) = XT + rxt = c T T + rt = X X Note that c is a costat ad ca go o either side. The fia aswer wi be the same regardess. We have a fuctio of t o the eft side ad a fuctio of x o the right side. The oy way both fuctios ca be equa is if they are equa to a costat. T + rt = X X = k Vaues of k for which X) = ad X) = are satisfied are caed the eigevaues, ad the otrivia fuctios Xx) associated with them are caed the eigefuctios. Determiatio of Positive Eigevaues: k = µ Assumig k is positive, the differetia equatio for X becomes X = µ = µ X

Strauss PDEs e: Sectio 4.1 - Exercise 4 Page of 5 The geera soutio ca be writte i terms of hyperboic sie ad hyperboic cosie. Xx) = C 1 cosh µx + C sih µx Now use the boudary coditios to determie C 1 ad C. X) = C 1 = X) = C 1 cosh µ + C sih µ = We see that C 1 = ad C =. Hece, oy the trivia soutio Xx) = resuts from cosiderig positive vaues for k, ad there are o positive eigevaues. Determiatio of the Zero Eigevaue: k = Assumig k is zero, the differetia equatio for X becomes X =. = The geera soutio is a iear fuctio. Xx) = C 3 x + C 4 Now use the boudary coditios to determie C 3 ad C 4. X) = C 4 = X) = C 3 + C 4 = We see that C 3 = ad C 4 =. Hece, oy the trivia soutio Xx) = resuts from cosiderig k =, ad zero is ot a eigevaue. Determiatio of Negative Eigevaues: k = λ Assumig k is egative, the differetia equatio for X becomes X = λ. = λ X The geera soutio ca be writte i terms of sie ad cosie. Xx) = C 5 cos λx + C 6 si λx Now use the boudary coditios to determie C 5 ad C 6. X) = C 5 = X) = C 5 cos λ + C 6 si λ =

Strauss PDEs e: Sectio 4.1 - Exercise 4 Page 3 of 5 The secod equatio simpifies to C 6 si λ =. To avoid gettig the trivia soutio, we isist that C 6. Doig so yieds a equatio for the eigevaues. Sove for λ. si λ = λ = π, = 1,,... So the λ = λ = π, = 1,,.... The eigefuctios associated with these eigevaues are Xx) = C 6 si λx X x) = si λ x, = 1,,.... Now sove the differetia equatio for T t). Mutipy both sides by. Brig c λ T to the eft side. T + rt = λ T + rt = c λ T T + rt + c λ T = This is a ODE with costat coefficiets, so its soutio is of the form Substitute this ito the ODE to determie s. Divide both sides by e st. T = e st. s e st + rse st + c λ e st = s + rs + c λ = This is a quadratic equatio for s, so use the quadratic formua to sove for it. s = r ± r 4c λ Now substitute λ = π/. s = r ± r 4c ) π = r ± r ) πc Sice < r < πc/, the quatity uder the square root is egative for every vaue that takes. Factor out 1 ad brig it out of the square root as i. s = r ± i πc ) r = r ± i π c r 4 = r 4 ± i π c r

Strauss PDEs e: Sectio 4.1 - Exercise 4 Page 4 of 5 Thus, the geera soutio to the ODE for T is ) ) T t) = C 7 e r 4 t cos π c r t + C 8 e r 4 t si π c r t. Accordig to the pricipe of iear superpositio, the soutio to the PDE for ux, t) is a iear combiatio of a products T t)x x) over a the eigevaues. [ ) ) ux, t) = A e r 4 t cos π c r t + B e r 4 t si π c r t si πx Therefore, ux, t) = =1 ) ) e r 4 [A t cos π c r 4 t + B si π c r t si πx. =1 The fia task is to use Fourier s method to express the coefficiets, A ad B, i terms of the provided iitia data, φx) ad ψx). ux, ) = =1 A si πx = φx) Mutipy both sides by si λ m x, where m is a positive iteger. ux, ) = =1 A si πx si mπx = φx) si mπx Itegrate both sides with respect to x over the domai the PDE is defied. ˆ =1 A si πx Brig the itegra iside the sum o the eft. ˆ A =1 Sice ad m are itegers, ˆ si πx si πx si mπx si mπx si mπx ˆ ˆ φx) si mπx φx) si mπx { = m m, as ca be verified with trigoometric idetities. This impies that every term i the ifiite series vaishes except for the = m term. Therefore, the first coefficiet is A = ˆ φx) si πx A = ˆ φx) si πx.

Strauss PDEs e: Sectio 4.1 - Exercise 4 Page 5 of 5 I order to use the secod iitia coditio, differetiate the boxed soutio for u with respect to t. u t x, t) = r ) ) ) e r 4 [A t cos π c r 4 t + B si π c r t si πx =1 [ ) + e r 4 t A π c r 4 si π c r t =1 ) 4 + B π c r 4 cos π c r t si πx Now pug i t =. u t x, ) = =1 =1 r ) A si πx =1 + =1 B 4 π c r Brig the costat i frot of the first sum. r ) A si πx 4 + B π c r The first sum is ux, ), ad ux, ) is equa to φx). r ) 4 φx) + B π c r =1 Brig the first term to the right side. 4 B π c r =1 =1 si πx si πx si πx si πx = ψx) = r φx) + ψx) = ψx) = ψx) Now mutipy both sides by si λ m x, where m is a positive iteger. 4 B π c r si πx si mπx = φx) + ψx) si mπx Itegrate both sides with respect to x over the domai the PDE is defied. ˆ 4 B π c r si πx si mπx ˆ φx) + ψx) si mπx =1 =1 Brig the itegra iside the sum o the eft. 4 B π c r ˆ si πx si mπx Every term i this ifiite series vaishes except for the oe whe = m. Therefore, B 4 π c r = ˆ ˆ 4 B = 4 π c r ˆ φx) + ψx) si mπx φx) + ψx) si πx φx) + ψx) si πx.