MateWiki - Contribuciones del usuario [es]

Ecuación del calor (ADMR)

2025-03-18T20:14:44Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Consideramos los parámetros, donde <math> n </math> es el término de la serie de Fourier donde truncamos:

*<math> a = b = 20 </math>
*<math> t_{max} = 1 </math>
*<math> n = 100 </math>

A continuación, ampliaremos <math>a, b</math> y <math>n</math> dejando fijo <math>t_{max}</math> y luego fijaremos <math>a</math> y <math>b</math> en vez del tiempo.
Podemos observar que al aumentar <math>a, b</math> y <math>n</math> el error entre las solucionnes acotadas y no acotada disminuye. Nótese que se estanca a partir de cierto <math>n</math>. De hecho, el error llega a aumentar a partir de ciertos <math> a,b </math> si truncamos demasiado pronto o demasiado tarde por la aproximación del trapecio para particiones del espacio demasiado bastas -no observado por pericia en código final-.

Al ampliar el tiempo con la <math>n</math> como los errores se concentran para tiempos pequeños, al hacer más grande <math>t_{max}</math> los errores máximos no aumentan. En cambio, si tomamos el máximo error para diferentes tiempos vemos con suficiente precisión no hay diferencia pero que para truncamientos precoces Neumann tiene ligeramente menor error.

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T19:13:13Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Consideramos los parámetros, donde <math> n </math> es el término de la serie de Fourier donde truncamos:

*<math> a = b = 20 </math>
*<math> t_{max} = 1 </math>
*<math> n = 100 </math>

A continuación, ampliaremos <math>a, b</math> y <math>n</math> dejando fijo <math>t_{max}</math> y luego fijaremos <math>a</math> y <math>b</math> en vez del tiempo.
Podemos observar que al aumentar <math>a, b</math> y <math>n</math> el error entre las solucionnes acotadas y no acotada disminuye. Nótese que se estanca a partir de cierto <math>n</math>. De hecho, el error llega a aumentar a partir de cierts <math> a,b </math> si truncamos demasiado pronto.

Al ampliar el tiempo con la <math>n</math> como los errores se concentran para tiempos pequeños, al hacer más grande <math>t_{max}</math> los errores máximos no aumentan. En cambio, si tomamos el máximo error para tiempos cercanos a <math>t_{max}</math> vemos que es mejor la Dirichlet.

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:24:57Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Consideramos los parámetros, donde <math> n </math> es el término de la serie de Fourier donde truncamos:

*<math> a = b = 20 </math>
*<math> t_{max} = 1 </math>
*<math> n = 100 </math>

A continuación, ampliaremos <math>a</math> y <math>b</math> dejando fijo <math>t_{max}</math> y luego ampliaremos al revés.

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:08:58Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Se consideran los parámetros:

*<math> a = b = 20 </math>
*<math> t_{max} = 1 </math>
*<math> n = 100 </math>

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:08:46Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Se consideran los parámetros:

*<math> a = b = 20 </math>
*<math> t_max = 1 </math>
*<math> n = 100 </math>

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:08:30Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Se consideran los parámetros:

*<math> \item a = b = 20 </math>
*<math> \item t_max = 1 </math>
*<math> \item n = 100 </math>

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:07:35Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Se consideran los parámetros:

<math>\begin{itemize}
\item a = b = 20
\item t_max = 1
\item n = 100
\end{itemize}</math>

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:06:09Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

Se consideran los parámetros:

*<math> a = b = 20
*t_max = 1
* n = 100
</math>

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:04:00Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>) y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T18:03:40Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| Representación gráfica de las soluciones del problema acotado (<math>a = b = 20</math>)y no acotado truncando las series de Fourier en su término <math>100</math>, para un tiempo máximo de <math>1</math>; y de sus errores. Nótese el orden de magnitud de los errores -pequeño-, los picos generados por los errores numéricos dados por la fórmula del trapecio y la suavidad en la frontera dada por las condiciones en esta.]]

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T17:57:00Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|550px|thumb|right| blabla]]

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T17:56:28Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|700px|thumb|right| blabla]]

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-18T17:56:14Z

Rafael Pascual: /* Análisis de los errores */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
[[Archivo:GrupoADMR Calor Errores1.jpg|200px|thumb|right| blabla]]

=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Archivo:GrupoADMR Calor Errores1.jpg

2025-03-18T17:54:11Z

