Diferencia entre revisiones de «ECUACIÓN DEL CALOR»

De MateWiki
Saltar a: navegación, buscar
(→‎. Deducción de la Solución estacionaria)
(→‎. Deducción de la Solución estacionaria)
Línea 54: Línea 54:
 
v(0)=0, & \quad 0 < t < T, \\
 
v(0)=0, & \quad 0 < t < T, \\
  
−
v(1)=1 & \quad 0 < t < T, \\  
+
v(1)=1 & \quad 0 < t < T. \\  
  
 
\end{array} \right.  
 
\end{array} \right.  

Revisión del 21:35 4 mar 2024

Trabajo realizado por estudiantes
Título Ecuación del calor. Grupo GRwM
Asignatura EDP
Curso 2023-24
Autores Arturo Barrena y Mario Ríos
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


1 . Introducción

La visualización de soluciones de la ecuación del calor en una dimensión constituye una herramienta poderosa para desentrañar los intricados patrones de propagación térmica en sistemas unidimensionales. En este estudio, nos sumergimos en la riqueza matemática que rodea la representación gráfica de soluciones, explorando la diversidad de comportamientos térmicos que emanan de este modelo fundamental. A través de la creación visual, buscamos capturar la esencia dinámica de la ecuación del calor y ofrecer una perspectiva única sobre cómo la temperatura evoluciona en función del tiempo y la posición.

La ecuación del calor en una dimensión sirve como un punto de partida teórico sólido para comprender la difusión térmica en estructuras lineales. Al trazar diferentes soluciones, nos sumergimos en el tejido mismo de la propagación térmica, revelando detalles intrincados que, de otra manera, podrían escapar a una descripción puramente analítica. Este enfoque gráfico no solo ilustra la dinámica temporal de la temperatura, sino que también destaca la influencia de condiciones iniciales y parámetros en la evolución térmica a lo largo del espacio unidimensional.

2 . Planteamiento del sistema

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 en la dirección longitudinal. La temperatura inicial de la varilla es 0 grados. En el extremo izquierdo se consigue mantener la temperatura a 0 grados mientras que en el derecho la temperatura es siempre de 1 grado. De acuerdo con todas estas indicaciones, el sistema de EDP que modeliza el comportamiento de la temperatura es el siguiente:

