1.
2 Linea con extremos a Temperaturas constantes $T_1$ y
$T_2$
Consideremos el PVI y de valores constantes $T_1$ y $T_2$ en los extremos de la linea
(barra). El problema para la distribución de temperatura $u(x,t)$ es
\begin{equation*} u_t=\alpha^2u_{xx} \end{equation*}\begin{equation*}
u(0,t)=T_1\hspace{2cm}u(L,t)=T_2\hspace{2cm}u(x,0)=f(x) \end{equation*}
Definido para $0<x<L$ y $t>0$. Se supone $\alpha>0$ constante
De la integración física del problema como una línea (aislada en su longitud $L$) que
disipa ó absorve calor por sus extremos, los cuales se encuentran a temperatura
constante $T_1$ y $T_2$ respectivamente. Se espera que despues de un tiempo
suficientemente largo $\tau$ la línea tenga una distribución de temperatura en equilibrio
la cual no depende del tiempo, es decir que se espera que $u(x,t)\longrightarrow V(x)$
para $t>\tau$
Si es así, entonces $u_t=\alpha^2u_{xx}$ se convierte en $V''(x)=0$ y las condiciones
de frontera $\quad u(0,t)=T_1=V(0)$, $\quad u(L,t)=T_2=V(L)$
La solución general de la ecuación $V''(x)=0$ es $V(x)=C_1x+C_2$
Ahora imponiendo las condiciones iniciales se obtiene
\begin{equation*} T_1=V(0)=C_2\hspace{2cm}T_2=V(L)=C_1L+C_2 \end{equation*}
De donde
\begin{equation*} C_1=\frac{T_2-T_1}{L}\hspace{2cm}C_2=T_1 \end{equation*}
Por tanto la solución de equilibrio es
\begin{equation*} V(x)=\left(\frac{T_2-T_1}{L}\right)x+T_1 \end{equation*}
Por otra parte, para tiempos cercanos a cero, $(t\longrightarrow 0)$, la temperatura $
u(x,t)$ es una función que depende de $x$ y $t$, luego $\forall x \in (0,L)$ y todo $t>0$
\begin{equation*} u(x,t)=\omega(x,t)+V(x)\hspace{2cm}(1) \end{equation*}
donde $\omega(x,t)\longrightarrow0$ para $t>\tau$
Calculando las derivadas parciales e incluyendolas en la ecuación de calor se tiene
\begin{equation*} \omega_t=\alpha^2(\omega_{xx}+V'') \end{equation*}
Usando la forma $V(x)=\left(\dfrac{T_2-T_1}{L}\right)x+T_1$ para la distribución de
temperaturas en equilibrio se tiene que $V''(x)=0$ $\forall t>0$ por tanto la ecuación
anterior se reduce
\begin{equation*} \omega_t=\alpha^2\omega_{xx}\hspace{1cm}\forall t>0 \end{equation*}
Esto es, si la función (1) es solución de la ecuación $u_t=\alpha^2u_{xx}$ entonces
$\omega(x,t)$ debe satisfacer la ecuación de calor, luego reemplazando las condiciones
de frontera en (1) se tiene
\begin{align*}
u(0,t)&=T_1=\omega(0,t)+V(0)\hspace{0.5cm}\text{pero}\hspace{0.5cm}V(0)=T_1\longrightar
(2)\\
u(L,t)&=T_2=\omega(L,t)+V(L)\hspace{0.5cm}\text{pero}\hspace{0.5cm}V(L)=T_2\longrightar
(3) \end{align*}
Es decir que $u(x,t)=\omega(x,t)+V(x)$ satisface las condiciones de frontera siempre
que $\omega(x,t)$ satisfaga las condiciones de frontera homogéneas (2) y (3)
Por tanto, el problema de valores homogéneos en la frontera para la ecuación de calor
es el que se resolvió en la sección anterior, por tanto
\begin{equation*} \omega(x,t)=\sum_{n=1}^{\infty}b_n\sin{\left(\frac{n\pi}{L}x\right)}e^{-
\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t} \end{equation*}
Ahora como
\begin{equation*} f(x)=u(x,0)=T_1+\left(\frac{T_2-T_1}
{L}\right)x+\sum_{n=1}^{\infty}b_n\sin{\left(\frac{n\pi}{L}x\right)} \end{equation*}
entonces
\begin{equation*} f(x)-T_1-\left(\frac{T_2-T_1}
{L}\right)x=\sum_{n=1}^{\infty}b_n\sin{\left(\frac{n\pi}{L}x\right)} \end{equation*}
Multiplicando por $\sin{\left(\frac{m\pi}{L}x\right)}$ e integrando a ambos lados sobre
$[0,L]$ y haciendo un procedimiento an\'alogo a la secci\'on anterior se encuentra que
$b_n$ son los coeficientes de la serie de fourier de senos para la funci\'on $f(x)-T_1-
\left(\frac{T_2-T_1}{L}\right)x$, es decir
\begin{equation*} b_n=\frac{2}{L}\int_0^L\left[f(x)-T_1-\left(\frac{T_2-T_1}
{L}\right)x\right]\sin{\left(\frac{n\pi}{L}x\right)}dx \end{equation*}
Por tanto la solución al problema de valores iniciales y valores constantes en la frontera
para la ecuación de calor es
\begin{align*} u(x,t)&=T_1+\left(\frac{T_2-T_1}
{L}\right)x+\sum_{n=1}^{\infty}b_n\sin{\left(\frac{n\pi}{L}x\right)}e^{-
\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t}\\ u(x,t)&=T_1+\left(\frac{T_2-T_1}
{L}\right)x+\sum_{n=1}^{\infty}\left[\frac{2}{L}\int_0^L\left[f(x)-T_1-\left(\frac{T_2-T_1}
{L}\right)x\right]\sin{\left(\frac{n\pi}{L}x\right)}dx\right]\sin{\left(\frac{n\pi}{L}x\right)}e^{-
\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t} \end{align*}
1.3. Temperatura de una barra con extremos
aislados
Consideremos el problema de conducción de calor en una linea sin perdida de energía
en los extremos. El problema con valores en la frontera para la distribución de
temperatura $u(x,t)$ es
\begin{align*} \frac{\partial u}{\partial t}=\alpha^2\frac{\partial^2u}{\partial
x^2}\hspace{2cm}&0<x<L,\hspace{2cm}t>0\\ \frac{\partial u}{\partial x}
(0,t)=0\hspace{1cm},\hspace{1cm}&\frac{\partial u}{\partial x}(L,t)=0\hspace{2cm}t>0\\
&u(x,0)=f(x) \end{align*}
Este problema puede interpretarse físicamente como la difusión de calor en el interior de
una linea (barra) de longitud $L$ que se encuentra completamente aislada y que parte
con una distribución inicial de temperatura $f(x)$
Usando el método de separación de variables y haciendo un análisis análogo al de la
sección (1.1) se encuentra la solución
\begin{equation*} u(x,t)=\frac{a_0}{2}+\sum_{n=1}^{\infty}a_n\cos{\left(\frac{n\pi}
{L}x\right)}e^{-\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t} \end{equation*}
Con los coeficientes $a_n$ dados por
\begin{equation*} a_n=\frac{2}{L}\int_0^Lf(x)\cos{\left(\frac{n\pi}{L}x\right)}dx
\end{equation*}
Para todo $n=0,1,2,...$
En particular, el término constante \begin{equation*} \frac{a_0}{L}=\frac{1}
{L}\int_0^Lf(x)dx \end{equation*}
Es decir, el valor promedio de la distribución inicial de temperatura $f(x)$ sobre toda la
longitud de la linea
Se puede observar que
\begin{align*} \lim\limits_{t \rightarrow \infty}u(x,t)&=\lim\limits_{t \rightarrow
\infty}\left[\frac{a_0}{L}+\sum_{n=1}^{\infty}a_n\cos{\left(\frac{n\pi}{L}x\right)}e^{-
\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t}\right]\\ &=\frac{a_0}{L} \end{align*}
Esto es, la distribución de temperatura de equilibrio de la barra completamente aislada
es el valor promedio de su distribución inicial.
Ejemplo:
Si $f(x)=x(L-x)$
\begin{align*} \frac{a_0}{L}=\frac{1}{L}&\int_0^L\xi(L-\xi)d\xi=\frac{1}{L}\int_0^L(L\xi-
\xi^2)d\xi=\left.\frac{\xi^2}{2}-\frac{\xi^3}{3L}\right|_0^L=\frac{L^2}{2}-\frac{L^3}
{3L}=\frac{L^2}{2}-\frac{L^2}{3}=\frac{L^2}{6}\\ a_n=\frac{2}{L}&\int_0^L\xi(L-
\xi)\cos{\left(\frac{n\pi}{L}\xi\right)}d\xi=\frac{2}{L}\left[L\int_0^L\xi\cos{\left(\frac{n\pi}
{L}\xi\right)}d\xi-\int_0^L\xi^2\cos{\left(\frac{n\pi}{L}\xi\right)}d\xi\right]\\ =\frac{2}
{L}&\left[L\left[\frac{L}{n\pi}\xi\sin{\left(\frac{n\pi}{L}\xi\right)}+\frac{L^2}
{n^2\pi^2}\cos{\left(\frac{n\pi}{L}\xi\right)}\right|_0^L\right]-\\ &\left[\left.\xi^2\frac{L}
{n\pi}\sin{\left(\frac{n\pi}{L}\xi\right)}+2\xi\frac{L^2}{n^2\pi^2}\cos{\left(\frac{n\pi}
{L}\xi\right)}+2\frac{L^3}{n^3\pi^3}\sin{\left(\frac{n\pi}{L}\xi\right)}\right|_0^L\right]\\
=\frac{2}{L}&\left[L\left[0+\frac{L^2}{n^2\pi^2}\cos{(n\pi)}-\frac{L^2}{n^2\pi^2}\right]-
\left[0+\frac{2L^3}{n^2\pi^2}\cos{(n\pi)}-0\right]\right]\\ =\frac{2}{L}&\left[\frac{L^3}
{n^2\pi^2}\cos{(n\pi)}-\frac{L^3}{n^2\pi^2}-\frac{2L^3}{n^2\pi^2}\cos{(n\pi)}\right]\\ =\frac{2}
{L}&\left[-\frac{L^3}{n^2\pi^2}-\frac{L^3}{n^2\pi^2}\cos{(n\pi)}\right]=-\frac{2L^2}{n^2\pi^2}
[\cos{(n\pi)}+1]=-\frac{2L^2}{n^2\pi^2}\left[(-1)^n+1\right] \end{align*}
La solución para este caso es
\begin{equation*} u(x,t)=\frac{1}{6}L^2-2L^2\sum_{n=1}^{\infty}\frac{1}
{n^2\pi^2}\left[(-1)^n+1\right]\cos{\left(\frac{n\pi}{L}x\right)}e^{-\left(\frac{n^2\pi^2\alpha^2}
{L^2}\right)t} \end{equation*}
pero como
\begin{equation*} (-1)^n+1=\left\{\begin{array}{cc} 2 & \text{Si n es par} \\ 0 & \text{Si n es
impar} \end{array}\right. \end{equation*}
La solución puede escribirse
\begin{align*} u(x,t)&=\frac{1}{6}L^2-4L^2\sum_{n=1}^{\infty}\frac{1}
{(2n)^2\pi^2}\cos{\left(\frac{2n\pi}{L}x\right)}e^{-\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t}\\
&=\frac{1}{6}L^2-\frac{L^2}{\pi^2}\sum_{n=1}^{\infty}\frac{1}{n^2}\cos{\left(\frac{2n\pi}
{L}x\right)}e^{-\left(\frac{n^2\pi^2\alpha^2}{L^2}\right)t} \end{align*}