Rafael Pascual:

Ecuación del calor (ADMR)

2025-03-18T17:31:44Z

Rafael Pascual:

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

El desarrollo de la ecuación del calor se atribuye principalmente al matemático y físico francés Joseph Fourier, quien en 1822 publicó su obra ''Théorie analytique de la chaleur''. En este trabajo formuló una ecuación diferencial en derivadas parciales que describe la propagación del calor por difusión en un medio continuo. Su análisis no sólo permitió comprender mejor este fenómeno físico, sino que también promovió el desarrollo de las Series de Fourier, que fueron objeto de estudio en nuestro trabajo anterior. En esta ocasión, profundicemos en dos aspectos de este modelo.

Desde un punto de vista físico, estas ecuaciones permiten la modelización de fenómenos como el calentamiento de materiales o la difusión de solutos en fluidos estacionarios, como ya hemos explorado en algunos ejercicios. Sin embargo, la influencia de la física es más profunda: de hecho, el coeficiente de difusión lleva su nombre por la interpretación física de su acción en el sistema. ¿Cuál es este efecto, y por qué se lo denota de esta forma?

Al tiempo que exploramos este efecto, planteemos también el paso entre problemas de dominio acotado y los de no acotado. Parecería evidente que estos últimos se extiendan naturalmente con los anteriores, pero en el desarrollo teórico de clase se obtenían mediante soluciones autosimilares en vez de problemas de autofunciones. ¿Convergen las soluciones del sistema acotado a la del no acotado cuando el dominio crece lo suficiente? Y de ser así, ¿qué condiciones de frontera en una dimensión funcionan mejor?

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

=Solución acotada vs Solución no acotada=
Demos por conocidos el uso de separación de variables para resolver problemas de autofunciones en caso acotado y el manejo de la solución fundamental junto con la convolución para el sistema no acotado. Para comparar todas las soluciones correspondientes, consideremos sistema <math> 1-</math>dimensional, con una gaussiana <math> e^{-x^2} </math> como condición inicial. La justificación de una curva así puede ser, en un experimento físico, el uso de un láser para calentar un punto de una tira fina de metal. Considerando que los fotones puedan desviarse como una distribución normal por el teorema central del límite, es un montaje experimental factible. Además, es una función continua y diferenciable, par y asintóticamente nula, buenas propiedades para un estado inicial del sistema. Pongamos ahora dos posibles condiciones en las fronteras.

Si consideramos de Dirichlet, a una tira de longitud <math> 2a </math> le colocaremos en los extremos un material que se mantenga siempre a la misma temperatura, pongamos que a una temperatura nula. El sistema resulta

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0, x\in [-a,a]\\
u(-a,t) = u(a,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2} \quad &\forall x\in [-a,a]
\end{cases}
\end{align}
</math>

La solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, con un flujo nulo que físicamente viene de poner extremos aislantes. Para cierto <math> b > 0,</math> y no necesariamente igual que con <math> a </math>. Consideraremos la ecuación del calor

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t > 0 , x\in [-b,b]\\
u_x(-b,t) = u_x(b,t) = 0 \quad &\forall t > 0 \\
u(x,0) = e^{-x^2}. \quad &\forall x\in[-b,b]
\end{cases}
\end{align}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, obteniendo

<math>\quad

u(x,t) = \sum_{m=0}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Un último detalle antes de proceder a los errores es la obtención simbólica de la solución no acotada, realizada a mano y mediante Maple.

<math>\quad
\begin{align}
\begin{cases} u_{t} - u_{xx} = 0 \quad &\forall t>0, x\in \mathbb{R}\\
u(x,0) = e^{-x^2} \quad &\forall x\in \mathbb{R}.
\end{cases}
\end{align}
\Longrightarrow
u(x,t) = \frac{1}{\sqrt{4\pi t}}e^{\frac{-x^2}{4t}} * e^{-x^2} = \frac{e^{\frac{-x^2}{4t+1}}}{\sqrt{4t+1}}
</math>

Ahora que pasamos a los errores, es relevante indicar la naturaleza numérica de los cálculos y la acumulación de errores en estos. Con esto en mente, queremos comprobar cuán bien aproximan las soluciones del caso acotado Dirichlet y Neumann a la obtenida por la solución fundamental. Para ello, escojamos un tiempo máximo <math>t_{max} </math> hasta el cuál tomaremos el error en norma <math> L^{\infty} </math>. Queremos estudiar

