Diferencia entre revisiones de «Ecuación del calor (Grupo ACIRV)»
(→Modelización de la ecuación del calor con una dimensión) |
|||
| Línea 4: | Línea 4: | ||
=Introducción= | =Introducción= | ||
| + | En esta página ilustraremos el uso del método de diferencias finitas para resolver la ecuación del calor en una dimensión. | ||
| + | |||
| + | Dicho método sirve para resolver de manera aproximada ecuaciones diferenciales parciales y ordinarias de manera numérica. Funciona sustituyendo las derivadas por cocientes de diferencias que vienen dados por valores de la función en puntos discretos de una malla. | ||
| + | |||
| + | La ecuación a resolver es | ||
| + | |||
| + | <center><math> | ||
| + | u_t = \alpha u_{xx}, \quad x \in [0,1], \ t > 0,\ \alpha = 1. | ||
| + | </math></center> | ||
| + | |||
| + | |||
| + | donde \u(x,t)\ representa la temperatura en el punto $x$ y en el instante $t$, y $\alpha$ es una constante que viene dada por la conductividad térmica del material. | ||
=Modelización de la ecuación del calor con una dimensión= | =Modelización de la ecuación del calor con una dimensión= | ||
Revisión del 12:45 17 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
En esta página ilustraremos el uso del método de diferencias finitas para resolver la ecuación del calor en una dimensión.
Dicho método sirve para resolver de manera aproximada ecuaciones diferenciales parciales y ordinarias de manera numérica. Funciona sustituyendo las derivadas por cocientes de diferencias que vienen dados por valores de la función en puntos discretos de una malla.
La ecuación a resolver es
donde \u(x,t)\ representa la temperatura en el punto $x$ y en el instante $t$, y $\alpha$ es una constante que viene dada por la conductividad térmica del material.
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:
