Heat Equation With Seperable Variables

Click For Summary

Homework Help Overview

The discussion revolves around solving the heat equation with an additional term, specifically du/dt = d²u/dx² + u, subject to boundary conditions u'(0) = u'(1) = 0 and initial condition u(x,0) = 1. Participants are exploring the implications of this modified equation on the solution process.

Discussion Character

  • Exploratory, Conceptual clarification, Mathematical reasoning, Problem interpretation

Approaches and Questions Raised

  • Participants discuss the separation of variables approach and the resulting equations, questioning the implications of the term 'u' in the heat equation. There is an exploration of eigenvalues using Sturm-Liouville theory, and some participants raise concerns about satisfying boundary conditions under certain conditions of λ.

Discussion Status

The discussion is active, with participants providing insights and raising questions about the implications of the boundary conditions and initial conditions. Some guidance has been offered regarding the nature of the heat equation and the behavior of the solution, but no consensus has been reached on the specific issues raised.

Contextual Notes

Participants note that the problem is complicated by the presence of the term 'u' in the equation, which diverges from the standard heat conduction equation. There is also mention of the orthogonality of basis functions and how this affects the coefficients in the Fourier expansion.

Rweir
Messages
1
Reaction score
0

Homework Statement


du/dt = d2u/dx2 + u
Bc: u'(0) = u'(1) = 0
Ic: u(x,0) = 1


Homework Equations



Using sturm-liouville to solve for eigenvalues.

The Attempt at a Solution



After first separating variables in the equation
we get G'/G - 1 = F'' = λ
after using Sturm-Liouville we find that
F(x) = Acos(n*Pi*x)
G(t) = Ae(-n2pi2-1)t

So after multiplying them together and then taking the initial condition of u(x,0) = 1
we get A*cos(n*pi*x) = 1 and thus the problem arises after using Fourier expansion we get A = 0 which makes everything 0. Any suggestions as to why it is coming out like this?
 
Physics news on Phys.org
first the u makes it a little different, its not your standard heat conduction equation, so you have
u_t = u_xx + u

so i get to
X'' + (λ-1)X = 0
T' + λT = 0

so you won't be able to satisfy the BCs for (λ-1)<0

now how about when λ = 1?
 
have a think about the normal heat equation
u_t = u_xx
the rate of change with time of temp is proportional to the spatial curvature of the temp, ie everything gets smoothed out...

u_t = u_xx + u
now the rate of change also has a component proportional to temp as well...

and note the IC is initially uniform everywhere

now think about your BCs
u'(0) = u'(1) = 0
which in the normal heat equation are equivalent to insulated ends
 
but also stepping back a bit, as the basis functions are orthogonal, shouldn't the coefficient in the sum be proportional to
a_n = \int dx.cos(n \pi x ) .u(x,0) = \int dx.cos(n \pi x ) .u(x,0)

which is different from every basis function having to satisfy the IC, only their sum has to
 
Last edited:

Similar threads

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