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 , t \in [0,1] [/math]

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

[math] u(x,0) = u_0(x) , x \in [0,1] [/math]

A continuación, realizaremos el mismo análisis que en el primer apartado de este documento para el nuevo sistema planteado.

4 Cambios en las condiciones de frontera y principio del máximo

Por último, repetiremos el apartado anterior pero esta vez suponiendo que la barra está aislada térmicamente en el extremo derecho. Esto quiere decir que pasamos a tener una condicion de frontera de tipo Neumann en [math] x = 1 [/math] :

[math] u_t(1,t) = 0 [/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

Ahora, vamos a dibujar la solución de la ecuación del calor en una dimensión en el semiespacio [math] x \gt 0 [/math] con condiciones:

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

Para ello, ...

A contuación veremos la relación de la ecuación del calor con la convolución.

Si tomamos la función [math]u_0(x) = 1_{[-1,1]} [/math] como condición inicial y denotamos la solución fundamental como [math] K(x,t) [/math] , entonces:

[math] u(x,t) = \int_{-\infty}^{\infty} K(x-y,t)·u_0(y) dy [/math]

es la solución de [math] u_t - u_{xx} = 0[/math] para [math] x \in \mathbb{R} [/math]

Si representamos esta solución para distintos valores de t obtenemos las siguientes gráficas:

Solución para t = 0.001

Solución para t = 0.01

Solución para t = 0.1

Podemos concluir que ...

Para terminar, tomamos el intervalo espacial [math] (x_1,x_2) = [-1,1]^2[/math] para representar la solución fundamental de la ecuación del calor en dimensión 2 para los mismos valores de t:

Solución n=2 para t = 0.001

Solución n=2 para t = 0.01

Archivo:343ot.png

Solución n=2 para t = 0.1

Observamos que cuanto más cerca está t de 0, más singular es la solución.