Verification of Steady and Unsteady Heat Equation Solutions on a Finite Interval

In summary, the conversation discusses a tricky heat equation question that involves solving for the steady and unsteady solutions. Part (a) involves showing that a specific function, theta_0, is the steady solution of the heat equation with given boundary conditions. Part (b) requires using separation of variables to find the unsteady solution, which is given as a sum of terms involving exponential and trigonometric functions. The person asking for help has already attempted part (a) and is seeking guidance for part (b).
  • #1
mathfied
16
0
Hi. Having problems with this tricky Heat Equation Question. Managed to do part (a) and would appreciate verification that it's right.
But I can't manage to finish off the second part. I've started it off so please do advice me. Thanks a lot!

QUESTION:
-----------------------------------------


(a)

Show that the steady solution (which is independent of t) of the heat equation,

[itex]
\frac{{\partial ^2 \theta }}
{{\partial x^2 }} = \frac{1}
{{\alpha ^2 }}\frac{{\partial \theta }}
{{\partial t}}
[/itex] where [itex]\alpha[/itex] is a constant, on the interval:
[itex]
- L \leqslant x \leqslant L
[/itex] with conditions: [itex]\begin{gathered} \theta ( - L,t) = 0 \hfill \\
\theta (L,t) = T \hfill \\
\end{gathered} [/itex] is: [itex]\theta = \theta _0 (x) = \frac{{T(1 + \frac{x}
{L})}}
{2}[/itex]

(b)
Use methods of separation of variables to show that the unsteady solution for
[itex]\theta = (x,t)[/itex] with conditions: [itex]\theta ( - L) = T,{\text{ }}\theta (L) = 0,{\text{ }}\theta (x,0) = \theta _0 (x){\text{ is}}[/itex]:

[itex]
\theta (x,t) = \frac{T}
{2}\left[ {1 - \frac{x}
{L}} \right] - 2T\sum\limits_{n = 1}^\infty {\left[ {e^{ - \alpha ^2 n^2 \pi ^2 t/L^2 } \frac{{( - 1)^n \sin (n\pi x/L)}}
{{n\pi }}} \right]}[/itex]




My attempt:

Part (A)::
-----------------------
Using separation of variables:
[itex]\begin{gathered}
Let:\theta = (x,t) = X(x)T(t) \hfill \\
X(x) = Ax + B \hfill \\
T(t) = C \hfill \\
so:X(x)T(t) = (Ax + B)(C) \hfill \\
\end{gathered} [/itex]

now to use the conditions:
----------------------------------
[itex]\begin{gathered}
when:\theta ( - L,t) = 0, \hfill \\
A(x + L) + B = 0 \hfill \\
A(0) + B = 0,so:B = 0 \hfill \\
\end{gathered} [/itex]


[itex]\begin{gathered}
when:\theta (L,t) = T \hfill \\
A(x + L) + 0 = T \hfill \\
A(2L) = T \hfill \\
A = \frac{T}
{{2L}} \hfill \\
\end{gathered}[/itex]


[itex]\begin{gathered}
so: \hfill \\
\theta _0 (x) = \frac{T}
{{2L}}(x + L) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
\end{gathered}[/itex]

That's part (A) done. is my method to approach the final answer correct?




Part (B)::
-----------------------
I'm having problems here. Can't quite finish the question off. Please could you guide me out here. Thanks.

Using separation of variables:
[itex]Let:\theta = (x,t) = X(x)T(t)[/itex]

[itex]\begin{gathered}
unsteady - solution: - \rho ^2 < 0 \hfill \\
X(x) = Ae^{i\rho x} + Be^{ - ipx} = A\cos px + B\sin px \hfill \\
T(t) = Ce^{ - \alpha ^2 \rho ^2 t} \hfill \\
so:X(x)T(t) = (A\cos \rho x + B\sin \rho x)(Ce^{ - \alpha ^2 \rho ^2 t} ) \hfill \\
\end{gathered}[/itex]

now to use the conditions:
------------------------------
[itex]\begin{gathered}
\theta ( - L) = T:so: \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos \rho ( - L + L) + B\sin \rho ( - L + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos (0) + B\sin (0))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
\left[ {A(e^{ - \alpha ^2 \rho ^2 t} )} \right] = T \hfill \\
\end{gathered}[/itex]


