Ecuación del calor. Yan, Otelo, Mika

Trabajo realizado por estudiantes
Título	Series de Fourier. Grupo 6-A
Asignatura	EDP
Curso	2023-24
Autores	Miguel Cazorla Pedraza Otelo Gallego Ayala Yan Wang
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

En este articulo primero estudiaremos la ecuación del calor en dominio acotado, y más tarde trataremos su solución fundamental en todo [math] \mathbb{R} [/math]

1 Ecuación del calor I

Vamos a dibujar diferentes soluciones de la ecuación del calor en una dimensión. Para ello, partimos del siguiente enunciado:

[math] \textit{Se considera una varilla metálica que ocupa el intervalo [0,1] y que se encuentra aislada por su superficie lateral, de manera que la conducción de calor sólo se produce}[/math] [math] \textit{en la dirección longitudinal. La temperatura inicial de la varilla es } [/math] [math] 0^oC [/math]. [math] \textit{ En el extremo izquierdo se consigue mantener la temperatura a } 0^oC \textit{mientras que en el} [/math] [math] \textit{ derecho la temperatura es de } 1^oC [/math].

Suponiendo que tanto la conductividad térmica como el calor específico tienen un valor constante [math] 1 [/math], tenemos que el anterior enunciado se ve modelizado por el siguiente sistema EDP:

[math] \left\{ \frac{dT}{dt}(x,t) - \frac{d^2T}{dx^2}(x,t)=0, \hspace{5mm} x\in[0,1],\hspace{3mm} t\gt0 \atop T(x,0)=0 , \hspace{5mm} T(0,t)=0, \hspace{5mm} T(1,t)=1 \hspace{5mm} x\in[0,1],\hspace{3mm} t\gt0 \right. [/math]

Para dar solución a este sistema, resolveremos primero el problema resultante de tomar [math] t \rightarrow \infty [/math], que implica [math] \frac{dT}{dt}(x,t) = 0 [/math] por lo que se denomina caso estacionario.

[math] \left\{ - \frac{d^2V}{dx^2}(x)=0, \hspace{5mm} x\in[0,1]\atop \hspace{5mm} V(0)=0, \hspace{5mm} V(1)=1 \hspace{5mm} x\in[0,1] \right. [/math]

La solución estacionaria es [math] V(x) = x [/math]. Veamos cómo es gráficamente:

Solución estacionaria

La solución del sistema original por separación de las variables que obtenemos es:

[math] u(x,t) = x + \sum_{k=1}^{\infty}{\frac{2(-1)^k}{k}\sin(k \pi x)e^{-k^2 \pi^2 t}} [/math]

Y su representación gráfica se aproxima bastante a la estacionaria:

Solución del problema original

2 Ecuación del calor II

En esta sección vamos a dibujar diferentes soluciones de la ecuacion del calor usando la solucion fundamental.

En primer lugar, tomamos [math]x \in [-1,1] [/math] y [math]x \in [10^{-2},1] [/math] para dibujar la solución fundamental de la ecuación del calor en dimensión 1, que como vimos en clase es:

[math] u(x,t) = \frac{1}{2 \sqrt{\pi t}}e^{-\frac{x^2}{4t}} [/math]