* El efecto de aumentar el rango de definición con <math> a,b </math>.
* El efecto de aumentar el número de términos en las series de los sistemas acotados.
* La relación entre estos parámetros, especialmente consecuencias por error acumulado.

=== Análisis de los errores ===
=Coeficiente de difusión=

[[Archivo:Fundamental.gif|300px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Vamos a interpretar como influye exactamente el coeficiente de difusión en la solución de la ecuación del calor. Imaginemos por un segundo que tenemos varias barras (cada una de un material diferente) y queremos estudiar como se difunde el calor a través de cada una de ellas. Usando Matlab podemos calcular y representar fácilmente la solución del problema para diferentes valores del coeficiente de difusión, que denotaremos por <math/> D <math/>.

[[Archivo:Fundamental.gif|400px|thumb|right|Solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión. La representación gráfica parte del tiempo inicial <math/> 0.01, <math/> esperando una delta de Dirac en el instante inicial. Nótese que en este tiempo, a mayor coeficiente de difusión más se ha difundido el calor en el sistema, lo que se nota especialmente en el valor máximo observado y su disminución cuanto más crece el coeficiente.]]

Como se puede observar en la animación, fijado un tiempo, a medida que crece el valor del coeficiente de difusión el pico de la solución fundamental se va haciendo más cada vez más pequeño, es decir, a mayor coeficiente de difusión, menor es el tiempo que tarda en expandirse el calor. De hecho, podemos observar este hecho directamente de la expresión de la solución fundamental del calor. La solución de la ecuación del calor unidimensional en una barra de longitud $ L = 1 $, con coeficiente de difusión $ D $ y condición inicial dada por una función gaussiana $ f(x) = e^{-x^2} $, se expresa como una serie de Fourier:

<math>
u(x,t) = \sum_{n=1}^{\infty} A_n e^{-D (n\pi)^2 t} \sin(n\pi x),
</math>

[Solución que depende de D]- [Viendo la solución del apartado 3 me estoy dando cuenta que también depende de la longitud de la barra, si estuviésemos trabajando con el problema acotado podríamos ver como afectan estos dos parámetros al tiempo de difusión]

Tomando el límite de D tendiendo a <math/> \infty <math/> vemos que u(x,t) tiende a 0, fijado un t determinado. Por tanto, podemos concluir que la solución del calor alcanzará el estado estacionario antes si su coeficiente del calor es mayor.

= Códigos =

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la solución fundamental de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,2,3,4,5,6,7,8,9,10];

% Nombre del archivo GIF de salida
output_gif = 'Fundamental.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for i = 1:length(DD)
D = DD(i);

% Función fundamental de la ecuación del calor
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);

% Definición de la malla
xx = -1:0.01:1; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t
[X, T] = meshgrid(xx, tt); % Malla de coordenadas
PHI_D = phi_D(X, T); % Evaluación de la función

% Graficar la superficie
clf;
surf(X, T, PHI_D);
shading flat
shading interp
colormap('jet')
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)$', 'Interpreter', 'latex');
title(['Solución fundamental de la ecuación del calor, cuando D=', num2str(D)])
axis([-1, 1, 0, 1, 0, 3])
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if i == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

[[Archivo:Convoluciones_teorica.gif|500px|thumb|right|Soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.]]

<source lang="matlab">
%%% Código generado junto con Chat GPT para la representación de la soluciones de la ecuación del calor unidimensional para distintos valores del coeficiente de difusión.

clc
clear all
close all

% Vector de valores de D (coeficiente de difusión)
DD = [1,10,20,30,40,50,60,70,80,90,100];

% Función gaussiana
gaus = @(x) (exp(-x.^2));

% Nombre del archivo GIF de salida
output_gif = 'Convoluciones_teorica.gif';

% Crear la figura
figura = figure(1);
grid on;
view(3);

% Bucle para iterar sobre cada valor de D
for h = 1:length(DD)
D = DD(h);

% Funciones fundamentales y convolución
phi_D = @(x, t) exp(-x.^2./(4*D*t))./sqrt(4*D*pi*t);
convo = @(x,t) (exp(-x.^2./(4*D.*t+1))/sqrt(4*D.*t+1));

% Definición de la malla
xx = -10:0.1:10; % Rango de x
tt = 10^(-2):0.01:1; % Rango de t

