Parte de Andrews y Lucía
Contenido
1 Introducción
En este documento se pretende mostrar al lector como la ecuación del calor en una dimensión describe el fujo de calor [math] u(x,t) [/math] ... Para ello estudiaremos distintas condiciones frontera e iniciales en una barra metálica que ocupa un intervalo [0,1]. CREO QUE AITANA TIENE QUE DECIR COSAS AQUI
2 Sistema no homogéneo
En este primer caso nos centraremos en una barra metálica que comienza estando a 0 °C y cuyas temperaturas al principio y al final de son constantes, pero distintas entre si. En concreto, consideraremos que la temperatura en la posición x = 0 es nula, y sin embargo, en x = 1 la sube un grado. Asimismo, estudiaremos la ecuación del calor cuya conductividad térmica, k, y calor específico se consideran 1. Todo esto se traduce en el sistema no homogéneo,
2.1 Solución Estacionaria
La resolución de este sistema se basa en el método de separación de variables perteneciente a la teoría de resolución ecuaciones diferenciales, lo cual carece de interés en este documento. Por limpieza en la lectura, este procedimiento no se incluirá, si embargo, hay un paso previo al método que cabe incluir. La resolución por separación de variables requiere que el sistema sea homogéneo. Esta modificación en el sistema original la conseguimos haciendo uso de la que se conoce como solución estacionaria. Esta se alcanza cuando ha pasado un tiempo infinito ([math] t \rightarrow \infty [/math]) y considerando por tanto, que el flujo del calor ha dejado de depender del tiempo, [math] u(t,x) \sim v(x) [/math]. Haciendo los respectivos cálculos es fácil llegar a que la solución estacionaria es [math] v(x) [/math] con [math] x\in[0,1] [/math]. La cual gráficamente muestra una bajada de temperatura en el espacio y tiempo.
2.2 Solución del sistema no homogéneo
Una vez hallada la solución estacionaria, el método de separación de variables nos devuelve el candidato a solución, [math] u_k(x,t)=x-c_k\cdot e^{-k^2\pi^2 t}\cdot \sin{(k\pi x)}[/math] Donde [math] c_{k} = \frac{2(-1)^k}{k\pi} [/math] son los coeficientes de Fourier hallados mediante aproximación impar. Al igual que con la solución estacionaria, podemos ver que la temperatura va decreciendo en el espacio y tiempo.
Más aún, es fácil observar como la tras el paso del tiempo, esta solución se va aproximando a la estacionaria. (NO SE SI PUEDO DECIR NADA MÁS).
2.3 Flujo del calor entrante y saliente
NI IDEA
2.4 ¿Qué pasa si la tomamos [math] k = \frac{1}{2} [/math]?
De momento en la ecuación del calor hemos considerado todas las constantes como 1 pero, ¿qué pasaría si la conductividad térmica disminuye?, ¿podremos apreciar algún cambio grande en la solución final?. Para ello tomamos [math] k = \frac{1}{2} [/math] y resolvemos el mismo sistema pero esta vez con ecuación,
Homogeneizar dicho sistema nos devuelve la solución estacionaria [math] v(x)_{\frac{1}{2}} = x[/math] ESCRIBIR SOLUCIÓN, que gráficamente tiene la forma,
DECIR ALGO DE LA FOTO De la misma forma que en el sistema con [math] k = 1 [/math], obtenemos la solución final, ESCRIBIR LA SOLUCIÓN FINAL La cual gráficamente se comporta de la siguiente manera
Comparando esta gráfica con la correspondiente de [math] u(x,t) [/math] no somos capaces a simple vista de encontrar ninguna diferencia. Una forma más visual de apreciar el efecto de la disminución de [math]k[/math] es comparando las soluciones de los sistemas con sus estacionarias correspondientes. Es decir, vamos a ver con cuanta "rapidez" cada sistema alcanzará su correspondiente solución estacionaria. Esto los haremos calculando las diferencias entre la solución y el estado estacionario para [math] x =\frac{1}{2} [/math] en ambos casos; [math] u(\frac{1}{2},t) = v(\frac{1}{2}) [/math] y [math] u_{\frac{1}{2}}(\frac{1}{2},t) = v_{\frac{1}{2}}(\frac{1}{2}) [/math].
Habiendo entendido el papel que juega la solución estacionaria resulta bastante intuitivo que esta se alcance en más tiempo en aquella barra con menor conductividad térmica. Esto es debido a que el calor pasará con mayor dificultad por la barra, ralentizando así la llegada a la solución estacionaria.
3 Sistema con cambios en la condición inicial
En este apartado vamos a estudiar el papel que juega la condición inicial imponiendo una condición que varíe con en el espacio. En concreto vamos a considerar la condición inicial [math] u(x,0) = \max\{0, 1-4·|x-1/2|\} - x [/math]. Para facilitar el estudio vamos a considerar directamente un sistema homogéneo, es decir, [math] u(0,t) = u(1,t) = 0 [/math]. El sistema completo sería,
Debido a que el sistema ya es homogéneo de partida, la solución estacionaria pasaría a ser [math] v(x) = 0 [/math] perdiendo así interés.
3.1 Solución del sistema
Tras aplicar el correspondiente método de separación de variables obtenemos la solución del sistema,
Donde los coeficientes de Fourier [math]c_k = 2\int_{0}^{1}\sin(k \pi x)u(x,0)[/math], se han calculado en Matlab aproximando las integrales por el método del trapecio,
clear
close all
format long
clc
% Intervalos con los que trabajaremos:
x=0:0.001:1; %longitud
t=0:0.001:1; %tiempo
% Cálculo de los coeficientes de Fourier:
K=10; % K primeros términos de la serie
f=@(x) max(0, 1-4*abs(x-1/2)); %condicion inicial
% Fórmula del trapecio para aproximación:
c= zeros(K,1); %matriz para almacenar los coeficientes de Fourier
for k=1:K
c(k)=2*trapz(x,f(x).*sin(k*pi*x));
end
% La solución del sistema de EDP viene dada por:
u=@(x,t)0;
for k=1:K
u=@(x,t) u(x,t)+c(k).*exp(-k^2*pi^2*t).*sin(k*pi*x);
end
% Creamos una malla de puntos:
[X,T]=meshgrid(x,t);
U=u(X,T);
% Representamos gráficamente:
figure
mesh(X,T,U)
xlabel('x')
ylabel('t')
zlabel('u(x,t)')
title('Solución de u(x,t)')Veamos a medida que pasa el tiempo nuestra solución se aproxima al estado estacionario [math]v(x) = 0[/math]. Otra cosa interesante a mencionar de esta gráfica es que la solución es estrictamente positiva si [math]t \gt 0[/math]. NO SE QUE DECIRLE AQUI
3.2 Flujo del calor entrante y saliente
NI IDEA
4 Sistema con cambios en las condiciones frontera
PARTE DE ANDREA
