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

Revisión del 19:15 5 mar 2024

1 Introducción

2 Ecuación del Calor en dimensión [math]n=1[/math]

2.1 Definición

Se considera una varilla metálica que ocupa el intervalo [math][0,L] [/math] 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. Inicialmente se planteará el problema de manera que la temperatura inicial de la varilla es [math]u1(x)[/math]. En el extremo izquierdo se establece una temperatura [math]u2(t)[/math] y en el extremos derecho [math]u3(t)[/math]. Dado esto, a partir del principio de conservación de la energía^[1], y la Ley de Fourier^[2], se plantea el problema de EDP que modeliza el comportamiento de la temperatura.

Donde la función [math]u(x,t)[/math] depende del espacio y del tiempo; donde [math]k[/math] es la conductividad térmica y [math]c[/math] el calor específico.

Este problema se ha planteado con condiciones frontera de Dirichlet.

2.2 Ejemplo

A continuación, se estudiará el comportamiento de este problema a partir de un caso particular.

Se considerará el intervalo [math][0,1] [/math], temperatura inicial de la varilla de 0º, y en los extremos, izquierdo y derecho, 0º y 1º respectivamente.

Además, se considerará que tanto la conductividad térmica, [math]k[/math], como el calor específico [math]c[/math] son 1. Entonces se plantea el siguiente problema a modelizar.

Una vez definido el problema, vamos a resolverlo para así observar su comportamiento.

2.2.1 Homogeneización y solución estacionaria

Para resolver esta EDP, en primer lugar se deben tener las condiciones frontera homogeneizadas, de manera que, [math] u(0,t)= 0, u(1,t)= 0 [/math]. Por tanto, en primer lugar debemos estudiar la solución estacionaria del problema. Para ello, supondremos que [math]{t \to \infty} [/math], en dicho caso, [math] u_{t} \approx 0[/math] y [math] u(x,t) \approx v(x)[/math].

Se obtiene entonces una EDO homogénea de orden 2, cuya solución es la siguiente.

A continuación, se muestra el código utilizado para representar la gráfica de la solución estacionaria, y la respectiva gráfica.

Representación de la solución estacionaria

%Definimos la solución estacionaria.
v=@(x) x;
%definimos el intervalo donde vamos a representar la solución. 
x=linspace(0,1,1000);

%Comando para representar la solución en R^2.
hold on 
plot(x, v(x), 'b-')
xlabel('x');
ylabel('v(x)');
title('Solución extacionaria.');
grid on; 
hold off

2.2.2 Separación de Variables

Posterioremente se utilizará el método de separación de variables^[3], a partir del cual, se obtiene la solución de a siguiente EDP, donde se ha realizado el cambio [math] w(x,t)=v(x) - u(x,t)[/math]

Una vez aplicado el método y obteniendo los coeficientes de Fourier correspondientes, tal y como se explicó en el anterior trabajo ^[4] para funciones impares, obtenemos la siguiente solución de [math] w(x,t)[/math].

Deshacemos el cambio de variable para obtener la solución del problema original, donde [math] u(x,t)=v(x) - w(x,t)[/math].

2.2.3 Representación y análisis

Una vez obtenida la solución, veremos la propagación del calor en función del tiempo y el espacio. Para ello, uilizaremos el siguiente programa que nos representará la solución obtenida tomando los 10 primeros términos de la serie.

Representación de la solución u(x,t)

% Definimos los términos de la serie
n=10;

%inicializamos la función en 0. 
f=@(x,t) 0;

%Tomamos un bucle para poder obetener la solución w con el nº de términos
%seleccionado.
for k=1:n
    f = @(x,t) exp(-k^2*pi^2*t).*sin(pi * k * x).*(2/(k*pi)).*((-1)^(k+1)) + f(x,t);
end

% Queremos representar u entonces definimos la solución deshaciendo el cambio
% w(x,t)=v(x)-u(x,t),donde x(x)=x.
g = @(x,t) x-f(x,t)

% Dibujar la solución en el intervalode [0,1].
[x,t] = meshgrid(0:.001:1);

%utilizamos el comando mesh para representar la solución en R^3.
mesh(x,t,g(x,t)) 
xlabel('x');
ylabel('t');
zlabel('u(x,t)');
title('Solución para k = 10');