[itex]\begin{gathered}
\theta (L) = 0:so: \hfill \\
(A\cos (2\rho L) + B\sin (2\rho L))(e^{ - \alpha ^2 \rho ^2 t} ) = 0 \hfill \\
\end{gathered}[/itex]


[itex]
\begin{gathered}
\theta (x,0) = \theta _0 (x) = \frac{{T(1 + \frac{x}
{L})}}
{2}:so: \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(e^{ - \alpha ^2 \rho ^2 (0)} ) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(1) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
\end{gathered}[/itex]

I'm stuck here. Any ideas on how to approach the final answer (as shown on the question)?

Thanks so much.
 
Physics news on Phys.org
  • #2
mathfied said:
Hi. Having problems with this tricky Heat Equation Question. Managed to do part (a) and would appreciate verification that it's right.
But I can't manage to finish off the second part. I've started it off so please do advice me. Thanks a lot!

QUESTION:
-----------------------------------------


(a)

Show that the steady solution (which is independent of t) of the heat equation,

[itex]
\frac{{\partial ^2 \theta }}
{{\partial x^2 }} = \frac{1}
{{\alpha ^2 }}\frac{{\partial \theta }}
{{\partial t}}
[/itex] where [itex]\alpha[/itex] is a constant, on the interval:
[itex]
- L \leqslant x \leqslant L
[/itex] with conditions: [itex]\begin{gathered} \theta ( - L,t) = 0 \hfill \\
\theta (L,t) = T \hfill \\
\end{gathered} [/itex] is: [itex]\theta = \theta _0 (x) = \frac{{T(1 + \frac{x}
{L})}}
{2}[/itex]
Why? You were not asked to solve the equation, only to show that that [itex]\theta_0[/itex] is the solution. Since it is linear in x, its second derivative with respect to x is 0 and since it is "steady state" (no dependence on t) its second derivative with respect to t is trivallly 0: the two partial derivatives are equal so it satisfies the differential equation.
By setting x= L and -L it is easy to see that it satisfies the boundary conditions.

(b)
Use methods of separation of variables to show that the unsteady solution for
[itex]\theta = (x,t)[/itex] with conditions: [itex]\theta ( - L) = T,{\text{ }}\theta (L) = 0,{\text{ }}\theta (x,0) = \theta _0 (x){\text{ is}}[/itex]:

[itex]
\theta (x,t) = \frac{T}
{2}\left[ {1 - \frac{x}
{L}} \right] - 2T\sum\limits_{n = 1}^\infty {\left[ {e^{ - \alpha ^2 n^2 \pi ^2 t/L^2 } \frac{{( - 1)^n \sin (n\pi x/L)}}
{{n\pi }}} \right]}[/itex]




My attempt:

Part (A)::
-----------------------
Using separation of variables:
[itex]\begin{gathered}
Let:\theta = (x,t) = X(x)T(t) \hfill \\
X(x) = Ax + B \hfill \\
T(t) = C \hfill \\
so:X(x)T(t) = (Ax + B)(C) \hfill \\
\end{gathered} [/itex]

now to use the conditions:
----------------------------------
[itex]\begin{gathered}
when:\theta ( - L,t) = 0, \hfill \\
A(x + L) + B = 0 \hfill \\
A(0) + B = 0,so:B = 0 \hfill \\
\end{gathered} [/itex]


[itex]\begin{gathered}
when:\theta (L,t) = T \hfill \\
A(x + L) + 0 = T \hfill \\
A(2L) = T \hfill \\
A = \frac{T}
{{2L}} \hfill \\
\end{gathered}[/itex]


[itex]\begin{gathered}
so: \hfill \\
\theta _0 (x) = \frac{T}
{{2L}}(x + L) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
\end{gathered}[/itex]

That's part (A) done. is my method to approach the final answer correct?




Part (B)::
-----------------------
I'm having problems here. Can't quite finish the question off. Please could you guide me out here. Thanks.

Using separation of variables:
[itex]Let:\theta = (x,t) = X(x)T(t)[/itex]

[itex]\begin{gathered}
unsteady - solution: - \rho ^2 < 0 \hfill \\
X(x) = Ae^{i\rho x} + Be^{ - ipx} = A\cos px + B\sin px \hfill \\
T(t) = Ce^{ - \alpha ^2 \rho ^2 t} \hfill \\
so:X(x)T(t) = (A\cos \rho x + B\sin \rho x)(Ce^{ - \alpha ^2 \rho ^2 t} ) \hfill \\
\end{gathered}[/itex]

