Diferencia entre revisiones de «Approximation of the heat equation: Fourier method»

De MateWiki
Saltar a: navegación, buscar
Línea 24: Línea 24:
  
 
The main point is to observe that, if <math>\varphi(x)</math> is solution of the eigenvalue problem
 
The main point is to observe that, if <math>\varphi(x)</math> is solution of the eigenvalue problem
 +
 
<math>
 
<math>
 
\left\{ \begin{array}{l}
 
\left\{ \begin{array}{l}
Línea 29: Línea 30:
 
\varphi(0)=0, \quad \varphi(L)=0, \end{array} \right.
 
\varphi(0)=0, \quad \varphi(L)=0, \end{array} \right.
 
</math>
 
</math>
 +
 
for some <math>\lambda</math>, then
 
for some <math>\lambda</math>, then
 +
 
<math>
 
<math>
 
u(x,t)=\varphi(x) e^{-\lambda t}
 
u(x,t)=\varphi(x) e^{-\lambda t}
 
</math>
 
</math>
 +
 
is a solution of the heat equation <math>u_t-u_{xx}=0</math> and the boundary conditions <math>u(0,t)=u(L,t)=0.</math>
 
is a solution of the heat equation <math>u_t-u_{xx}=0</math> and the boundary conditions <math>u(0,t)=u(L,t)=0.</math>
 +
 +
The eigenvalues of the above eigenvalue problem are <math>\lambda_k=\frac{k^2 \pi^2}{L^2}</math> with <math>k=1,2,3,...</math> and the associated eigenfunctions are
 +
 +
<math>
 +
\varphi_k=\sin(\frac{k\pi x}{L}).
 +
</math>
 +
 +
Therefore, all these functions are solutions of the heat equation and boundary conditions
 +
 +
<math>
 +
u_k (x,t)= e^{-k^2\pi^2 /L^2 t}\sin(\frac{k\pi x}{L})
 +
</math>
  
  

Revisión del 10:17 26 abr 2016

Let [math]{\bf u(x_1,x_2,x_3,t)}[/math] be the temperature at the point [math](x_1,x_2,x_3)\in [/math] bar, and time [math]t\gt0[/math].

Assume that the cross section is so small that we can consider the bar as an unidimensional object in the interval [math]x\in [0,L][/math], [math] u(x_1,x_2,x_3,t)=u(x_1,t)=u(x,t). [/math]

If we assume that the extreme are at zero temperature, the system of equations for u is given by [math] \left\{ \begin{array}{l} u_t-u_{xx}=0, \qquad x\in(0,L), \qquad t\gt0, \\ u(0,t)=0, \qquad t\gt0, \\ u(L,t)=0, \qquad t\gt0, \\ u(x,0)=u^0(x), \qquad x\in(0,L). \end{array} \right. [/math] where [math]u_0(x)[/math] is a function that describes the initial temperature of the bar.

We describe below the Fourier method to approximate the solutions of this system.

The main point is to observe that, if [math]\varphi(x)[/math] is solution of the eigenvalue problem

[math] \left\{ \begin{array}{l} \varphi''(x)+\lambda \varphi(x)=0, \\ \varphi(0)=0, \quad \varphi(L)=0, \end{array} \right. [/math]

for some [math]\lambda[/math], then

[math] u(x,t)=\varphi(x) e^{-\lambda t} [/math]

is a solution of the heat equation [math]u_t-u_{xx}=0[/math] and the boundary conditions [math]u(0,t)=u(L,t)=0.[/math]

The eigenvalues of the above eigenvalue problem are [math]\lambda_k=\frac{k^2 \pi^2}{L^2}[/math] with [math]k=1,2,3,...[/math] and the associated eigenfunctions are

[math] \varphi_k=\sin(\frac{k\pi x}{L}). [/math]

Therefore, all these functions are solutions of the heat equation and boundary conditions

[math] u_k (x,t)= e^{-k^2\pi^2 /L^2 t}\sin(\frac{k\pi x}{L}) [/math]