sol9 - CAAM 336 DIFFERENTIAL EQUATIONS Problem Set 9...

Info iconThis preview shows pages 1–4. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: CAAM 336 DIFFERENTIAL EQUATIONS Problem Set 9 Solutions Posted Thursday 28 October 2010. Due Wednesday 3 November 2010, 5pm. 1. [40 points: 12 points each for (a) and (b); 6 points for (c); 10 points for (d)] (a) Consider the function u ( x ) = 1 , x [0 , 1 / 3]; , x (1 / 3 , 2 / 3); 1 , x [2 / 3 , 1] . Recall that the eigenvalues of the operator L : C 2 N [0 , 1] C [0 , 1], Lu =- u 00 are n = n 2 2 for n = 0 , 1 ,... with associated (normalized) eigenfunctions ( x ) = 1 and n ( x ) = 2cos( nx ) , n = 1 , 2 ,.... We wish to write u ( x ) as a series of the form u ( x ) = X n =0 a n (0) n ( x ) , where a n (0) = ( u , n ). Compute these inner products a n (0) = ( u , n ) by hand and simplify as much as possible. For m = 0 , 2 , 4 , 80, plot the partial sums u ,m ( x ) = m X n =0 a n (0) n ( x ) . (You may superimpose these on one single, well-labeled plot if you like.) (b) Write down a series solution to the homogeneous heat equation u t ( x,t ) = u xx ( x,t ) , < x < 1 , t with Neumann boundary conditions u x (0 ,t ) = u x (1 ,t ) = 0 and initial condition u ( x, 0) = u ( x ). Create a plot showing the solution at times t = 0 , . 002 , . 05 , . 1. You will need to truncate your infinite series to show this plot. Discuss how the number of terms you use in this infinite series affects the accuracy of your plots. (c) Describe the behavior of your solution as t . (To do so, write down a formula for the solution in the limit t .) (d) How would you expect the solution to the inhomogeneous heat equation u t ( x,t ) = u xx + 1 , < x < 1 , t with Neumann boundary conditions u x (0 ,t ) = u x (1 ,t ) = 0 to behave as t ? Solution. (a) To expand u ( x ) in the form u ( x ) = X n =0 a n (0) n ( x ) , we must compute the coefficients a n (0). For n = 0 we compute a (0) = Z 1 u ( x ) 1d x = Z 1 / 3 1d x + Z 1 2 / 3 1d x = 2 / 3 . For n > 0 we have a n (0) = 2 Z 1 u ( x )cos( nx ) dx = 2 Z 1 / 3 cos( nx ) dx + Z 1 2 / 3 cos( nx ) dx = 2 h sin( nx ) n i 1 / 3 + h sin( nx ) n i 1 2 / 3 = 2(sin( n/ 3)- sin(2 n/ 3)) n . [GRADERS: this last expression is sufficiently simplified for full credit.] Note that sin(2 n/ 3) = 2sin( n/ 3)cos( n/ 3), and hence sin( n/ 3)- sin(2 n/ 3) = sin( n/ 3)(1- 2cos( n/ 3)) . Thus we have a n (0) = 0 in two cases: if n is a multiple of 3, or if cos( n/ 3) = 1 / 2. The former occurs when n = 3 , 6 , 9 , 12 , 15 ,... , whiles the latter occurs when n/ 3(mod 2 ) = / 3 or 5 / 3, and hence a n (0) = 0 when n = 1 + 6 p for integers p 0 or n =- 1 + 6 p for integers p 1. Together, this implies that for all odd integers n , a n (0) = 0. We end up with the partial sums shown in the following figure. (MATLAB code follows at the end of this solution.) 0.2 0.4 0.6 0.8 1-0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 x m (x) m=0 m=2 m=4 m=80 (b) We seek a series solution of the form u ( x,t...
View Full Document

This note was uploaded on 01/21/2012 for the course CAAM 330 taught by Professor Tompson during the Fall '09 term at UVA.

Page1 / 11

sol9 - CAAM 336 DIFFERENTIAL EQUATIONS Problem Set 9...

This preview shows document pages 1 - 4. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online