How Do You Solve Heat Conduction in a Sphere with Fixed Boundary Temperature?

Click For Summary
SUMMARY

The discussion focuses on solving the heat conduction equation in a solid sphere with fixed boundary temperature T_0 and initial temperature T_1. The equation is expressed as dT/dt = D(d^2T/dr^2), where the solution is approached using separation of variables, leading to the form T = X(t)R(r). The resulting ordinary differential equations yield X(t) = C*[exp(-t/(D*k^2))] and R'' + k^2*R = 0. The challenge lies in applying the boundary conditions to determine the constants A and B, as the Laplacian in spherical coordinates complicates the solution process.

PREREQUISITES
  • Understanding of partial differential equations (PDEs)
  • Familiarity with separation of variables technique
  • Knowledge of boundary value problems in heat conduction
  • Basic concepts of spherical coordinates and Laplacians
NEXT STEPS
  • Study the method of separation of variables in depth
  • Learn about boundary value problems and their solutions
  • Explore Fourier series applications in heat conduction
  • Investigate the Laplacian operator in spherical coordinates
USEFUL FOR

Mathematicians, physicists, and engineers working on heat transfer problems, particularly those dealing with spherical geometries and boundary value conditions.

genius2687
Messages
11
Reaction score
0
A solid sphere of radius a is immersed in a vat of fluid at a temperature T_0. Heat is conducted into the sphere according to

dT/dt = D(d^2T/dr^2)
(d-> partial derivative btw)

If the temperature at the boundary is fixed at T_0 and the initial temperature of the sphere is T_1, find the temperature within the sphere as a function of time.

My reasoning

Ok. Here's my reasoning. Use a solution of the form T=X(t)R(r), and plug into the above equation to get

R''/R=X/(DX')=-k^2.

I get X(t) = C*[exp(-t/(D*k^2))]. (k^2>0 for convergence)

Then I have R'' +k^2*R = 0

so R= Acos(kr) + Bsin(kr), since k^2>0.

The Problem

Assuming that the above steps are right, we could set the boundary conditions and have Acos(ka)+Bsin(ka)=T_0. That only eliminates one variable, either A or B. I don't know where to go from there however. Should this be like a Fourier series or something?
 
Physics news on Phys.org
The Laplacian in spherical coords is not d^2T/dr^2.
It is \frac{1}{r}\left[\partial_r^2(rT)\right]
 

Similar threads

Replies
13
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
6
Views
3K
  • · Replies 6 ·
Replies
6
Views
4K
Replies
5
Views
15K
  • · Replies 29 ·
Replies
29
Views
4K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K