Ver el código de «Ecuación del calor (PPAD)»

{{ TrabajoED | Ecuación del calor (Grupo PPAD). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Pablo Vidal Nacle, Pablo Maestro Fernández, Alex Heredero Santamaría, Diego Moñino Vizmanos}}

=. 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.

=. 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:

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

Su solución estacionaria es:

<center><math> v(x)=9x+1 </math></center>

y la solución general es:

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

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>:

[[Archivo:Sol11.jpeg|800px|marco|izquierda|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. 

=. Paradoja de la velocidad de propagación=
Observamos que, en nuestra solución, para <math>t>0</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>0</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. 

==. 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:

<center><math>q=-k\nabla u \rightarrow q+ \tau \frac{\partial q}{\partial t} =k\nabla u</math></center>

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:

<center><math> \frac{\partial u}{\partial t} + \tau \frac{\partial^{2}u}{\partial t^{2}} = \alpha \Delta u </math></center>

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>:

<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 =

* Amin Moosaie. Non-Fourier heat conduction in a finite medium with insulated boundaries and arbitrary initial conditions ([https://www.sciencedirect.com/science/article/pii/S0735193307001637 [1<nowiki>]</nowiki>]).

* Francisco R. Villatoro. La ciencia de la mula Francis: La velocidad de la propagación del calor, entre la paradoja y la entropía([https://francis.naukas.com/2008/10/22/la-velocidad-de-la-propagacion-del-calor-entre-la-paradoja-y-la-entropia/ [2<nowiki>]</nowiki>])

[[Categoría:EDP]]
[[Categoría:EDP24/25]]