% Inicialización de la matriz de convolución
conv = zeros(length(xx), length(tt));

% Cálculo de la convolución
for j = 1:length(tt)
auxt = tt(j);
for i = 1:length(xx)
auxx = xx(i);
val_conv(i,j) = convo(auxx, auxt);
end
end

% Graficar la superficie
clf;
hold on;
view(3);
surf(xx, tt, val_conv.');
xlabel('x');
ylabel('t');
zlabel('$\Phi_D(x,t)*N(x)$', 'Interpreter', 'latex');
shading flat;
shading interp;
colormap('jet');
title(['Soluciones en dominio no acotado, cuando D=', num2str(D)]);
axis([-10, 10, 0, 1, 0, 1]);
hold off;
drawnow;

% Captura de la imagen para el GIF
frame = getframe(figura);
img = frame2im(frame);
[imind, cm] = rgb2ind(img, 256);

% Escribir la imagen en el archivo GIF
if h == 1
imwrite(imind, cm, output_gif, 'gif', 'LoopCount', Inf, 'DelayTime', 0.7);
else
imwrite(imind, cm, output_gif, 'gif', 'WriteMode', 'append', 'DelayTime', 0.7);
end

end
</source>

= Referencias =

* Sandro Salsa, Gianmaria Verzini. Partial Differential Equations in Action From Modelling to Theory ([https://link.springer.com/book/10.1007/978-3-031-21853-8 [1<nowiki>]</nowiki>]).

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

Ecuación del calor (ADMR)

2025-03-16T17:11:17Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:09:11Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:06:34Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:04:27Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:04:07Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:03:12Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T17:02:36Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T16:52:21Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T16:40:26Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T16:32:51Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-16T16:32:26Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-15T17:55:24Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, con <math> b \in \mathbb{R}</math> y no necesariamente igual que con <math> a </math>. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Veamos si, para valores de <math> a \text{y} b </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores y las diferencias que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:54:32Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, con <math> b \in \mathbb{R}</math> y no necesariamente igual que con <math> a </math>. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Notemos que no necesariamente <math> a = b </math>.
Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:54:10Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables, con <math> b \in \mathbb{R}</math> no necesariamente igual que con <math> a </math>. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Notemos que no necesariamente <math> a = b </math>.
Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:51:35Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{b^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:51:03Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{b} x\right) e^{-\dfrac{m^2\pi^2}{a^2} t}\quad \text{ donde }
B_m = \frac{\int_{-b}^{b} u(x,0) \cos\left(\frac{m\pi}{b} x\right)dx}{\int_{-b}^{b} \cos^2\left(\frac{m\pi}{b}x\right)dx}.

</math>

Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:50:08Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-b,t) = u_x(b,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables. Obteniendo

<math>\quad

u(x,t) = \sum_{m=1}^\infty B_{m} cos\left(\frac{m\pi}{a} x\right) e^{-\dfrac{m^2\pi^2}{a^2} t}\quad \text{ donde }
B_m = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{m\pi}{a} x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{m\pi}{a}x\right)dx}.

</math>

Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:44:42Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-a,t) = u_x(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables. Obteniendo

<math>\quad

u(x,t) = \sum_{n=1}^\infty B_{n} cos\left(\frac{n\pi}{a} x\right) e^{-\dfrac{n^2\pi^2}{a^2} t}\quad \text{ donde }
B_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{n\pi}{a} x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{n\pi}{a}x\right)dx}.

</math>

Veamos si, para valores de <math> a </math> grandes, las soluciones previamente halladas distan mucho entre sí, es decir, vamos a estudiar los errores que hay entre ambas soluciones.

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:30:51Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-a,t) = u_x(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables. Obteniendo

<math>\quad

u(x,t) = \sum_{n=1}^\infty B_{n} cos\left(\frac{n\pi}{a} x\right) e^{-\dfrac{n^2\pi^2}{a^2} t}\quad \text{ donde }
B_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{n\pi}{a} x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{n\pi}{a}x\right)dx}.

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:30:02Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

Consideremos ahora condiciones de frontera Neumann, de forma que el sistema quede aislado. Para cierto <math> a \in \mathbb{R}</math> consideraremos la ecuación del calor

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u_x(-a,t) = u_x(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus.
\end{cases}
</math>

De nuevo podemos resolver el sistema, por ser homogéneo, mediante separación de variables.

<math>\quad

u(x,t) = \sum_{n=1}^\infty B_{n} cos\left(\frac{n\pi}{a} x\right) e^{-\dfrac{n^2\pi^2}{a^2} t}\quad \text{ donde }
B_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{n\pi}{a} x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{n\pi}{a}x\right)dx}.

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:29:55Z

