Heat Equation with Periodic Boundary Conditions

Click For Summary
SUMMARY

The discussion focuses on solving the heat equation on a ring of radius R using separation of variables, yielding the general solution U(s,t) = (Acos(λs) + Bsin(λs))exp(-λ²kt). The periodicity condition U(s',t) = U(s' + 2πR, t) leads to the conclusion that λ must equal k/R, where k is an integer. The confusion arises when considering the periodicity for n complete revolutions, which introduces the condition λ = m/(nR), leading to non-eigenvalue solutions. The key takeaway is that while 2π periodic functions are inherently 4π periodic, the reverse is not true, clarifying the constraints on λ.

PREREQUISITES
  • Understanding of the heat equation and its solutions
  • Familiarity with separation of variables technique
  • Knowledge of periodic boundary conditions in differential equations
  • Basic concepts of eigenvalues and eigenfunctions
NEXT STEPS
  • Study the implications of periodic boundary conditions in partial differential equations
  • Explore the concept of eigenvalues in the context of Sturm-Liouville problems
  • Learn about the diffusion equation and its applications in physics
  • Investigate the mathematical properties of Fourier series and their relation to periodic functions
USEFUL FOR

Mathematicians, physicists, and engineering students interested in solving differential equations, particularly those involving periodic boundary conditions and eigenvalue problems.

StretchySurface
Messages
5
Reaction score
0
Homework Statement
Find General Solution of Heat Equation on a Ring.
Relevant Equations
See below
I'm solving the heat equation on a ring of radius ##R##. The ring is parameterised by ##s##, the arc-length from the 3 o'clock position. Using separation of variables I have found the general solution to be:

$$U(s,t) = S(s)T(t) = (A\cos(\lambda s)+B\sin(\lambda s))*\exp(-\lambda^2 kt)$$.

Here is my confusion:
Periodicity demands that ##U(s',t) = U(s' + 2\pi R, t).## This means that
$$\cos(\lambda s') = cos(\lambda s' + \lambda*2\pi R)$.$

This finally gives condition that ##\lambda = k/R%## where ##k## is an integer. This is the correct answer as told by my professor/textbook.

However, periodicity can also demands ##U(s',t) = U(s' + 2\pi nR, t)##, i.e revolving around the ring ##n## complete times. But this leads to the condition ##2\pi \lambda nR = 2 \pi m## where ##m## and ##n## are arbitrary integers. This gives us ##\lambda = m/(nR)##. This doesn't seem to be the correct answer but I don't understand why and rational number divided by ##R## is not an eigenvalue.
 
Physics news on Phys.org
Something with a period of ##4\pi## does not necessarily satisfy the condition that ##f(\theta)=f(\theta+2\pi)## !
 
Does that mean the following is untrue: We require ##U(s′,t)=U(s′+4πR,t)##?
 
It is a true statement, but:
It's a necessary but insufficient condition.
 
Okay, so I don't really see the point you're trying to convey and how it relates to my question. Below, I'll write out a set steps in reasoning to see if you can pinpoint the error I'm making.

Let's only consider the case ##n = 2## i.e traveling around the ring twice.
Periodicity demands that ##U(s',t)=U(s'+4πR,t)##.
This means that ##2\lambda R = k##, for any integer ##k##.
Thus ##\lambda_{k} = k/(2R)## which includes some values not mentioned in the notes I'm working from.
 
StretchySurface said:
some values not mentioned
Yes. And those do NOT satisfy the condition ##U(s′,t)=U(s′+2\pi R,t)##. Do you propose to drop that condition ?
 
Put slightly differently, 2pi periodic functions are automatically 4pi periodic but 4pi periodic functions are not necessarily 2pi periodic.
 
Ah okay, now I understand! Thank you very much!
 
Btw if ##\lambda## is the inverse of a length, and ##k## is an integer (I suppose dimensionless), that exponential doesn't look good to me.

Edit: ok the ##k## in the exponential is the diffusion coefficient (I suppose)... my bad
 
  • #10
dRic2 said:
Btw if ##\lambda## is the inverse of a length, and ##k## is an integer (I suppose dimensionless), that exponential doesn't look good to me.

Edit: ok the ##k## in the exponential is the diffusion coefficient (I suppose)... my bad
I think you can be excused as the OP used ##k## both for the diffusion coefficient as well as an integer ... :wink:
 

Similar threads

Replies
6
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
9
Views
2K
Replies
0
Views
1K
  • · Replies 1 ·
Replies
1
Views
3K
Replies
6
Views
2K