[math]\left \{ \begin{array}{ll} \frac{\partial u}{\partial t}- \frac{\kappa}{c} \frac{\partial ^2 u}{\partial x^2}=0 & \quad 0 \lt x \lt 1, 0 \lt t \lt T, \\ u(x,0)=0, & \quad 0 \lt x \lt 1, \\ u(0,t)=0, & \quad 0 \lt t \lt T, \\ u(1,t)=1 & \quad 0 \lt t \lt T, \end{array} \right. [/math]

donde [math]u(x,t)[/math] es la temperatura en la posición [math]x[/math] en el tiempo [math]t[/math], [math]\kappa\gt0[/math] es la constante de conductividad térmica y [math]c[/math] es el calor específico del material, el cual también se considera constante. Cabe destacar que por simplificar el problema se suelen tomar los valores [math]\kappa[/math] y [math]c[/math] como constantes tal y como se ha hecho en este caso, ya que en bastantes casos de interés es razonable ignorar su variación. Si simplificamos aún más, suponiendo que ambas constantes son igual a 1, se tiene el siguiente sistema:

[math]\left \{ \begin{array}{ll} \frac{\partial u}{\partial t}-\frac{\partial ^2 u}{\partial x^2}=0 & \quad 0 \lt x \lt 1, 0 \lt t \lt T, \\ u(x,0)=0, & \quad 0 \lt x \lt 1, \\ u(0,t)=0, & \quad 0 \lt t \lt T, \\ u(1,t)=1 & \quad 0 \lt t \lt T. \end{array} \right. [/math]

3 . Deducción de la Solución estacionaria

Antes de tratar de resolver el sistema, intentemos conjeturar qué podría suceder. Dado que la temperatura en el lado derecho es mayor que la del lado izquierdo según las condiciones frontera, el calor comenzará a fluir desde el extremo más caliente, aumentando la temperatura dentro de la barra y causando una salida de calor hacia el lado frío. Por otro lado, el aumento interior de la temperatura hace que el flujo de entrada de calor disminuya con el tiempo (ya que el flujo de entrada es exactamente la derivada de la temperatura con respecto a x cambiada de signo), mientras que el flujo de salida aumenta. Esperamos que tarde o temprano los dos flujos se equilibren entre sí y que la temperatura eventualmente alcance una distribución de estado estacionario. Por tanto, es razonable suponer que si se hace el límite cuando el tiempo tiende a infinito de la función temperatura [math]u[/math], se obtendrá una función que tan sólo depende de [math]x[/math], ya que la temperatura debe haber llegado al estado estacionario mencionado en algún momento. Por tanto, bajo esta suposición, la siguiente función [math]v[/math] está bien definida (depende sólo de x):

[math] v(x)=\lim_{{t \to \infty}} u(x,t). [/math]

Además, cumple el sistema de EDO:

[math]\left \{ \begin{array}{ll} \frac{\partial ^2 v}{\partial x^2}=0 & \quad 0 \lt x \lt 1, \\ v(0)=0, & \quad 0 \lt t \lt T, \\ v(1)=1 & \quad 0 \lt t \lt T. \\ \end{array} \right. [/math]

Si resolvemos este problema de valor inicial obtendremos:

[math]v(x)=x. [/math]

Es decir, que una vez que la temperatura se estabilice y deje de variar con el tiempo, seguirá el comportamiento de esta función lineal [math]v[/math].

[math][/math]

4 Resolución de la parte no estacionaria

Homogeneizamos las condiciones frontera utilizando la solución estacionaria que hemos obtenido en el anterior apartado, es decir plantear un problema equivalente con condiciones de frontera homogéneas. Por aclarar la expresion final de la solución se compondrá de dos partes la parte estacionaria y la parte no estacionaria, es decir que lo podemos escribir de la siguiente forma:

[math]u(x,t)= u_{s}(x) + v(x,t) [/math]

donde v(x,t) corresponde a la parte no estacionaria de la solución que es la abordaremos en esta sección que despejando de la anterior expresión será:

[math] u(x,t)-u_{s}(x)= v(x,t) [/math]

De esta forma podemos llegar al sistema homogéneo que será de la forma:

[math]\left \{ \begin{array}{ll} \frac{\partial c}{\partial t}-\frac{\partial ^2 v}{\partial x^2}=0 & \quad 0 \lt x \lt 1, 0 \lt t \lt T, \\ v(x,0)=-x, & \quad 0 \lt x \lt 1, \\ v(0,t)=0, & \quad 0 \lt t \lt T, \\ v(1,t)=1 & \quad 0 \lt t \lt T. \end{array} \right. [/math]

Ahora procedemmos a la resolución del sistema por medio del método de separación de variables, es decir consideramos que la solución es de la forma v(x,t)=X(x)T(t). Resolviendo obtendremos finalmente:

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

De modo que la expresión para la solución final será:

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

5 Cambio de condiciones de frontera

Ahora cambiamos las condiciones en la frontera, en lugar de suponer que la temperatura es 0 grados en el extremo derecho, suponemos que la barra está aislada térmicamente en dicho extremo. De este modo obtendremos el siguiente sistema

[math]\left \{ \begin{array}{ll} \frac{\partial u}{\partial t}- \frac{\partial ^2 u}{\partial x^2}=0 & \quad 0 \lt x \lt 1, 0 \lt t \lt T, \\ u(0,t)=0, & \quad 0 \lt t \lt T, \\ u_{x}(1,t)=0 & \quad 0 \lt t \lt T, \\ u(x,0)=max \{0,1-4|x-1/2|\}, & \quad 0 \lt x \lt 1. \end{array} \right. [/math]

Se puede observar que el sistema ya está homogeneizado ya que las condiciones de frontera son nulas, por tanto solo haría falta resolver el problema por el método de separación de variables:

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

Véase que no podemos trabajar con la base trigonométrica usual [math]\{ 1, sen(\pi kx), cos((\pi kx)\}_{k=1}^{\infty} [/math] ya que no nos permite describir adecuadamente la solución. En su lugar, vamos a trabajar con la base [math]\{ \frac{1}{2}, sen((\frac{\pi}{2}+k\pi)x), cos((\frac{\pi}{2}+k\pi)x) \}_{k=1}^{\infty} [/math] que ha sido modificada de tal forma que podamos representar adecuadamente la solución .Extendiendo de forma impar la función [math]max \{0,1-4|x-1/2|\}[/math], vamos a tener la base [math]\{sen((\frac{\pi}{2}+k\pi)x)\}_{k=1}^{\infty} [/math] pues forma un conjunto completo en [math][0,1][/math] y es ortogonal. Verificamos que la base es ortonormal:

[math] \lt sen((\frac{\pi}{2}+n\pi)x), sen((\frac{\pi}{2}+m\pi)x)\gt_{L^2(0,1)} = \int_{0}^{1}sen((\frac{\pi}{2}+n\pi)x) \cdot sen((\frac{\pi}{2}+m\pi)x) dx = 0[/math]

Gracias a esta comprobación podemos ver que la base es ortogonal y completa por tanto podemos calcular los [math] c_{k} [/math] correspondientes a la solución que son los coeficientes de Fourier en esta nueva base. Planteamos la expresión de los [math] c_{k} [/math]:

[math] c_{k} = \int_{0}^{1}sen((\frac{\pi}{2}+n\pi)x) \cdot f(x) dx [/math]

Si realizamos la integral obtenemos:

[math] c_{k} = (8 \pi n cos(\pi n) - 8 sen(\pi n))/(\pi^2 n^2) + 2[/math]

También aportaremos una aproximación de estos coeficientes usando el método numérico del trapecio, para poder representar adecuadamente la solución.

6 Principio del Máximo

El teorema del principio del máximo nos permite asegurar la unicidad del sistema. Dicho teorema nos dice lo siguiente:

Sea [math] u \in C^{2,1} (Q_T) \cap C(\overline{Q_T})[/math] que verifica [math]u_t - \Delta u\leq 0 [/math] en [math]Q_T[/math].Entonces u alcanza su máximo en la frontera parabólica:

[math]\max \limits_{(x,t) \in \overline{Q_T}} u = \max \limits_{(x,t) \in \partial _P Q_T} u [/math]

Nuestra solución correspondiente al apartado anterior cumplen todas las hipótesis del teorema del máximo, en primer lugar [math] u \in C^{2,1} (Q_T) \cap C(\overline{Q_T})[/math] además por el apartado anterior sabemos que la u que es la solución correspondiente al sistema cumple la expresión [math]u_t - \Delta u\leq 0 [/math]. Además la condición inicial converge uniformemente de manera que la condición que necesitamos para la continuidad se cumple a la perfección. (HAY QUE JUSTIFICAR POR QUE LA CONDICION INICIAL CONVERGE UNIFORMEMENTE).

De esta forma si consideramos dos soluciones u y v de tal forma que definimos w=u-v que satisface:

[math]\left \{ \begin{array}{ll} w_t - \Delta w =0 & \quad x \in (0,1), T \gt 0 \\ w(0,t)=0, & \quad x \in (0,1) , \\ w(x,T)=0, & \quad x \in \{0,1\}, T \gt 0. \end{array} \right. [/math]

Aplicando el principio del máximo obtenemos: [math] \min \limits_{(x,t) \in \partial _P Q_T}w \leq w \leq \max \limits_{(x,t) \in \partial _P Q_T}w [/math] de aqui deducimos que [math]w \equiv 0[/math] por tanto podemos concluir que u=v, quedando mostrado que la solución del sistema es única.

Obtenido de «https://mat.caminos.upm.es/w/index.php?title=ECUACIÓN_DEL_CALOR&oldid=68692»