Partial Differential Equation using separation of variables

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
1 reply · 3K views
Hendrick
Messages
41
Reaction score
0

Homework Statement


Solve the heat flow problem using the method of separation of variables:


Homework Equations


PDE:[tex]\frac{\partial u}{\partial t}=k\frac{\partial^{2} u}{\partial t^{2}}[/tex]
for 0<x<L, 0<t<[tex]\infty[/tex]

BC's:[tex]\frac{\partial u}{\partial x}(0,t)=0,\frac{\partial u}{\partial x}(L,t)=0[/tex]
for 0<t<[tex]\infty[/tex]

IC's: [tex]u(x,0)=[/tex]
{0, [tex]0<x<L/4[/tex]
{1, [tex]L/4<x<3L/4[/tex]
{0, [tex]3L/4<x<L[/tex]
(Piecewise IC)

The Attempt at a Solution


I have separated the variables, then applied the boundary conditions. I am stuck on applying the initial conditions.

I have come up with a general product solution of [tex]u_{n}=F_{n}cos(\frac{n \pi x}{L}) e^{-k(\frac{n \pi x}{L})^{2}t}[/tex]

Trying to combine all product solutions and match the initial data:
[tex]u(x,0)=f(x)[/tex]
[tex]\sum^{\infty}_{n=1}u_{n}(x,0)=f(x)[/tex]
[tex]\sum^{\infty}_{n=1}F_{n}cos(\frac{n \pi x}{L})=f(x)[/tex]


I don't know how to apply the piecewise initial condition, any help would be appreciated. Thank you
 
Physics news on Phys.org
Once you have
[tex]\sum^{\infty}_{n=1}F_{n}cos(\frac{n \pi x}{L})=f(x)[/tex]
you do it exactly the same way you would if f were not "piecewise". Find the Fourier cosine coefficients by doing the appropriate integrals for C0 and Cn for n> 0. The only difference "piecewise" makes is that instead of integrating a single formula from 0 to L, you integrate using the given formulas from 0 to L/4, L/4 to 3L/4, 3L/4 to L and adding those integrals. (Which I notice now is just integrating from L/4 to 3L/4 since outside that the function is 0. The difference between this and just doing the problem on the interval from L/4 to 3L/4 is that you use the whole interval, of length L in determining the "normalization" constant.)