Diferencia entre revisiones de «Ecuación del calor (Grupo ACIRV)»
(→Modelización de la ecuación del calor con una dimensión) |
(→Modelización de la ecuación del calor con una dimensión) |
||
| Línea 27: | Línea 27: | ||
Puesto que no es un sistema homogéneo, debemos encontrar primero la solución estacionaria. Es decir, suponemos que \(t \rightarrow \infty\) (la solución ya no varía en el tiempo). | Puesto que no es un sistema homogéneo, debemos encontrar primero la solución estacionaria. Es decir, suponemos que \(t \rightarrow \infty\) (la solución ya no varía en el tiempo). | ||
| − | La solución estacionaria obtenida es | + | La solución estacionaria obtenida es \(v(x) = 9x +1\) . Al representarla gráficamente en 3D obtenemos: |
[[Archivo:Sol estacionaria.png|600px|thumb|right|Solución estacionaria]] | [[Archivo:Sol estacionaria.png|600px|thumb|right|Solución estacionaria]] | ||
Revisión del 19:52 16 mar 2025
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Ecuación del calor (Grupo ACIRV). |
| Asignatura | EDP |
| Curso | 2024-25 |
| Autores | Ángela Sotelo Fernández, Carmen Doñoro Molina, Inés Torres Gómez, Rubén Gutiérrez Hernández, Violeta Luján Barrios. |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
1 Introducción
2 Modelización de la ecuación del calor con una dimensión
En primer lugar vamos a resolver el siguiente problema de manera analítica, para después comparar su solución con la obtenida numéricamente.
Se considera una varilla metálica que ocupa el intervalo \([0,1]\) y que se encuentra aislada por su superficie lateral, de manera que la conducción de calor solo se produce en la dirección longitudinal. La temperatura inicial de la varilla es \(10^{\circ}C\). En el extremo derecho se consigue mantener la temperatura a \(10^{\circ}C\) mientras que en el izquierdo la temperatura es siempre de \(1^{\circ}C\).
El sistema que modeliza el comportamiento de la temperatura, representada por la función \(u(x,t)\), es el siguiente:
Puesto que no es un sistema homogéneo, debemos encontrar primero la solución estacionaria. Es decir, suponemos que \(t \rightarrow \infty\) (la solución ya no varía en el tiempo).
La solución estacionaria obtenida es \(v(x) = 9x +1\) . Al representarla gráficamente en 3D obtenemos:
%Creamos una matriz que representa una malla de puntos (tiempo,
% espacio) y en cada una de las columnas introducimos el valor de 9x+1 para todos los tiempos.
% Matriz
evaluaciones = zeros(100,100);
% Límite temporal y vectores
T = 1;
t = linspace(0,T,100);
x = linspace(0,1,100);
% Rellenamos la matriz
for j = 1:100
evaluaciones(:,j) = (9*x(j) + 1) * ones(100,1);
end
% Representación gráfica
surf(t,x,evaluaciones')
title('Solución estacionaria: v(x) = 9x + 1')
xlabel('Tiempo')
ylabel('Espacio')
zlabel('Temperatura (°C)')
Una vez calculada la solución estacionaria, homogeneizamos el sistema definiendo la ecuación \(w(x,t) = u(x,t) - v(x)\).
Resolviendo mediante el método de separación de variables y por superposición de soluciones, obtenemos la siguiente solución:
