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 su representación gráfica en 3D:

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

A continuación observamos en una gráfica el flujo de calor saliente y entrnte en ambos extremos de la varilla a lo largo del tiempo:

YAN

2 Cambio en el coeficiente de conductividad térmica

En este apartado vamos a variar el coeficiente de conductividad térmica de 1 a 1/2. Para empezar, vamos qué implica este cambio en el sistema de ecuaciones:

MIKA

Por tanto, la solución para este caso es:

[math] u(x,t) = [/math]

Y su representación gráfica es la siguiente:

Analicemos también la diferencia entre la solución y el estado estacionario en [math] x = 1/2 [/math] para [math] t \in [0,1] [/math] considerando ambos coeficientes de conductividad:

Para [math] K = 1 [/math]

Y para [math] K = 1/2 [/math], se obtiene:

Por tanto, la conclusión obtenida es

3 Cambio en la condición inicial

Pasamos ahora a variar la condición inicial. Para ello, consideremos que la temperatura en el extremo derecho es también 0 °C:

[math] u(1,t) = 0 [/math]

y consideramos la función [math] u_0(x) = \max{0, 1-4·|(x-1/2|} [/math]

5 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{1}{4t}x^2} [/math]

Solución fundamental