now to use the conditions:
------------------------------
[itex]\begin{gathered}
\theta ( - L) = T:so: \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos \rho ( - L + L) + B\sin \rho ( - L + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos (0) + B\sin (0))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
\left[ {A(e^{ - \alpha ^2 \rho ^2 t} )} \right] = T \hfill \\
\end{gathered}[/itex]


[itex]\begin{gathered}
\theta (L) = 0:so: \hfill \\
(A\cos (2\rho L) + B\sin (2\rho L))(e^{ - \alpha ^2 \rho ^2 t} ) = 0 \hfill \\
\end{gathered}[/itex]


[itex]
\begin{gathered}
\theta (x,0) = \theta _0 (x) = \frac{{T(1 + \frac{x}
{L})}}
{2}:so: \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(e^{ - \alpha ^2 \rho ^2 (0)} ) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(1) = \frac{{T(1 + \frac{x}
{L})}}
{2} \hfill \\
\end{gathered}[/itex]

I'm stuck here. Any ideas on how to approach the final answer (as shown on the question)?
In general you can't write the solution as you have shown. What you can do is write the solutions as a sum (possibly infinite) of such things- a Fouriere series.

I would recommend a standard method: define a new function
[tex]\phi(x,t)= \theta(x,t)+ \frac{T(1+ \frac{x}{L}}{2}- T[/itex]
Since the added function, [itex]\theta_0(x,t)- T[/itex] has second derivatives with respect to both x and t 0, this function satisfies the same differential equation as [itex]\theta(x,t)[/itex] but satisifies the boundary conditions [itex]\phi(L, t)= 0[/itex], [itex]\phi(-L,t)= 0[/itex]. Then, taking [itex]\rho[/itex] so that [itex]\rho(-L)= -\pi[/itex] and [itex]\rho(L)= \pi[/itex], the solution can be written as a sum in [itex]sin(n\rho x)[/itex]

Thanks so much.
 
  • #3
Hi. Thanks.
I understand that I have to get the general solution in a Summation form so that I can apply the Fourier series.

But I have a little difficulty on applying that new function you told me to introduce.


This is how I've been taught to work out the solution: I'll try to explain in as much detail as possible. Hopefully we can take it from there..

- First I try to simplify the general form of the equation and try and turn it into a "summation" kind of form so that it is ready to have the FOURIER SERIES applied to it. So here is what i have so far:


First I have been taught to always use the general equation form for an unsteady solution which is:
-------------------------------------------------------
[itex]\theta(x,t) = X(x)T(t) = (A\cos \rho x + B\sin \rho x)(Ce^{ - \alpha ^2 \rho ^2 t} )[/itex]

We then apply the 3 conditions on the general form to try and convert the general solution to something that resembles a summation form, so that we can apply the Fourier series to it.

now to use the First condition. Please note that I have absorbed the "C " from the general equation into the other constants:
--------------------------------------------------------------

[itex]\begin{gathered}
\theta ( - L) = T:so: \hfill \\
(A\cos \rho (x + L) + B\sin \rho (x + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos \rho ( - L + L) + B\sin \rho ( - L + L))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
(A\cos (0) + B\sin (0))(e^{ - \alpha ^2 \rho ^2 t} ) = T \hfill \\
\left[ {A(e^{ - \alpha ^2 \rho ^2 t} )} \right] = T \hfill \\
\end{gathered}[/itex]
[itex]A = \frac{T}
{{e^{ - \alpha ^2 \rho ^2 t} }}[/itex]


- What usually happens with these questions is that the first condition says [itex] \theta ( - L) = 0[/itex] instead of = T.
- If it was = 0, then the whole solution would have resulted in A=0.
- So going back to the GENERAL SOLUTION, I would have been able to cancel the term with the A becauase "0" makes the whole term go to "0".
- Instead we have A as [itex]A = \frac{T}
{{e^{ - \alpha ^2 \rho ^2 t} }}[/itex] so I can't cancel the term out and instead have to work with both terms of the general solution when applying the second condition.


Using 2nd Condition:
-------------------------------------------------------------

[itex]
\begin{gathered}
\theta (L) = 0:so: \hfill \\
(\frac{T}
{{e^{ - \alpha ^2 \rho ^2 t} }}\cos (2\rho L) + B\sin (2\rho L))(e^{ - \alpha ^2 \rho ^2 t} ) = 0 \hfill \\
\end{gathered} [/itex]

sin(2pl)=0 when [itex]2\rho L = n\pi[/itex] so [itex]\rho = \frac{n\pi}{2L}[/itex]

so now subsituting "p" back into the general equation, we get:
[itex]\left( {\sum\limits_{{\text{odd n}}}^{} {\left[ {\frac{T}
{{e^{\frac{{ - \alpha ^2 n^2 \pi ^2 t}}
{{4L^2 }}} }}_n \cos \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} + \sum\limits_{{\text{even n}}}^{} {\left[ {B_n \sin \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} } \right)e^{\frac{{ - \alpha ^2 n^2 \pi ^2 t}}
{{4L^2 }}}=0[/itex]

Now using the 3rd and final condition:
--------------------------------


[itex]\theta (x,0) = \theta _0 (x) = \frac{{T(1 + \frac{x}
{L})}}
{2}:so:[/itex]

substituting t=0 into the general form:

[itex]\left( {\sum\limits_{{\text{odd n}}}^{} {\left[ {\frac{T}
{{e^{\frac{{ - \alpha ^2 n^2 \pi ^2 (0)}}
{{4L^2 }}} }}_n \cos \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} + \sum\limits_{{\text{even n}}}^{} {\left[ {B_n \sin \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} } \right)e^{\frac{{ - \alpha ^2 n^2 \pi ^2 (0)}}
{{4L^2 }}} = \frac{{T(1 + \frac{x}
{L})}}
{2}[/itex]


anything to the power of 0 becomes "1", so the e's turn into 1.
so we have:

[itex]
\left( {\sum\limits_{{\text{odd n}}}^{} {\left[ {T_n \cos \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} + \sum\limits_{{\text{even n}}}^{} {\left[ {B_n \sin \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} } \right) = \frac{{T(1 + \frac{x}
{L})}}
{2}[/itex]


[itex]
\sum\limits_{{\text{n = 1}}}^\infty {\left[ {T_n \cos \left[ {\frac{{n\pi x}}
{{2L}}} \right] + \left[ {B_n \sin \left[ {\frac{{n\pi x}}
{{2L}}} \right]} \right]} \right]} = \frac{{T(1 + \frac{x}
{L})}}
{2}[/itex]

-------------------------------------------------------------

I now seem to have the general solution in a form that can have the Fourier series applied to it. but I've tried this all day and just can't seem to get to the final answer which is:

[itex]\theta (x,t) = \frac{T}
{2}\left[ {1 - \frac{x}
{L}} \right] - 2T\sum\limits_{n = 1}^\infty {\left[ {e^{ - \alpha ^2 n^2 \pi ^2 t/L^2 } \frac{{( - 1)^n \sin (n\pi x/L)}}
{{n\pi }}} \right]}[/itex]

any idea where to take this ? thanks so much..
 

FAQ: Verification of Steady and Unsteady Heat Equation Solutions on a Finite Interval

What is the Heat Equation Problem?

The Heat Equation Problem is a mathematical model used to describe the flow of heat in a given space over time. It is a partial differential equation that takes into account factors such as temperature, thermal conductivity, and heat sources or sinks.

What is the significance of the Heat Equation Problem?

The Heat Equation Problem is important in many areas of science and engineering, as it allows us to predict and understand how heat will flow in various systems. It is used in fields such as thermodynamics, materials science, and fluid mechanics.

What are the boundary conditions for the Heat Equation Problem?

The boundary conditions for the Heat Equation Problem typically include the initial temperature distribution, the temperature at the boundaries of the system, and any external sources or sinks of heat. These conditions are necessary for solving the equation and obtaining a solution.

What are some applications of the Heat Equation Problem?

The Heat Equation Problem has many practical applications, such as predicting the temperature distribution in buildings, designing heating and cooling systems, and understanding the flow of heat in industrial processes. It is also used in fields such as geophysics to study the temperature distribution in the Earth's crust.

What are the limitations of the Heat Equation Problem?

The Heat Equation Problem is a simplified model and does not take into account factors such as convection and radiation, which can have a significant impact on the flow of heat in certain systems. It also assumes that the material properties and boundary conditions remain constant over time, which may not always be the case in real-world situations.

Back
Top