Diferencia entre revisiones de «Ecuación del calor (PPAD)»

De MateWiki
Saltar a: navegación, buscar
Línea 50: Línea 50:
  
 
Esta nueva ecuación garantiza que la velocidad de propagación del calor se produzca con velocidad finita, por lo que se utiliza en modelos donde el efecto de propagación instantánea tiene importancia.  
 
Esta nueva ecuación garantiza que la velocidad de propagación del calor se produzca con velocidad finita, por lo que se utiliza en modelos donde el efecto de propagación instantánea tiene importancia.  
 +
 +
Si resolvemos esta ecuación por separación de variables, tenemos las siguientes ecuaciones para las función espacial <math>X(x)</math> y temporal <math> T(t)</math>:
 +
 +
<center><math> X'' + \frac{\lambda}{\alpha}X=0, \hspace{5mm} X(0)=X(L)=0 </math></center> 
 +
 +
<center><math> \tau T'' + T' +\lambda T =0, \hspace{5mm}  </math></center>
 +
 +
La ecuación espacial tiene como solución general <math> X(x) = Asen(kx) + Bcos(kx) </math>, con <math> k^{2}=\frac{\lambda}{\alpha} <math>. Usando las condiciones de Dirichlet, tenemos que <math>B=0</math> y nos queda la colección de soluciones
 +
 +
<center><math> X_{n}(x)=A_{n}sen(\frac{n\pi}{L}x) </math></center> 
  
 
=. Referencias =
 
=. Referencias =

Revisión del 14:33 16 mar 2025

Trabajo realizado por estudiantes
Título Ecuación del calor (Grupo PPAD).
Asignatura EDP
Curso 2024-25
Autores Pablo Vidal Nacle, Pablo Maestro Fernández, Alex Heredero Santamaría, Diego Moñino Vizmanos
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura


1 . Introducción

La ecuación del calor clásica es una herramienta crucial para los ingenieros y físicos. Ahora bien, esta ecuación plantea un serio problema conocido como la paradoja de la velocidad de propagación. En este trabajo estudiaremos estefenómeno y plantearemos una solución, apoyándonos en un ejemplo concreto.

2 . Problema concreto

Consideramos el siguiente problema: 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. En el extremo derecho se consigue mantener la temperatura a 10°C mientras que en el izquierdo la temperatura es siempre de 1°C. Además, la temperatura en el instante inicial viene dada por la función [math] u_0 (x)= 10-10 \cdot 1_{[1/3,2/3]}(x)[/math] . Así, el sistema que modeliza este problema es el siguiente:

[math] \begin{cases} u_t-u_{xx} = 0 & x \in [0,1], t\gt0 \\ u(0,t)=10 & t\gt0 \\ u(1,t)=1 & t\gt0 \\ u(x,0)=u_0(x) & x \in [0,1] \end{cases} [/math]

Su solución estacionaria es:

[math] v(x)=9x+1 [/math]

y la solución general es:

[math] u(x,t)= 9x+1+ \sum_{n=1}^{\infty}b_n sen(n \pi x)e^{n^2\pi ^2 t}[/math]
,

donde [math] b_n [/math] es una constante de la forma:

AQUI VEMOS QUE HACEMOS CON BN

Si visualizamos estas soluciones con [math] t \in [0,1] [/math] y [math] x \in [0,1] [/math]:

Solución de la ecuación del calor

Observamos que la solución u comienza en la temperatura correcta (salvo por oscilaciones) y para un tiempo suficientemente pequeño se aproxima considerablemente a la solución de equilibrio. Además, en el siguiente apartado se comentará un detalle de la solución que se debe estudiar en mayor profundidad.

3 . Paradoja de la velocidad de propagación

Observamos que, en nuestra solución, para [math]t\gt0[/math] la solución es estrictamente positiva. Sin embargo, en [math]t=0[/math] observamos que para el intervalo [math] x \in [1/3,2/3][/math] la temperatura es cero. Esto contradice la teoría de la relatividad de Einstein, ya que podemos tomar un valor [math]t\gt0[/math] arbitrariamente pequeño de manera que la temperatura sea no nula, lo cual no es posible bajo la teoría de la relatividad. Este fenómeno se conoce como la paradoja de la velocidad de propagación del calor.

3.1 . Ecuación de Cattaneo-Vernotte

Para solucionar esta aparente falla en nuestro modelo, Cattaneo propuso una solución modificando la Ley de Fourier de la siguiente manera:

[math]q=-k\nabla u \rightarrow q+ \tau \frac{\partial q}{\partial t} =k\nabla u[/math]

Añadimos el factor [math] \tau \frac{\partial q}{\partial t} [/math] a la Ley de Fourier, donde [math]\tau[/math] se conoce como el término de relajación térmica. De esta manera, se modela un retardo en la respuesta de flujo de calor a los cambios de la temperatura. Si ahora aplicamos la conservación de la energía, nos queda la conocida ecuación hiperbólica:

[math] \frac{\partial u}{\partial t} + \tau \frac{\partial^{2}u}{\partial t^{2}} = \alpha \Delta u [/math]

donde [math] \alpha [/math] es la difusividad térmica.

Esta nueva ecuación garantiza que la velocidad de propagación del calor se produzca con velocidad finita, por lo que se utiliza en modelos donde el efecto de propagación instantánea tiene importancia.

Si resolvemos esta ecuación por separación de variables, tenemos las siguientes ecuaciones para las función espacial [math]X(x)[/math] y temporal [math] T(t)[/math]:

[math] X'' + \frac{\lambda}{\alpha}X=0, \hspace{5mm} X(0)=X(L)=0 [/math]
[math] \tau T'' + T' +\lambda T =0, \hspace{5mm} [/math]

La ecuación espacial tiene como solución general [math] X(x) = Asen(kx) + Bcos(kx) [/math], con [math] k^{2}=\frac{\lambda}{\alpha} \ltmath\gt. Usando las condiciones de Dirichlet, tenemos que \ltmath\gtB=0[/math] y nos queda la colección de soluciones

[math] X_{n}(x)=A_{n}sen(\frac{n\pi}{L}x) [/math]

4 . Referencias

  • Amin Moosaie. Non-Fourier heat conduction in a finite medium with insulated boundaries and arbitrary initial conditions ([1]).
  • Francisco R. Villatoro. La ciencia de la mula Francis: La velocidad de la propagación del calor, entre la paradoja y la entropía([2])