Quantcast Heat equation...Fourier Series Text - Physics Forums Library

PDA

View Full Version : Heat equation...Fourier Series


Link-
Jul5-08, 01:29 AM
I'm trying to solve a basic heat equation \frac{\partial T}{\partial t}=\alpha \frac{\partial^2 T}{\partial x^2} I manage to get T=X(x)\tau(t)
Then \tau(t)=A*e^{-\alpha*\lambda^2*t} and X(x)=C*sin(\lambda x) where \lambda=\pi/Ln n=1,2,3,...

From here I don't know how or why I get to a Fourier Serie. Like this u(t,x) = \sum_{n = 1}^{+\infty} D_n \left(\sin \frac{n\pi x}{L}\right) e^{-\frac{n^2 \pi^2 kt}{L^2}} \\ and \\ D_n = \frac{2}{L} \int_0^L f(x) \sin \frac{n\pi x}{L} \, dx.

Can someone explain why I have this Fourier Serie?

Thanks Link

Hootenanny
Jul5-08, 06:25 AM
\tau(t)=A*e^{-\alpha*\lambda^2*t}

X(x)=C*sin(\lambda x) where \lambda=\pi/Ln n=1,2,3,...

From here, recall that when you seperated the variables you stated that:

T=X(x)\tau(t)

Which leads to:

T = u_n(x,t) = A_n*e^{-\alpha\left(\frac{n\pi}{L}\right)^2 t}\left(C_n*sin(\frac{n\pi}{L} x)\right)

This is the general solution for each value of n\in\mathbb{Z}^+. However, the general solution to the heat equation is a superposition of all possible modes (or values of n). In otherwords, to find the general solution to the original PDE you need to sum over all possible n:

u(x,t) = \sum^{\infty}_{n=1}\left\{ A_n*e^{-\alpha\left(\frac{n\pi}{L}\right)^2 t}\left[C_n*sin\left(\frac{n\pi}{L} x\right)\right]\right\}


Combining the coefficients such that A_n\cdot C_n = D_n:

u(x,t) = \sum^{\infty}_{n=1}D_ne^{-\alpha\left(\frac{n\pi}{L}\right)^2 t}\sin\left(\frac{n\pi}{L} x\right)


In general the Fourier series of a function my be written as the sum of sines and cosines, but since your function is odd (i.e. only has sines and exponents) the coefficients of the cosine component are identically equal to zero. Therefore, you are simply left with the coefficients of the sine component which are given by Euler's formula:

b_n = \frac{2}{L}\int_0^Lf(x)\sin\left(\frac{2n\pi}{T}x\right)

Where T is the period, which in this case is 2L. Hence:

D_n = \frac{2}{L}\int_0^Lf(x)\sin\left(\frac{n\pi}{L}x\right)

I hope this helps.

Link-
Jul6-08, 02:07 PM
Thanks a lot Hootenanny.

klondike
Jul7-08, 10:53 PM
I'm trying to solve a basic heat equation \frac{\partial T}{\partial t}=\alpha \frac{\partial^2 T}{\partial x^2} I manage to get T=X(x)\tau(t)
Then \tau(t)=A*e^{-\alpha*\lambda^2*t} and X(x)=C*sin(\lambda x) where \lambda=\pi/Ln n=1,2,3,...

From here I don't know how or why I get to a Fourier Serie. Like this u(t,x) = \sum_{n = 1}^{+\infty} D_n \left(\sin \frac{n\pi x}{L}\right) e^{-\frac{n^2 \pi^2 kt}{L^2}} \\ and \\ D_n = \frac{2}{L} \int_0^L f(x) \sin \frac{n\pi x}{L} \, dx.

Can someone explain why I have this Fourier Serie?

Thanks Link
I'm assuming f(x)=u(0,x) which is the initial condition of the PDE. So the solution u(t,x) must satsify IC (well I'm more used to u(x,t) but what the heck...)


u(t,x) = \sum_{n = 1}^{+\infty} D_n \left(\sin \frac{n\pi x}{L}\right) e^{-\frac{n^2 \pi^2 kt}{L^2}}

so by setting t=0, you got

u(0,x) = f(x)=\sum_{n = 1}^{+\infty} D_n \left(\sin \frac{n\pi x}{L}\right) }

Multiply both sides by sin(mPI/L)

f(x) \sin \frac{m\pi x}{L}=\sum_{n = 1}^{+\infty} D_n \left(\sin \frac{m\pi x}{L} \sin \frac{n\pi x}{L}\right) }

perform a term by term integration on RHS over 0=<x=<L

\int _{0}^{L}f(x) \sin \frac{m\pi x}{L}dx= \sum_{n = 1}^{+\infty} \int _{0}^{L}D_n \left(\sin \frac{m\pi x}{L} \sin \frac{n\pi x}{L}\right) }dx

because

\int _{0}^{L}\left(\sin \frac{m\pi x}{L} \sin \frac{n\pi x}{L}\right) }dx = 0 \quad when \quad m \neq n

and

\int _{0}^{L}\sin ^{2} \left(\frac{m\pi x}{L}\right) }dx = \frac{L}{2} \quad when \quad m=n \neq0

It reduced to

\int _{0}^{L}f(x) \sin \frac{m\pi x}{L}dx=D_{n}\frac{L}{2}

now replace m with n (because m=n)

D_{n}=\frac{2}{L}\int _{0}^{L}f(x) \sin \frac{n\pi x}{L}dx


Hope that helps, at least partially.