AMAT416/516 Partial Differential Equations
MWF 9:20-10:15 ES 146
Instructor: Professor Edward Thomas, ES 132F
Email : firstname.lastname@example.org
Office hours: TBA
Text: “Applied PDEs with Fourier Series, etc”, Fifth Edition, by Richard Haberman
2) a) compute Δ( x2 y z ) (I can’t get a “del” symbol to print so I had to use a delta instead…sorry) b) compute Δ2 ( x + y + z ) c) find several
solutions of Δ2u=0. d) Grad students: If there is only one variable,
Δ2 u=0 reduces to an ODE. What is the general solution of this equation?
3) Solve the steady state heat equation for a thin rod of length L=42 ( inches) with u(0) =100, u(42) = 14 ( temps in degrees). What is the temperature at x = 36.8 ? at x=.03? at x=3 ?Hint: when x =15, u (x) is approximately equal to 69.3 degrees.
4) page 51 2.3.1 (b), (c) and (f) and 2.3.2 (a) and (b) ( with lambda >0 )
5) page 52 # 2.3.5 ( orthogonality of the sines) plus, for you 516 students, compute the general Fourier coefficient, Bm , for F(x) = x.
Now… in the Senility Department…. when I was verifying that the Fourier series actually adds up to the initial temperature f(x) = 100, I forgot to include the sine terms in the infinite series ( whata dope!). As you will see next class, this makes the series add up to 100 exactly.
6) First, in the series expansion 100=400/ Pi( 1 – 1/3 + 1/5…), compute the RHS going out to 7 terms and then 8 terms.
Second, take the full solution for u( x,t). take L = Pi and take x = Pi/2. Using the first two terms of the series, compute the temperature when t = 1/10 and then when t = ½. ( the temp should drop considerably in that time)
By the way, back about two lectures ago, we started taking k=1 for simplicity. Thanks to Cody for pointing this out.
7) Take k = 1.Solve the heat equation with insulated boundaries and initial temperature f(x) = ½ for 0< x <L/2 and 1 for L/2 < x < L. Expand f(x) in a cosine series, write out the complete solution, and say what happens as x approaches infinity ( steady state)OOPS I should have said “as t approaches infinity”.
8) A handout on expanding f(x) , -L<x<+L, in a Fourier sine /cosine series.I’ll leave a couple on my door.
9) Another handout on shifting boundary conditions around. I’ll leave some copies on my door. Added Thursday AM: Couple of remarks about this problem. First, it might help to rephrase things like this. Given a function u(x,y) that satisfies Laplace’s equation, define a new function v(x,y) = u( L-x,y). Then v also satisfies Laplace’s equation. ( When we say a function satisfies Laplace’s equation, that means it satisfies the equation for all values of the two independent variables.) The second question is simply to figure out what the boundary condition on the function u becomes for the function v.
10) We saw that the temperature at the midpoint of our plate does not vary linearly across the plate. Taking L=1 and using the first two terms of the series, estimate the temp at the midpoint when H=1, then when H=5. Grad students, analyze the temperature as a function of H as H goes to infinity.
On another tack, if you fix a value of H, someone might want to plot the temperature across the plate at height H/2 as a function of x…that would be cool.
11) Fix L=H=1 and estimate values of u(x,1/2) as x=.5,.6,.7, etc. approaching x=1. ( You could use the average of the first two values of the series as the estimator.) Then, in the polar form of Laplace’s equation in a disk, substitute u(r,theta)= G(r)Phi(theta) and separate variables as discussed in class.
[ I’m sorry that I can’t get Greek letters to appear when I go to post these notes, so I have to write “theta” , etc.)
12) Grads: Verify that the eigenvalues for the disk are lambda = n2 assuming perfect thermal contact along the seam theta = +/- Pi
Everybody: Verify that the solutions for G(r) are a rn + b r-n if r is not equal to 0, and c + d ln(r) if r = 0. Suppose we are given the boundary function f(theta) =0 for –Pi < theta<-Pi/2, 150o from there to+ Pi/2, 0 from there to +Pi. Calculate the temperature values at the origin and at r=a/2, theta = 0, Pi/2. Pi, 3Pi/2.
13) Compute the Fourier ( sine and cosine) series for the function f(x) = 0 for –L<x<0, 15 for 0<x<L ( I hope I remembered that 15 right. In particular check the jump condition.
Secondly, compute the Fourier series for the even extension of f(x) =sin(Pi*x/L) on the interval –Pi < x < Pi. Then take L=1 and add up the first four nonzero terms at x =1/2.
Grad special) Grad students: compute the integral of sin(nPix/L) times cos( m Pix/L) on the interval [0,L]Oops…I might have said to take L=1 in this problem. If you want to do it with L=1, go ahead.
14) The Wave Equation In “plucked string”, show that the second initial condition implies that all Bn = 0. Then you want to expand f(x) in a sine series, so your next step is to compute the formula for f….. go ahead and do this.
This will set you up to compute An. But, you may assume ( as a gift) that all the even Ans are 0 ( Grads prove this) and the formula for the odd An is what I gave in class.
Actually, let me change that formula:
An = 4/L times the integral from 0 to L/2 times( 2h/L)times x sin(nPi/L)x dx
Finally write out the formula for u(x,t) …you only need to write the first few terms so the pattern is obvious.
( It goes without saying that you didn’t want to miss class on this particular day.)
15) page 142 # 4.4.3( b)
16) Make the change of variables z = sqrt( lambda) times r to obtain Bessel’s equation in the form 7.7.25 on page 299. Be explicit about how the chain rule comes into the computation.
17) Find the first 5 terms of the power series expansion for J0 (z), showing the details of the computation. Grads do this plus the same thing for J1 (z)
18) Draw the nodal patterns for : m =3, n=2, sine; m=1,n=3, cosine; and m=4,n=2, sine
EXAM I Wednesday October 23
19) A variety of easy problems to get comfy with complex variables ( you kinda had to be there) plus 516 students be ready to present a derivation of Euler’s Formula based on Maclaurin series.
20) i) Verify the orthogonality relationship for complex exponentials ( page 129 in the text) ii) compute cm for the step function given in class and iii) grad students be ready to prove that if f(x) is real then, in the complex Fourier series for f, we must have c-n equal to the conjugate of cn