Rafael Pascual: /* Códigos */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:20:40Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}.

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:20:15Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\quad \text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:19:35Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frontera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}
\text{ donde }
A_n = \frac{\int_{-a}^{a} u(x,0) \cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}{\int_{-a}^{a} \cos^2\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right)dx}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:15:08Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty A_{n} cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:14:09Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

u(x,t) = \sum_{n=1}^\infty cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:13:36Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

\frac{\pi}{a}(n - \frac{1}{2})

u(x,t) = \sum_{n=1}^\infty cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}\left(n - \frac{1}{2}\right)^2 t}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:12:41Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

\frac{\pi}{a}(n - \frac{1}{2})

u(x,t) = \sum_{n=1}^\infty cos\left(\frac{\pi}{a}(n - \frac{1}{2}) x\right) e^{-\dfrac{\pi^2}{a^2}(n - \frac{1}{2})^2t}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:10:50Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

\frac{\pi}{a}(n - \frac{1}{2})
u(x,t) = \sum_{n=1}^\infty cos(\frac{\pi}{a}(n - \frac{1}{2}) x) e^{-\frac{\pi^2}{a^2}(n - \frac{1}{2})^2t}

</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T17:10:33Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] |
*Ángel De Lucas Miranda,
*Daniel Rodríguez Calderón,
*Marcos Cabellos Hernández,
*Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\
u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \\
u(x,0) = Gaus
\end{cases}
</math>

cuya solución al resolver el problema de autofunciones viene dada por

<math>\quad

\frac{\pi}{a}(n - \frac{1}{2})
u(x,t) = \sum_{n=1}^\infty cos(\frac{\pi}{a}(n - \frac{1}{2}) x) e^{-\frac{\pi^2}{a^2}(n - \frac{1}{2})^2t}
}
</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T16:52:28Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Ángel De Lucas Miranda, Daniel Rodríguez Calderón, Marcos Cabellos Hernández, Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\ u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \end{cases}
</math>

cuya solución viene dada por

<math>\quad
u(x,t) =
</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T16:52:18Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

Ecuación del calor (ADMR)

2025-03-15T16:52:05Z

Rafael Pascual: /* To do list */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Ángel De Lucas Miranda, Daniel Rodríguez Calderón, Marcos Cabellos Hernández, Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. y error numérico para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\ u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \end{cases}
</math>
cuya solución viene dada por

<math>\quad
u(x,t) =
</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T16:51:31Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Ángel De Lucas Miranda, Daniel Rodríguez Calderón, Marcos Cabellos Hernández, Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\ u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \end{cases}
</math>
cuya solución viene dada por

<math>\quad
u(x,t) =
</math>

=Códigos=

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

Ecuación del calor (ADMR)

2025-03-15T16:35:18Z

Rafael Pascual: /* Solución acotada vs Solución no acotada */

{{ TrabajoED | Ecuación del calor (Grupo ADMR). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Ángel De Lucas Miranda, Daniel Rodríguez Calderón, Marcos Cabellos Hernández, Rafael Pascual Ortega.}}

=Introducción=

=To do list=

1) (Problema acot y sol)Damos la sol. con frontera 0 --> para extender mejor a no acot
2) Dibujinchis
3) (Problema no acot) Dar el problema, sol, dibujinchi
4) Comparar ambas sol. error cuadrático, y el unif. para intervalos crecientes [-a,a]

5) Pasar de Diric. a Neum. y dar la explicacion fis. de Neum. Dibujinchis comparación, solapadas
6) Ver cómo afecta la cte calorifica, densidad y k = relacion flujo gradiente (la D) en ambos casos.

7) Introducción con cosas reales y prácticas

=Solución acotada vs Solución no acotada=

Vamos a considerar el siguiente problema de calor con condiciones Dirichlet cero en la frotera con cierto <math> a \in \mathbb{R} </math>

<math>\quad
\begin{cases} u_{t} - u_{xx} = 0 \\ u(-a,t) = u(a,t) = 0 \quad \forall t > 0 \end{cases}
</math>

=Códigos=

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