BOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
|
|
- Octavia Gibbs
- 5 years ago
- Views:
Transcription
1 BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann, Robin; all homogeneous) and standard initital conditions (ICs), one can employ the separation of variables method. Let Then u = T(t)X(x, y, z). T c T or T kt = X X = λ, and the series solutions have the form (for λ n 0): u(x, t) = [ A n cos( λ n ct) + B n sin( ] λ n ct) X n (x) n or u(x, t) = n A n e λ nkt X n (x), respectively. X n are eigenfunctions of the eigenvalue problem: X = λx with one of the three classical BCs. (1) 1
2 Eigenvalue problems of on D Rewrite Green s Second Identity as: ( (u v v u)dv = u v ) n v u ds. () n D For any pair of u, v satisfying the classical BCs, the right hand side of this equation is 0. A boundary condition is called symmetric for the operator on D if it can make the right hand side of Eq.() zero for all pairs of functions u, v that satisfy the boundary condition. The classical BCs are symmetric so that the Laplacian satisfies u v = ( u)v. D D D where D (...) is a shorthand for D (...)dv. In terms of the inner product defined by (f, g) = f(x)g(x), D this can be written as (u, v) = ( u, v). Associated with the inner product, the norm can be defined as f = (f,f). The set C n (D) of real (or complex) functions having continuous n th -order derivatives in the domain D is a vector space
3 over the real (or complex) number field. It is also an inner product space and normed space. Similar to the familiar vector spaces R n (or C n ), there are a number of theorems that describe the properties of the eigenvalues and eigenfunctions of the eigenvalue problem (1) associated with the symmetric (self-adjoint) operator. Theorem: Consider any of the problem (1). All the eigenvalues are real. The eigenfunctions can be chosen to be real valued. The eigenfunctions that correspond to distinct eigenvalues are necessarily orthogonal. All the eigenfunctions can be chosen to be real and orthogonal. use self-adjoint property An eigenvalue has multiplicity m if it has m linearly independent eigenfunctions. The eigenfunctions can be orthogonalized by the Graham-Schmidt orthogonalization method. Theorem: All the eigenvalues are positive in the Dirichlet case. All the eigenvalues are positive or zero in the Neumann case, as well as in the Rubin case u/ n + au = 0 (a is a real number ) provided that a 0. use Green s first identity 3
4 Eigenfunction expansions (general Fourier series) A general Fourier series is a series form by all orthogonal eigenfunctions X n of (1) ordered in increasing eigenvalues (i.e. λ 1 λ λ 3...): f(x) = n A n X n (x) with the expansion coefficients defined by A n = (f,x n )/(X n, X n ). The series n A nx n (x) is said to converge in the mean-square sense to f(x) if N f(x) f n (x) 0 as N. D The quantity n=1 f g = (f g, f g) = [ D f g ] 1/ measures the distance between two functions f and g. It is sometimes called the L metric. Theorem: (Least-Square Approximation) Let {X n } be any orthogonal set of functions. Let f <. Let N be a fixed positive integer. Among all possible choices of N constants c 1, c,..., c N, the choice that minimizes 4
5 N f c n X n n=1 is c 1 = A 1,..., c n = A N, where A n = (f,x n )/ X n (1 n N). Proof: Recall that the magnitude (or modulus) of a complex number a is usually written as a (not a ), and a = a a. Denote the error by E N, so that E N = f c n X n = (f,f) (f, c n X n ) ( c n X n, f) + ( c n X n, c m X m ) = f (c n A n + c n A n ) X n + c n X n = ( f ) A n X n + c n A n X n where the sum for n (or m) is over the range 1 to N. The two terms inside the large bracket do not depend on c 1,..., c N. The last term is made 0 by choosing c 1 = A 1,..., c N = A N, the choice that minimizes E N. With the choice c n = A n, N EN = f A n X n 0. n=1 5
6 As this is true for arbitrary N, one gets the inequality: A n X n f. n=1 This is called Bessel s inequality. It is valid as long as f is finite. Theorem: The series A n X n converges to f in the mean-square n=1 sense if and only if A n X n = f. n=1 The equation is called Parseval s equality. Proof: Mean-square convergence means that E N 0 as N. Thus, the equality holds. Definition: The infinite orthogonal set of functions X 1, X,... is called complete if f A n X n = 0 n for all f with f <. Therefore, X 1, X,... is complete if and only if Parseval s equality is true for all f with f <. 6
7 The following theorem is stated without proof. Theorem: (L Convergence & Completeness) The set of eigenfunctions of a symmetric BVP as defined earlier is complete. The eigenfunction expansion converges to f in the mean-square sense in D provided only that f(x) is any function for which f = f(x) is finite. D Note: Due to Parseval s equality, all the Fourier coefficients A n = 0 f = 0. This does not mean that f 0. f can be non-zero at some points. However, we can write f = 0 to label it as belonging to an equivalence class of functions satisfying f = 0. Expansions for special geometric regions In two or three dimensions, it is often more convenient to write the index n as a double or triple index, one for each coordinate, e.g. (l, m, n). The series can be expressed as a double or triple series. Eigenvalues in a rectangular region Consider the domain D = {0 < x < a, 0 < y < b} on the plane. The eigenvalue problem (1) with Dirichlet BC takes the form 7
8 X xx + X yy = λx in D X = 0 on D. With the substitution X = F(x)G(y), the pde can be written as F /F + G /G = λ, and can be split as F = λ x F F(0) = F(a) = 0, G = λ y G G(0) = G(b) = 0, with λ = λ x + λ y. From previous discussions, it is clear that λ x = (mπ/a) (m = 1,,...) and λ y = (nπ/b) (n = 1,,...), so that ( mπ ) ( nπ ) λ = λ mn = +. a b The corresponding eigenfunction is X mn (x, y) = sin(mπx/a) sin(nπy/b). Example: Consider the case a = b = π for the wave equation u tt = c (u xx + u yy ) for 0 < x < a, 0 < y < b with u(x, y, 0) = 0, u t (x, y, 0) = xy(π x)(π y). B nm = 1 c 16 m + n n 3 m 3 π [( 1)n 1][( 1) m 1] Vibration of a drumhead - a disk region If u represents the small displacement of the drum surface from the equilibrium position (u = 0), the problem can be formulated as 8
9 u tt = c u in D = {x + y < a } u = 0 on D u, u t are given functions at t = 0. Following the separation procedure discussed earlier, one can split the temporal and spatial parts u = T(t)X(x). The spatial function X(x) needs to satisfy the eigenvalue problem: X rr + 1 r X r + 1 r X θθ = λx X = 0 on D. Further separation can be made by writing X = R(r)Θ(θ). An ODE for each spatial coordinate is obtained. For θ, the problem is Θ + αθ = 0 with periodic BC. The solution is Θ(θ) = A n cos nθ + B n sin nθ (n = 0, 1,,...), and α = n. For r, the problem is R rr + 1 r R r + (λ n r)r = 0 R(0) finite, R(a) = 0. With the rescaling ρ = λr, the equation can be written as R ρρ + 1 ρ R ρ + (1 n ρ)r = 0. This is Bessel s differential equation of order n. It s solutions need to be found by the power series method. As n is a non-negative integer, the two independent solutions are 9
10 and J n (ρ) = Y n (ρ) = π ( 1) j (ρ/)n+j j!(n + j)!, j=0 ( γ + ln ρ ) J n (ρ) 1 π 1 π n 1 j= n j=0 n j (n j 1)! j!ρ n j ( 1) j (H j + H n+j ) j!(n + j)! ( ρ ) n+j where γ is called Euler s constant and H j = n. J n and Y n are respectively called Bessel function of order n of the first and second kind. So R(ρ) = c 1 J n (ρ) + c Y n (ρ). From the expression for Y n, it is clear that Y n as ρ 0. In order that R be finite at r = 0, the Y n term needs to be dropped. For fixed n, the separated solution is J n ( λr)(a n cos nθ + B n sinnθ). To satisfy the boundary condition at r = a, it is necessary that λ r be a root of J n, namely J n ( λa) = 0. Asymptotically, 10
11 J n (ρ) (ρ πρ cos nπ π ) + O(ρ 3/ ). 4 For each fixed n, J n has an infinite number of roots. The roots ρ nm = λ nm a need to be counted by a different index m, with the arrangement 0 < λ n1 < λ n <... < λ nm <... m = 1,,... The solution for the wave equation can then be written as u(r, θ, t) = (C 0m cos λ 0m ct + D 0m sin λ 0m ct)j 0 ( λ 0m r) m=1 + m,n=1 (C nm cos λ nm ct + D nm sin λ nm ct) J n ( λ nm r)(a nm cos nθ + B nm sinnθ). The initial conditions require that u(r, θ, 0) = φ(r, θ) = C 0m J 0 ( λ 0m r) and + + m,n=1 m=1 C nm J n ( λ nm r)(a nm cos nθ + B nm sin nθ) u t (r, θ, 0) = ψ(r, θ) = m,n=1 D 0m λ0m cj 0 ( λ 0m r) m=1 D nm λnm cj n ( λ nm r)(a nm cos nθ + B nm sinnθ). 11
12 Using the formulas a 0 J n ( λ nm r)j n ( λ np r)rdr = 0 for m p and a [J n ( λ nm r)] rdr = 1 a [J n( λ nm a)], 0 one can find the coefficients through the standard procedure. Note that rdr comes from rdrdθ. Vibrations in a solid ball The Dirichlet eigenvalue problem for a ball with radius a is X rr + r X r + 1 [ 1 r sin θ (sin θx θ) θ + 1 ] sin θ X φφ = λx X(a, θ, φ) = 0. Separating the angular coordinates from the radial coordinate X = Y (θ, φ) R(r), one gets the split problems 1 sin θ (sin θy θ) θ + 1 sin θ Y φφ + αy = 0 and R rr + r R r + (λ α ) R = 0 R(0) finite, R(a) = 0. r The Y equation can be further split with Y (θ, φ) = p(θ)q(φ). 1
13 The equation for q is q φφ + βq = 0, The eigenfunctions are q(φ) = A cosmφ + B sin mφ so that β = m (m = 0, 1,,...). with periodic BC on [ π, π]. The p equation is ( 1 d sin θ d ) ) sin θdθ dθ p + (α m sin p = 0 θ with the conditions p finite at θ = 0, π. With the substitution χ = cos θ, the equation can be cast in the from ( d (1 χ ) d ) ) dχ dχ p + (α m p = 0 1 χ with p finite at χ = ±1. The singular behavior of the equation at χ = ±1 is due to the degeneracy of the coordinate system at the poles and is the reason for the special BC. This equation is known as the associated Legendre equation. It also has to be handled by the power series method. The eigenvalues are α = l(l + 1) where l is an integer m, and the eigenfunctions are of the form P m l (χ) = ( 1)m l l! (1 χ m/ dl+m ) 1) l, dχ l+m(χ 13
14 called the associated Legendre functions. Example: P 0 0 (x) = 1 P 0 1 (x) = x, P 1 1 (x) = (1 x ) 1/ P 0 (x) = 1 (3x 1), P 1 (x) = 3x(1 x ) 1/, P (x) = 3(1 x ) These functions satisfy the orthogonal relationship 1 1 P m l (χ)p m l (χ) dχ = 0 if l l, and 1 [Pl m (χ)] dχ = 1 (l + m)! (l + 1) (l m)!. Combining the θ and φ eigenfunctions, one gets the spherical harmonics Y m l (θ, φ) = P m (cos θ) e imφ l where l = 0, 1,,... and m = l,..., 0,..., l, The sines and cosines have been replaced by the complex exponentials. These functions satisfy the orthogonal relationship π π 0 0 Y m l (θ, φ)y m l (θ, φ) sinθdθdφ = δ ll δ mm 4π (l + m)! l + 1(l m)!. Given α = l(l + 1), the R equation can be written as R rr + ( ) r R l(l + 1) r + λ R = 0. r 14
15 It is similar to the Bessel equation. The difference is in the factor r (instead of 1 r ) in front of the derivative R r. It is possible to convert the equation to the Bessel form with a substitution w(r) = rr(r), i.e. R(r) = r 1/ w(r) and a rescaling ρ = λ r. The result is w ρρ + 1 ) ρ w ρ + (1 s w = 0. ρ where s = l(l + 1) = l + 1. As s is not an integer, the two independent solutions are ( 1) j ( ρ ) j±s J ±s (ρ) = Γ(j + 1)Γ(j + s + 1) where Γ(s) = j=0 0 e t t s 1 dt is the Γ function which has the property Γ(s + 1) = sγ(s). If s is an integer Γ(s + 1) = s!. Since J (l+ 1 ) is singular at ρ = 0, only the solution J l+ 1 should be retained. As ρ, J s again has the asymptotic formula J s (ρ) (ρ πρ cos sπ π ) + O(ρ 3/ ). 4 This formula turns out to be exact when s = ±1/ so that J 1/ = πρ sin ρ and J 1/ = cos ρ. πρ The asymptotic formula shows that J s has infinitely many roots. If the index j (= 1,,...) is used to count the roots 15
16 ρ lj = λ lj a of J l+ 1, the eigenfunctions associated with the eigenvalue λ lj can be expressed as X lmj (r,θ, φ) = J l+ 1 ( λ lj r) r Y m l (θ, φ) l m l. Therefore, the eigenvalue λ lj has multiplicity l + 1. The series solution of the wave equation has the form u(r, θ, φ, t) l = l=0 m= l (A lmj cos λ lj ct + B lmj sin λ lj ct) X lmj. j=1 Note that Y m l = Yl m. To ensure that u is real, the coefficients satisfy A l( m)j = A lmj (same for the B coefficients). Example: Consider a firecracker set off at the center of a spherical enclosure with hard walls. What kind of sound waves will be generated? Suppose that this can be model by the following mathematical problem u tt = c u in D = {x + y + z < a } { 1 for r < ɛ u(x, 0) = 0 u t (x, 0) = δ ɛ (r) = 0 for r > ɛ. Steps: 16
17 The series solution is of the general form u = (A lmj cos β lj ct + B lmj sinβ lj ct)x lmj lmj where β lj = λ lj and X lmj = 1 J l+1/ (β lj r) P m r l (cos θ)e imφ. First, note that the initial conditions do not depend on θ and φ. So the solution has the same property (the expansion coefficients are determined by the initial conditions). This means that the only non-zero terms are those with l = m = 0. The sum is therefore only over j = 1,... For this reason, only the j index will be displayed. Second, u(x, 0) = 0 = j A jx j A j = 0. It is only necessary to determine B j from the condition δ ɛ u t (x, 0) = j B jβ j cx j, as B j = (δ ɛ, X j ) β j c X j = a 0 β j c δ ɛ (r)(1/ r)j 1/ (β j r)r dr a 0 (1/r)J 1/ (β jr)r dr. J 1/ has a very nice connection with the sine function, J 1/ (β j r) = πβ j r sin β jr. As β j a has to be a root of J 1/ (β j r), it is clear from the sine function that β j a = jπ with j = 1,,... The 17 (3)
18 integral in the numerator of Eq.(3) becomes ɛ sin(β j r)rdr = 0 πβ j 1 πβ j βj [sin(β j ɛ) (β j ɛ) cos(β j ɛ)]. The integral in the denominator is a sin (β j r)dr = β j a 0 πβ j πβj = a jπ. Therefore, πβ j B j = and the solution is a c(jπ) 3[sin(β jɛ) (β j ɛ) cos(β j ɛ)], u = a cπ 3 [sin(β j ɛ) (β j ɛ) cos(β j ɛ)] rj 3 j=1 ( ) ( jπct jπr sin sin a a ). Problem Set 11 18
19 Justification for separation of variables Given that the set of eigenfunctions {v m } of in an interval (a, b) is complete, one can prove the following. Theorem: Let D = D 1 D be a rectangle in the xy plane. The set of products {v m (x)w n (y) : v m is an eigenfunction in D 1, w n is an eigenfunction in D, m, n = 1,,...} is a complete set of eigenfunctions for with the given boundary conditions. Proof: Note that each product is an eigenfunction since (v m w n ) = (α m + β n )v m w n. It is also clear that the eigenfunctions of are mutually orthogonal. Suppose that there were an eigenfunction u(x, y) of D other than these products. Then, for some λ, u = λu. If λ were different from every one of the α m + β n, then u is orthogonal to all the products v m w n. So 0 = (u, v m w n ) = [ u(x, y)v m (x)dx]w n (y)dy. By the completeness of w n, u(x, y)v m (x) = 0 y. By the completeness of v m, u(x, y) = 0 (x, y). It is not a valid eigenfunction. 19
20 If λ = α m + β n for some m, n (this could be true for one or several such pairs), consider the difference φ = u c mn v m w n where the sum is over all the m, n pairs for which λ = α m + β n and c mn = (u, v m w n )/ v m w n. As φ is then orthogonal to all the v m w n, it is again 0 according to the earlier agrument. So, u is no more than a linear combination of the old products. 0
LAPLACE EQUATION. = 2 is call the Laplacian or del-square operator. In two dimensions, the expression of in rectangular and polar coordinates are
LAPLACE EQUATION If a diffusion or wave problem is stationary (time independent), the pde reduces to the Laplace equation u = u =, an archetype of second order elliptic pde. = 2 is call the Laplacian or
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
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 informationd Wave Equation. Rectangular membrane.
1 ecture1 1.1 2-d Wave Equation. Rectangular membrane. The first problem is for the wave equation on a rectangular domain. You can interpret this as a problem for determining the displacement of a flexible
More informationSAMPLE FINAL EXAM SOLUTIONS
LAST (family) NAME: FIRST (given) NAME: ID # : MATHEMATICS 3FF3 McMaster University Final Examination Day Class Duration of Examination: 3 hours Dr. J.-P. Gabardo THIS EXAMINATION PAPER INCLUDES 22 PAGES
More informationMath 2930 Worksheet Final Exam Review
Math 293 Worksheet Final Exam Review Week 14 November 3th, 217 Question 1. (* Solve the initial value problem y y = 2xe x, y( = 1 Question 2. (* Consider the differential equation: y = y y 3. (a Find the
More informationSeparation of Variables
Separation of Variables A typical starting point to study differential equations is to guess solutions of a certain form. Since we will deal with linear PDEs, the superposition principle will allow us
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More informationQuantum Mechanics in 3-Dimensions
Quantum Mechanics in 3-Dimensions Pavithran S Iyer, 2nd yr BSc Physics, Chennai Mathematical Institute Email: pavithra@cmi.ac.in August 28 th, 2009 1 Schrodinger equation in Spherical Coordinates 1.1 Transforming
More informationIn what follows, we examine the two-dimensional wave equation, since it leads to some interesting and quite visualizable solutions.
ecture 22 igher-dimensional PDEs Relevant section of text: Chapter 7 We now examine some PDEs in higher dimensions, i.e., R 2 and R 3. In general, the heat and wave equations in higher dimensions are given
More informationMath 316/202: Solutions to Assignment 7
Math 316/22: Solutions to Assignment 7 1.8.6(a) Using separation of variables, we write u(r, θ) = R(r)Θ(θ), where Θ() = Θ(π) =. The Laplace equation in polar coordinates (equation 19) becomes R Θ + 1 r
More informationMATH 241 Practice Second Midterm Exam - Fall 2012
MATH 41 Practice Second Midterm Exam - Fall 1 1. Let f(x = { 1 x for x 1 for 1 x (a Compute the Fourier sine series of f(x. The Fourier sine series is b n sin where b n = f(x sin dx = 1 = (1 x cos = 4
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationThe Laplacian in Polar Coordinates
The Laplacian in Polar Coordinates R. C. Trinity University Partial Differential Equations March 17, 15 To solve boundary value problems on circular regions, it is convenient to switch from rectangular
More informationChapter 5.8: Bessel s equation
Chapter 5.8: Bessel s equation Bessel s equation of order ν is: x 2 y + xy + (x 2 ν 2 )y = 0. It has a regular singular point at x = 0. When ν = 0,, 2,..., this equation comes up when separating variables
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 informationTHE METHOD OF SEPARATION OF VARIABLES
THE METHOD OF SEPARATION OF VARIABES To solve the BVPs that we have encountered so far, we will use separation of variables on the homogeneous part of the BVP. This separation of variables leads to problems
More informationSpherical Harmonics on S 2
7 August 00 1 Spherical Harmonics on 1 The Laplace-Beltrami Operator In what follows, we describe points on using the parametrization x = cos ϕ sin θ, y = sin ϕ sin θ, z = cos θ, where θ is the colatitude
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationDiffusion on the half-line. The Dirichlet problem
Diffusion on the half-line The Dirichlet problem Consider the initial boundary value problem (IBVP) on the half line (, ): v t kv xx = v(x, ) = φ(x) v(, t) =. The solution will be obtained by the reflection
More informationIntroduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series
CHAPTER 5 Introduction to Sturm-Liouville Theory and the Theory of Generalized Fourier Series We start with some introductory examples. 5.. Cauchy s equation The homogeneous Euler-Cauchy equation (Leonhard
More informationWaves on 2 and 3 dimensional domains
Chapter 14 Waves on 2 and 3 dimensional domains We now turn to the studying the initial boundary value problem for the wave equation in two and three dimensions. In this chapter we focus on the situation
More informationStriking a Beat. Ashley Martin PHY 495. Spring Striking a Beat. Ashley Martin PHY 495. Introduction. Outline. Cartesian Coordinates
Spring 2012 Polar Where is it optimal to strike a circular drum? Polar Daniel Bernoulli (1700-1782) - introduced concept of Bessel functions Leonhard Euler (1707-1783) - used Bessel funtions of both zero
More informationDUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs
DUHAMEL S PRINCIPLE FOR THE WAVE EQUATION HEAT EQUATION WITH EXPONENTIAL GROWTH or DECAY COOLING OF A SPHERE DIFFUSION IN A DISK SUMMARY of PDEs MATH 4354 Fall 2005 December 5, 2005 1 Duhamel s Principle
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationHomework for Math , Fall 2016
Homework for Math 5440 1, Fall 2016 A. Treibergs, Instructor November 22, 2016 Our text is by Walter A. Strauss, Introduction to Partial Differential Equations 2nd ed., Wiley, 2007. Please read the relevant
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More informationc2 2 x2. (1) t = c2 2 u, (2) 2 = 2 x x 2, (3)
ecture 13 The wave equation - final comments Sections 4.2-4.6 of text by Haberman u(x,t), In the previous lecture, we studied the so-called wave equation in one-dimension, i.e., for a function It was derived
More information# Points Score Total 100
Name: PennID: Math 241 Make-Up Final Exam January 19, 2016 Instructions: Turn off and put away your cell phone. Please write your Name and PennID on the top of this page. Please sign and date the pledge
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationPartial Differential Equations
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
More informationDepartment of Mathematics
INDIAN INSTITUTE OF TECHNOLOGY, BOMBAY Department of Mathematics MA 04 - Complex Analysis & PDE s Solutions to Tutorial No.13 Q. 1 (T) Assuming that term-wise differentiation is permissible, show that
More informationEigenfunctions on the surface of a sphere. In spherical coordinates, the Laplacian is. u = u rr + 2 r u r + 1 r 2. sin 2 (θ) + 1
Eigenfunctions on the surface of a sphere In spherical coordinates, the Laplacian is u = u rr + 2 r u r + 1 r 2 [ uφφ sin 2 (θ) + 1 sin θ (sin θ u θ) θ ]. Eigenfunctions on the surface of a sphere In spherical
More informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationBoundary-value Problems in Rectangular Coordinates
Boundary-value Problems in Rectangular Coordinates 2009 Outline Separation of Variables: Heat Equation on a Slab Separation of Variables: Vibrating String Separation of Variables: Laplace Equation Review
More informationBessel Functions. A Touch of Magic. Fayez Karoji 1 Casey Tsai 1 Rachel Weyrens 2. SMILE REU Summer Louisiana State University
Bessel s Function A Touch of Magic Fayez Karoji 1 Casey Tsai 1 Rachel Weyrens 2 1 Department of Mathematics Louisiana State University 2 Department of Mathematics University of Arkansas SMILE REU Summer
More information4.10 Dirichlet problem in the circle and the Poisson kernel
220 CHAPTER 4. FOURIER SERIES AND PDES 4.10 Dirichlet problem in the circle and the Poisson kernel Note: 2 lectures,, 9.7 in [EP], 10.8 in [BD] 4.10.1 Laplace in polar coordinates Perhaps a more natural
More informationPhysics 116C Solutions to the Practice Final Exam Fall The method of Frobenius fails because it works only for Fuchsian equations, of the type
Physics 6C Solutions to the Practice Final Exam Fall 2 Consider the differential equation x y +x(x+)y y = () (a) Explain why the method of Frobenius method fails for this problem. The method of Frobenius
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationMath Assignment 14
Math 2280 - Assignment 14 Dylan Zwick Spring 2014 Section 9.5-1, 3, 5, 7, 9 Section 9.6-1, 3, 5, 7, 14 Section 9.7-1, 2, 3, 4 1 Section 9.5 - Heat Conduction and Separation of Variables 9.5.1 - Solve the
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationMidterm 2: Sample solutions Math 118A, Fall 2013
Midterm 2: Sample solutions Math 118A, Fall 213 1. Find all separated solutions u(r,t = F(rG(t of the radially symmetric heat equation u t = k ( r u. r r r Solve for G(t explicitly. Write down an ODE for
More informationPartial Differential Equations
M3M3 Partial Differential Equations Solutions to problem sheet 3/4 1* (i) Show that the second order linear differential operators L and M, defined in some domain Ω R n, and given by Mφ = Lφ = j=1 j=1
More informationLaplace s equation in polar coordinates. Boundary value problem for disk: u = u rr + u r r. r 2
Laplace s equation in polar coordinates Boundary value problem for disk: u = u rr + u r r + u θθ = 0, u(a, θ) = h(θ). r 2 Laplace s equation in polar coordinates Boundary value problem for disk: u = u
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationMATH 220: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY
MATH 22: INNER PRODUCT SPACES, SYMMETRIC OPERATORS, ORTHOGONALITY When discussing separation of variables, we noted that at the last step we need to express the inhomogeneous initial or boundary data as
More informationLecture Notes for MAE 3100: Introduction to Applied Mathematics
ecture Notes for MAE 31: Introduction to Applied Mathematics Richard H. Rand Cornell University Ithaca NY 14853 rhr2@cornell.edu http://audiophile.tam.cornell.edu/randdocs/ version 17 Copyright 215 by
More informationCircular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431
The Problem Polar coordinates Solving the Problem by Separation of Variables Circular Membranes Farlow, Lesson 30 Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 November
More informationM412 Assignment 5 Solutions
M41 Assignment 5 Solutions 1. [1 pts] Haberman.5.1 (a). Solution. Setting u(x, y) =X(x)Y(y), we find that X and Y satisfy X + λx = Y λy =, X () = X () = Y() =. In this case, we find from the X equation
More informationPartial Differential Equations
Partial Differential Equations Xu Chen Assistant Professor United Technologies Engineering Build, Rm. 382 Department of Mechanical Engineering University of Connecticut xchen@engr.uconn.edu Contents 1
More informationThe Wave Equation on a Disk
R. C. Trinity University Partial Differential Equations April 3, 214 The vibrating circular membrane Goal: Model the motion of an elastic membrane stretched over a circular frame of radius a. Set-up: Center
More informationFinal Examination Solutions
Math. 42 Fulling) 4 December 25 Final Examination Solutions Calculators may be used for simple arithmetic operations only! Laplacian operator in polar coordinates: Some possibly useful information 2 u
More informationDifferential Equations
Differential Equations Problem Sheet 1 3 rd November 2011 First-Order Ordinary Differential Equations 1. Find the general solutions of the following separable differential equations. Which equations are
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 informationGREEN S IDENTITIES AND GREEN S FUNCTIONS
GREEN S IENTITIES AN GREEN S FUNCTIONS Green s first identity First, recall the following theorem. Theorem: (ivergence Theorem) Let be a bounded solid region with a piecewise C 1 boundary surface. Let
More informationConsequences of Orthogonality
Consequences of Orthogonality Philippe B. Laval KSU Today Philippe B. Laval (KSU) Consequences of Orthogonality Today 1 / 23 Introduction The three kind of examples we did above involved Dirichlet, Neumann
More informationBranch: Name of the Student: Unit I (Fourier Series) Fourier Series in the interval (0,2 l) Engineering Mathematics Material SUBJECT NAME
13 SUBJECT NAME SUBJECT CODE MATERIAL NAME MATERIAL CODE UPDATED ON : Transforms and Partial Differential Equation : MA11 : University Questions :SKMA13 : May June 13 Name of the Student: Branch: Unit
More informationSolving the Heat Equation (Sect. 10.5).
Solving the Heat Equation Sect. 1.5. Review: The Stationary Heat Equation. The Heat Equation. The Initial-Boundary Value Problem. The separation of variables method. An example of separation of variables.
More information1. Partial differential equations. Chapter 12: Partial Differential Equations. Examples. 2. The one-dimensional wave equation
1. Partial differential equations Definitions Examples A partial differential equation PDE is an equation giving a relation between a function of two or more variables u and its partial derivatives. The
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationFINAL EXAM, MATH 353 SUMMER I 2015
FINAL EXAM, MATH 353 SUMMER I 25 9:am-2:pm, Thursday, June 25 I have neither given nor received any unauthorized help on this exam and I have conducted myself within the guidelines of the Duke Community
More informationProblem (p.613) Determine all solutions, if any, to the boundary value problem. y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6,
Problem 10.2.4 (p.613) Determine all solutions, if any, to the boundary value problem y + 9y = 0; 0 < x < π, y(0) = 0, y (π) = 6, by first finding a general solution to the differential equation. Solution.
More informationMATH FALL 2014 HOMEWORK 10 SOLUTIONS
Problem 1. MATH 241-2 FA 214 HOMEWORK 1 SOUTIONS Note that u E (x) 1 ( x) is an equilibrium distribution for the homogeneous pde that satisfies the given boundary conditions. We therefore want to find
More information3 Green s functions in 2 and 3D
William J. Parnell: MT34032. Section 3: Green s functions in 2 and 3 57 3 Green s functions in 2 and 3 Unlike the one dimensional case where Green s functions can be found explicitly for a number of different
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationSolutions VI. MAE294A/SIO203A: Methods in Applied Mechanics Fall Quarter
MAE94A/SIO3A: Methods in Applied Mechanics Fall Quarter 17 http://web.eng.ucsd.edu/~sgls/mae94a_17/ Solutions VI 1 We let Tx, t) = Xx)τt). We have two homogeneous conditions in x, so let X = λ X to obtain
More informationProcess of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion
Brian England Objective Process of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion Cylindrical Coordinates Initial and Boundary Conditions Methodologies and
More informationAn Introduction to Bessel Functions
An Introduction to R. C. Trinity University Partial Differential Equations Lecture 17 Bessel s equation Given p 0, the ordinary differential equation x 2 y + xy + (x 2 p 2 )y = 0, x > 0 is known as Bessel
More informationFinal Examination Linear Partial Differential Equations. Matthew J. Hancock. Feb. 3, 2006
Final Examination 8.303 Linear Partial ifferential Equations Matthew J. Hancock Feb. 3, 006 Total points: 00 Rules [requires student signature!]. I will use only pencils, pens, erasers, and straight edges
More informationHomework 7 Math 309 Spring 2016
Homework 7 Math 309 Spring 2016 Due May 27th Name: Solution: KEY: Do not distribute! Directions: No late homework will be accepted. The homework can be turned in during class or in the math lounge in Pedelford
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More information24 Solving planar heat and wave equations in polar coordinates
24 Solving planar heat and wave equations in polar coordinates Now that all the preparations are done, I can return to solving the planar heat and wave equations in domains with rotational symmetry. 24.1
More informationAngular momentum. Quantum mechanics. Orbital angular momentum
Angular momentum 1 Orbital angular momentum Consider a particle described by the Cartesian coordinates (x, y, z r and their conjugate momenta (p x, p y, p z p. The classical definition of the orbital angular
More informationMath 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
Name Section Math 51 December 14, 5 Answer Key to Final Exam There are 1 questions on this exam. Many of them have multiple parts. The point value of each question is indicated either at the beginning
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 informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More information14 Fourier analysis. Read: Boas Ch. 7.
14 Fourier analysis Read: Boas Ch. 7. 14.1 Function spaces A function can be thought of as an element of a kind of vector space. After all, a function f(x) is merely a set of numbers, one for each point
More informationENGI 9420 Lecture Notes 8 - PDEs Page 8.01
ENGI 940 Lecture Notes 8 - PDEs Page 8.01 8. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationCHAPTER 4. Introduction to the. Heat Conduction Model
A SERIES OF CLASS NOTES FOR 005-006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 4 A COLLECTION OF HANDOUTS ON PARTIAL DIFFERENTIAL EQUATIONS
More informationMA 201: Method of Separation of Variables Finite Vibrating String Problem Lecture - 11 MA201(2016): PDE
MA 201: Method of Separation of Variables Finite Vibrating String Problem ecture - 11 IBVP for Vibrating string with no external forces We consider the problem in a computational domain (x,t) [0,] [0,
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
More informationProblem 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.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationMATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).
MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). Bessel s differential equation of order m 0: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0 The equation is considered on the
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
More informationSession 5 Heat Conduction in Cylindrical and Spherical Coordinates I
Session 5 Heat Conduction in Cylindrical and Spherical Coordinates I 1 Introduction The method of separation of variables is also useful in the determination of solutions to heat conduction problems in
More informationswapneel/207
Partial differential equations Swapneel Mahajan www.math.iitb.ac.in/ swapneel/207 1 1 Power series For a real number x 0 and a sequence (a n ) of real numbers, consider the expression a n (x x 0 ) n =
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More information