Heat conduction in bars: varying the boundary condition. Harmonics problem

[Benzoate]

1. The problem statement, all variables and given/known data

A problem with odd harmonics only. Show that the solution of the heat equation du/dt=c2*(d2u)/(dx2), subject to boundary conditions u(0,t)=0 and ux(L,t)=0, and the initial condition u(x,0)=f(x) , is

u(x,t)= \sum Bnsin[(\pi/2L)(2n+1)x]e-((c*\pi/2L)*(2n+1))^2

where n extends from 0 to \infty

Bn=2/L\int f(x) sin[(\pi/2L)*(2n+1)x]dx

where the limits extend from 0 to L

2. Relevant equations

Seperation of variables

3. The attempt at a solution
u(x,t)=X(x)T(t)
u(x,0)=f(x)
u2(x,t)=0
T'-kc2T=0
x''-kx=0
x'(0)=0
x'(L)=0
u2(0,t)=0

if u(0,t)=0 then u2(0,t)=0

u2(0,t)=0= u2(L,t)

not sure how to finish this derivation.
1. The problem statement, all variables and given/known data



2. Relevant equations



3. The attempt at a solution
1. The problem statement, all variables and given/known data



2. Relevant equations



3. The attempt at a solution

[gabbagabbahey]

u(x,t)=X(x)T(t)
u(x,0)=f(x)
T'-kc2T=0
x''-kx=0


This part is correct. What does that make the general solution for T(t) and X(t) and hence u(x,t)? Remember, you will need to examine 3 different cases: (1)k=0, (2)k<0 and (3)k>0.


u2(x,t)=0


Huh? What is this supposed to mean?


x'(0)=0
x'(L)=0
u2(0,t)=0

if u(0,t)=0 then u2(0,t)=0

u2(0,t)=0= u2(L,t)

Where are you getting these boundary conditions? They are not equal to the conditions given in the question.

[Benzoate]

This part is correct. What does that make the general solution for T(t) and X(t) and hence u(x,t)? Remember, you will need to examine 3 different cases: (1)k=0, (2)k<0 and (3)k>0.

3 cases : you mean when

case 1: lambda =alpha^2 ==> X-alpha^2*X

X''-alphaX=0

r^2-alpha^2=0
(r-alpha)(r+alpha)=0 ==> r= +-alpha

X=c1e^alpha*x + c2e^alphae^-alpha
X'(0)=0 ==> c1*alpha-c2*alpha
X'(L)=alpha(e^alpha*L - e^(-alpha*L))c1=0 ==> c1=0

case 2 lamda =0

r=0

X=c1*x+c2

X'(0)=0 ==> c1=0

X'(L) =0 ==>0=0

X(x)=c2


case 3

lambda = alpha^2 >0

r^2+alpha^2 =0 ==> r=+- alpha*i

X1=cos(alpha*x) X2=sin(alpha*x)

X=X1+X2= cos(alpha*x)+sin(alpha*x)



Huh? What is this supposed to mean?
Ignore all the functions with u2. Those are all mistakes. ; I was trying to say that if u(0,t)=0 then u(0,t)=0==> ux(0,t)=0

Therefore, ux(0,t)=0= ux(L,t)

[klondike]

1. After SOV, you'll get 2 ODEs. One for spatial the other for time, then you select proper eigenvalue such that the PDE has nontrivial and physically possible solution considering the given BCs. In the case of
\dfrac {T'}{T}= k \dfrac {X''}{X}=- \lambda with Dirichlet boundary conditions, the possible eigenvalue is \lambda >0 (easy to find out by plugging the BCs in).

2. Solve the spatial ODE, determine the coefficients and eigenvalues to satisfy the BCs.
3. Solve the T, and plug the eigenvalue into the T.
4. Write the u(x,t) in series form by applying the principle of superposition.
3. Apply initial condition to u(x,t) by setting t=0 such that u(x,0)=f(x).
4. determine the Fourier coeffs.

[Benzoate]

1. After SOV, you'll get 2 ODEs. One for spatial the other for time, then you select proper eigenvalue such that the PDE has nontrivial and physically possible solution considering the given BCs. In the case of
\dfrac {T'}{T}= k \dfrac {X''}{X}=- \lambda with Dirichlet boundary conditions, the possible eigenvalue is \lambda >0 (easy to find out by plugging the BCs in).


Should I let X= c1 cos (\mu x) + c2 sin (\mux)
and k=-\mu2 ; should I write T/T' in terms of c and k ?


2. Solve the spatial ODE, determine the coefficients and eigenvalues to satisfy the BCs.

What is the spatial ODE?

3. Solve the T, and plug the eigenvalue into the T.


Not sure how to determined the eigenvalues

[klondike]

Let's say you have a 1D heat diffusion problem

\dfrac{\partial u}{\partial t} = k \dfrac{\partial ^2u}{\partial x^2}

where 0<=x<=L, and with Dirichlet boundary conditions

u(0,t)=u(L,t)=0

and initial temperature distribution

u(x,0)=f(x)


To SOV set:

u(x,t)=\phi(x)T(t)

We could set the separation constant -\lambda
you get

\dfrac{\phi^''}{\phi}=\dfrac{1}{k} \dfrac {T^'}{t}=- \lambda

you end up having 2 ODEs

\phi^{''}+\lambda \phi = 0

and

T^{'}+\lambda kt = 0


\phi^{''}+\lambda \phi = 0 is the spatial ODE. it's the simplest form of Sturm-Liouville Eigenvalue problem.
we can prove that the only way to get nontrivial solution is when \lambda>0, and the general solution to the spatial problem is

\phi(x)=c1 sin(\sqrt{\lambda}x)+c2 cos(\sqrt{\lambda}x)

To satisfy the Dirichlet BC, it requires c2=0 and

\lambda=(\dfrac{n \pi}{L})^2

now you have

\phi(x)= cn \sin(\dfrac{n\pi x}{L})

the solution to T is

T(t)=Bne^{-k \lambda t}

Plug lambda into T(t) and apply principle of superposition you have

u(x,t)=\sum _{n=1} ^{\infty} A_{n}sin \dfrac{n \pi x}{L}e^{-k(n \pi/L)^2t}


Now you need to make u(x,t) satisfy the initial condition and determine the A_{n}. Let us know if you have further problems.


Should I let X= c1 cos (\mu x) + c2 sin (\mux)
and k=-\mu2 ; should I write T/T' in terms of c and k ?



What is the spatial ODE?


Not sure how to determined the eigenvalues