Se puede observar como la representación de la solución final es muy similar a la que sería la representación en [math]R^{3} [/math] de la solución estacionaria, en el intervalo[math][0.2,1] [/math]. Fuera de este, la solución no se asemeja a la solución estacionaria. Debido a esto, debemos señalar que la solución obtenida del problema modificado, la [math]w(x,t) [/math], será similar al plano generado en [math]x=0[/math], fuera del entorno de [math]t=0 [/math], ya que como hemos visto, en el intervalo [math][0,0.2] [/math] la gráfica de [math]u(x,t) [/math] no se asemeja a la solución estacionaria.

Representación de la comparativa entre u(x,t) y v(x)

En este gráfica se puede observar lo anteriormente mencionado, donde la solución estacionaria y la original en gran parte del intervalo son prácticamente idénticas, formando prácticamente el mismo plano, tal y como se observa en la foto de la derecha.

2.2.4 Flujo

3 Ecuación del Calor en diferentes dimensiones

En este apartado, se verá que aspecto toma la función solución del calor en diferentes dimensiones y diferentes condiciones.

3.1 Solución fundamental de la ecuación del calor en dimensión 1

La solución fundamental del calor en una dimensión es la siguiente:

Se representará para [math] x \in [-1,1], t\in [10^{-2},1] [/math]

% Definimos el recinto en el que vamos a representar la función
[X, T] = meshgrid(-1:0.05:1, 10^(-2):0.05:1);
%Definimos la solución fundamental del calor en una dimensión
u= 1./sqrt(2*pi*T).*exp(-X.^(2)./(4*T));
% Graficamos la función
figure;
surf(X, T, u);
title('Gráfico de la función');
xlabel('X');
ylabel('T');
zlabel('u(x,t)');

Considerando las siguientes condiciones iniciales y frontera:

ecuaci´on del calor en una dimensi´on en el semiespacio x > 0 con la condici´on inicial u(x, 0) = 0 y la condici´on frontera u(0, t) = 1. ...

3.2 Solución de la ecuación del calor en diferentes situaciones

En el siguiente caso, se considera el siguiente problema:

La solución, viene dada por la convolución respecto la solución fundamental del calor y la condición inicial, es decir,

Se debe notar que la función anterior aplicada a la condición inicial, en este caso resulta

Para obtener una mayor visión de dicha solución, se representará en diferentes instantes de tiempo, [math] t \in (0.001, 0.01, 0.1) [/math] y en el espacio [math] x \in [-1,1] [/math]

Solución en diferentes instantes de tiempo

clc
clear 
format long
close all

t_vector = [0.001, 0.01, 0.1]; % Parámetros de t
X= -1:10^(-3):1; %Intervalo para representar
Y= -1:10^(-3):1; %Intervalo para integración
[fila,columna]= size(X);

%Hacemos el bucle para los 3 valores de t
for j=1:3
    t= t_vector(j);
    F = ones(size(X));
    for i=1:columna
        %Definimos cada elemento dado por la solución
        F(i) = F(i)* 1/sqrt(4 * pi * t) * trapz(Y,exp(-(X(i)-Y).^2 ./ (4*t)));
    end
    %Graficamos la función
    subplot(1,3,j);
    plot(X,F);
    xlabel('x');
    ylabel('u(x,t)');
    title(['t = ',num2str(t)]);
end

Por otro lado, la solución del calor en dimensión [math] 2[/math] es la siguiente:

De igual forma, se representará la solución para instantes de tiempo [math] t \in (0.001, 0.01, 0.1) [/math], y para [math] (x,y) \in [-1,1]\times[-1,1] [/math]

Solución en diferentes instantes de tiempo

clc
clear 
format long
close all

% Definimos el recinto en el que vamos a representar la función
[X1, X2] = meshgrid(-1:0.05:1, -1:0.05:1);
% Definimos el vector de tiempos
t_vector=[0.001,0.01,0.1];
%Preparamos la representación
tiledlayout(1,3);
%Hacemos un bucle para los tiempos
for i=1:3
t= t_vector(i);
%Definimos la solución fundamental del calor en dos dimensiones
u= 1./sqrt(2*pi*t).*exp((-X1.^(2)-X2.^(2))./(4*t));
%Graficamos la solución
nexttile
surf(X1, X2, u);
title(['t = ',num2str(t)]);
xlabel('X');
ylabel('X');
zlabel('u(x,t)');
end

En la gráfica se puede observas como la solución se hace más estrecha a menor [math] t[/math]. Esto se debe a que la solución no está definida para [math] t=0 [/math] y, por ello, tiende a estrecharse.

4 Referencias

1. ↑ Conservación de la energía
2. ↑ Ley de Fourier
3. ↑ Método de separación de variables
4. ↑ Series de Fourier