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

[Benzoate]

Homework Statement

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

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

where n extends from 0 to \infty

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

where the limits extend from 0 to L

Homework Equations

separation of variables

The Attempt at a Solution

u(x,t)=X(x)T(t)
u(x,0)=f(x)
u₂(x,t)=0
T'-kc²T=0
x''-kx=0
x'(0)=0
x'(L)=0
u₂(0,t)=0

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

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

not sure how to finish this derivation.

 

[gabbagabbahey]

u(x,t)=X(x)T(t)
u(x,0)=f(x)
T'-kc²T=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.

u₂(x,t)=0

Huh? What is this supposed to mean?

x'(0)=0
x'(L)=0
u₂(0,t)=0

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

u₂(0,t)=0= u₂(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)=c2case 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 u₂. Those are all mistakes. ; I was trying to say that if u(0,t)=0 then u(0,t)=0==> u_x(0,t)=0

Therefore, u_x(0,t)=0= u_x(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&#039;}{T}= k \dfrac {X&#039;&#039;}{X}=- \lambda with Dirichlet boundary conditions, the possible eigenvalue is \lambda &gt;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&#039;}{T}= k \dfrac {X&#039;&#039;}{X}=- \lambda with Dirichlet boundary conditions, the possible eigenvalue is \lambda &gt;0 (easy to find out by plugging the BCs in).

Should I let X= c1 cos (\mu x) + c2 sin (\mux)
and k=-\mu² ; 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
<br /> \dfrac{\partial u}{\partial t} = k \dfrac{\partial ^2u}{\partial x^2}<br />
where 0<=x<=L, and with Dirichlet boundary conditions
<br /> u(0,t)=u(L,t)=0<br />
and initial temperature distribution
<br /> u(x,0)=f(x)<br />

To SOV set:
<br /> u(x,t)=\phi(x)T(t)<br />
We could set the separation constant -\lambda
you get
<br /> \dfrac{\phi^&#039;&#039;}{\phi}=\dfrac{1}{k} \dfrac {T^&#039;}{t}=- \lambda<br />
you end up having 2 ODEs
<br /> \phi^{&#039;&#039;}+\lambda \phi = 0<br />
and
<br /> T^{&#039;}+\lambda kt = 0<br />

\phi^{&#039;&#039;}+\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&gt;0, and the general solution to the spatial problem is
<br /> \phi(x)=c1 sin(\sqrt{\lambda}x)+c2 cos(\sqrt{\lambda}x)<br />
To satisfy the Dirichlet BC, it requires c2=0 and
<br /> \lambda=(\dfrac{n \pi}{L})^2<br />
now you have
<br /> \phi(x)= cn \sin(\dfrac{n\pi x}{L})<br />
the solution to T is
<br /> T(t)=Bne^{-k \lambda t}<br />
Plug lambda into T(t) and apply principle of superposition you have
<br /> u(x,t)=\sum _{n=1} ^{\infty} A_{n}sin \dfrac{n \pi x}{L}e^{-k(n \pi/L)^2t}<br />

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=-\mu² ; should I write T/T' in terms of c and k ?



What is the spatial ODE?


Not sure how to determined the eigenvalues