Solving the Heat Equation: Investigating an Error

Click For Summary
SUMMARY

The forum discussion centers on solving the heat equation given by the formula \( \frac{1}{a^2} \theta_t = \theta_x \) with specific boundary conditions \( \theta(0,t) = \theta(L,t) = 0 \) and initial conditions \( \theta(x,0) = T_0 \sin\left(\frac{x\pi}{L}\right) \). The user derived the general solution as \( \theta(x,t) = \sum_{n=1}^\infty B_n \sin\left(\frac{n\pi x}{L}\right) e^{-\frac{a^2 n \pi^2}{L} t} \) but encountered an issue determining the coefficient \( B_1 \) after separating it from the summation. The solution suggests using Fourier sine series to expand \( T_0 \sin\left(\frac{\pi x}{L}\right) \) and equate coefficients for a valid approach.

PREREQUISITES
  • Understanding of partial differential equations, specifically the heat equation.
  • Familiarity with Fourier series, particularly Fourier sine series.
  • Knowledge of orthogonality conditions for sine functions.
  • Basic integration techniques involving trigonometric identities.
NEXT STEPS
  • Learn how to derive Fourier sine series expansions for functions.
  • Study the orthogonality properties of sine functions in the context of Fourier series.
  • Explore integration techniques for products of sine functions, including relevant trigonometric identities.
  • Investigate methods for determining coefficients in series solutions of partial differential equations.
USEFUL FOR

Students studying applied mathematics, particularly those focusing on partial differential equations and Fourier analysis, as well as educators looking for examples of solving the heat equation.

catcherintherye
Messages
47
Reaction score
0

Homework Statement



I am solving the heat equation [tex]1/a^2\theta_t = \theta_x[/tex]

boundary conditions are [tex]\theta(0,t) = \theta(L,t) = 0 t > 0[/tex]

initial conditions are [tex]\theta(x,0) = T_0sin(x\pi/L)[/tex]


now I have derived the steady solution to be 0 and I have derived that the general solution will be of the form

[tex]\theta(x,t) = \sum_{n=1}^\infty\B_n sin\frac{n\pix}{L}exp\frac{-a^2n^\pi^2}{L}t[/tex]

I am next required to determine the sequence [tex]{B_n} \\<br /> <br /> \mbox{now usually I would proceed in the following manner, using the initial conditions}\\<br /> \theta(x,0) = \sum_{n=1}^\infty B_n sinn\pix/L = T_0sin(\frac{\pix}{L})[/tex]

I would then multiply both sides by [tex]\int_{0}^{L}sin\frac{mx\pi}{L}[/tex]

but instead i decided first to separted term [tex]B_1[/tex] and cancel term sin(mx\pi/L)

so I proceeded using the cancellation by the orthogonality conditions and derived the following

[tex]B_mL/2 = (T_0 -B_1)\int_{0}^{L}sin\frac{mx\pi}{L}\\<br /> <br /> = \frac{(B_1 -T_0)}{m\pix}\left[cosmx\pi/L\right]_{0}^{L}\\<br /> = \frac{B_1-T_0}{m\pi}\left[(-1)^m -1][/tex]

you see the problem is I end up with [tex]B_1[/tex] undetermined
[tex]B_1 = \frac{-4T_0}{m\pi -4}[/tex]

where did I go wrong/ what should I have done differently??:frown:





Homework Equations





The Attempt at a Solution

 
Last edited:
Physics news on Phys.org
You didn't do anything wrong except that you guarenteed you would not be able to determine B1 when you "decided first to separted term [itex]B_1[/itex] and cancel term sin(mx\pi/L)"!

There are many different ways to do this but I would probably expand [itex]T_0 sin(\pi x/L)[/itex] in a Fourier sine series in [itex]\pi x/L[/itex] and equate coefficients.
 
mmm...I am not familiar with this method, you see the reason I took out term B1 was because I knew it would allow me to proceed by multplying by

[tex]sin\frac{m\pix}{L}[/tex]and then integrating from 0 to L. Without taking out term B1 this would eventually lead me to do the following integral
[tex]L/2B_m = T_0\int_{0}^{L} sin(\frac{m\pix}{L})sin(\frac{pix}{L})dx[/tex]

is this also a valid method of solution and if so How do I do the integral? Do I have to use some sort of trig identity
 
Last edited:

Similar threads

  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 16 ·
Replies
16
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K