How can I solve the heat equation with fixed and varying temperatures?

Click For Summary
SUMMARY

The discussion focuses on solving the heat equation, specifically the transient heat conduction equation represented by the partial differential equation (PDE) \(\frac{\partial T}{\partial t} = K\frac{\partial ^2 T}{\partial x^2}\). The user explores solutions under fixed and varying temperature conditions, noting the lack of non-decaying traveling wave solutions and the behavior of temperature equalization over time. Key references for further understanding include "Conduction of Heat in Solids" by Carslaw and Jaeger and "Transport Phenomena" by Bird, Stewart, and Lightfoot, which provide insights into boundary conditions and solution techniques.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with thermodynamics concepts
  • Knowledge of boundary conditions in heat conduction problems
  • Basic mathematical skills for manipulating exponential functions and complex numbers
NEXT STEPS
  • Study the methods for solving the heat equation in finite and seminfinite domains
  • Learn about boundary conditions and their implications on heat conduction solutions
  • Explore the application of Fourier series in solving heat conduction problems
  • Investigate numerical methods for approximating solutions to PDEs
USEFUL FOR

Students and professionals in thermodynamics, applied mathematics, and engineering disciplines, particularly those interested in heat transfer and solving partial differential equations.

MisterX
Messages
758
Reaction score
71
I've been teaching myself some thermodynamics, and I've been thinking about solving the heat equation.

\frac{\partial T}{\partial t} = K\frac{\partial ^2 T}{\partial x^2}

I haven't taken a course in PDEs.

I have noticed that if I assume an exponential solution, there are not non-decaying traveling wave solutions; the phase and group velocities are imaginary.

i\omega = k^2K
\omega = -ik^2K

non-traveling space oscillations decay in time, which isn't surprising since temperature tends to equalize:
e^{i(kx - \omega t) } = e^{ikx }e^{-k^2K t}

If we wanted to make omega real, we could have

k = \pm a(1+i)

e^{\mp ax} e^{\pm iax }e^{-i2a^2K t}

This indicates a phase speed of 2aK and a group speed of 4aK for the non-decaying factor, if I have done everything properly.

I have also seen there are solutions like
\frac{A}{\sqrt{t}}e^{-x^2/4Kt}

Also one can add a constant to any solution.

The problem I would like to solve is this: x = 0 is fixed at some temperature Ts. For x > 0 T(x, 0) is initially some other temperature Ti. If it made the solution simpler, the drop off needn't be so sharp. The point is that I would expect to see a solution so that every point with x>0 would become arbitrarily close to Ts if enough time was passed. Also, it should move along, so that it would take longer for a place with larger x to reach a given temperature than a place with smaller x.

The purpose of this thread is to solicit help in solving this problem, or any thoughts on my ideas about solutions to the heat equation.
 
Last edited:
Science news on Phys.org
Is there another boundary, say at x = L, or is the system seminfinfinite? If there is another boundary, what is the boundary condition at x = L.

In most cases, unless you've studied PDEs, you are going to have trouble solving the transient heat conduction equation. However, how to solve these specific cases can be explained.
 
I was considering a system that is seminfinite.

Actually in the system that inspired this there is a boundary but it has no fixed value (for example the end of a rod which is free to have any temperature). But I was interested in how heat would travel through the rod, so the seminfinite case would be fine I think.
 
MisterX said:
I was considering a system that is seminfinite.

Actually in the system that inspired this there is a boundary but it has no fixed value (for example the end of a rod which is free to have any temperature). But I was interested in how heat would travel through the rod, so the seminfinite case would be fine I think.

If the rod is of finite length, then that definitely matters. An appropriate boundary condition at the far end would probably be zero heat flux ( zero temperature gradient).

If you want to find out how to solve problems like these, and want to see some solutions to standard problems (such as yours), see Conduction of Heat in Solids by Carslaw and Jaeger. Another reference that shows similar results is Transport Phenomena by Bird, Stewart, and Lightfoot.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 5 ·
Replies
5
Views
16K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 27 ·
Replies
27
Views
5K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 35 ·
2
Replies
35
Views
5K
  • · Replies 4 ·
Replies
4
Views
4K
Replies
13
Views
3K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 2 ·
Replies
2
Views
1K