Diferencia entre revisiones de «Ecuación del calor (PPAD)»
| Línea 1: | Línea 1: | ||
{{ TrabajoED | Ecuación del calor (Grupo PPAD). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Pablo Vidal Nacle, Pablo Maestro Fernández, Alex Heredero Santamaría, Diego Moñino Vizmanos}} | {{ TrabajoED | Ecuación del calor (Grupo PPAD). | [[:Categoría:EDP|EDP]]|[[:Categoría:EDP24/25|2024-25]] | Pablo Vidal Nacle, Pablo Maestro Fernández, Alex Heredero Santamaría, Diego Moñino Vizmanos}} | ||
| + | = Problema = | ||
| + | |||
| + | Se tiene la ecuación de difusión: | ||
| + | |||
| + | <center><math> u_t - u_{xx} = 0, \quad x \in (0,1), \quad t > 0 </math></center> | ||
| + | |||
| + | con las siguientes condiciones de frontera y condiciones iniciales: | ||
| + | |||
| + | <center><math> u(0,t) = 1, \quad t > 0 </math></center> | ||
| + | <center><math> u(1,t) = 10, \quad t > 0 </math></center> | ||
| + | <center><math> u(x,0) = 10, \quad x \in (0,1) </math></center> | ||
| + | |||
| + | = Cálculo de la solución estacionaria = | ||
| + | |||
| + | Para encontrar la solución estacionaria <math> v(x) </math>, se asume que cuando <math> t \to \infty </math>, la solución <math> u(x,t) </math> se estabiliza y tiende a una función <math> v(x) </math> que satisface: | ||
| + | |||
| + | <center><math> v''(x) = 0 </math></center> | ||
| + | |||
| + | con las condiciones de frontera: | ||
| + | |||
| + | <center><math> v(0) = 1 </math></center> | ||
| + | <center><math> v(1) = 10 </math></center> | ||
| + | |||
| + | Resolviendo la ecuación diferencial, tomamos <math> v(x) = ax + b </math>. Sustituyendo las condiciones de frontera: | ||
| + | |||
| + | <center><math> v(0) = b = 1 </math></center> | ||
| + | <center><math> v(1) = a + 1 = 10 \Rightarrow a = 9 </math></center> | ||
| + | |||
| + | Por lo tanto, la solución estacionaria es: | ||
| + | |||
| + | <center><math> v(x) = 9x + 1, \quad x \in (0,1) </math></center> | ||
| + | |||
| + | Esta función representa el estado estable de la temperatura en la barra cuando el tiempo tiende a infinito. | ||
| + | |||
| + | = Resolviendo la ecuación homogénea = | ||
| + | |||
| + | Se introduce el cambio de variable: | ||
| + | |||
| + | <center><math> w(x,t) = u(x,t) - v(x) </math></center> | ||
| + | |||
| + | De esta manera, <math> w(x,t) </math> satisface la ecuación de calor homogénea: | ||
| + | |||
| + | <center><math> w_t - w_{xx} = 0, \quad x \in (0,1), \quad t > 0 </math></center> | ||
| + | |||
| + | Las condiciones de frontera e iniciales se transforman en: | ||
| + | |||
| + | <center><math> w(0,t) = 0, \quad t > 0 </math></center> | ||
| + | <center><math> w(1,t) = 0, \quad t > 0 </math></center> | ||
| + | <center><math> w(x,0) = u(x,0) - v(x) = 10 - (9x + 1) = -9x + 9, \quad x \in (0,1) </math></center> | ||
| + | |||
| + | = Resolviendo por separación de variables = | ||
| + | |||
| + | Aplicando separación de variables y la serie de Fourier, se obtiene la solución para <math> w(x,t) </math>: | ||
| + | |||
| + | <center><math> w(x,t) = \sum_{n=1}^{\infty} B_n e^{-n^2 \pi^2 t} \sin(n\pi x) </math></center> | ||
| + | |||
| + | Donde los coeficientes de Fourier son: | ||
| + | |||
| + | <center><math> B_n = \frac{18}{n\pi} </math></center> | ||
| + | |||
| + | Finalmente, la solución general para <math> u(x,t) </math> es: | ||
| + | |||
| + | <center><math> u(x,t) = 9x + 1 + \sum_{n=1}^{\infty} \frac{18}{n\pi} e^{-n^2 \pi^2 t} \sin(n\pi x) </math></center> | ||
| + | |||
| + | = Nota sobre la función inicial = | ||
| + | |||
| + | Para calcular la serie de Fourier, se considera la extensión impar de la función inicial: | ||
| + | |||
| + | <center><math> | ||
| + | \tilde{g}(x) = | ||
| + | \begin{cases} | ||
| + | -9x + 9, & x \in (0,1) \\ | ||
| + | -9x - 9, & x \in (-1,0) | ||
| + | \end{cases} | ||
| + | </math></center> | ||
| + | |||
| + | Dado que <math> \tilde{g}(x) </math> es impar, su desarrollo de Fourier se mantiene en <math> x=0 </math>, asegurando la convergencia de la solución. | ||
| + | |||
| + | Por lo tanto, la solución final de la ecuación de calor con las condiciones dadas es: | ||
| + | |||
| + | <center><math> u(x,t) = 9x + 1 + \sum_{n=1}^{\infty} \frac{18}{n\pi} e^{-n^2 \pi^2 t} \sin(n\pi x) </math></center> | ||
Revisión del 13:33 15 mar 2025
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Ecuación del calor (Grupo PPAD). |
| Asignatura | EDP |
| Curso | 2024-25 |
| Autores | Pablo Vidal Nacle, Pablo Maestro Fernández, Alex Heredero Santamaría, Diego Moñino Vizmanos |
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Problema
Se tiene la ecuación de difusión:
con las siguientes condiciones de frontera y condiciones iniciales:
2 Cálculo de la solución estacionaria
Para encontrar la solución estacionaria [math] v(x) [/math], se asume que cuando [math] t \to \infty [/math], la solución [math] u(x,t) [/math] se estabiliza y tiende a una función [math] v(x) [/math] que satisface:
con las condiciones de frontera:
Resolviendo la ecuación diferencial, tomamos [math] v(x) = ax + b [/math]. Sustituyendo las condiciones de frontera:
Por lo tanto, la solución estacionaria es:
Esta función representa el estado estable de la temperatura en la barra cuando el tiempo tiende a infinito.
3 Resolviendo la ecuación homogénea
Se introduce el cambio de variable:
De esta manera, [math] w(x,t) [/math] satisface la ecuación de calor homogénea:
Las condiciones de frontera e iniciales se transforman en:
4 Resolviendo por separación de variables
Aplicando separación de variables y la serie de Fourier, se obtiene la solución para [math] w(x,t) [/math]:
Donde los coeficientes de Fourier son:
Finalmente, la solución general para [math] u(x,t) [/math] es:
5 Nota sobre la función inicial
Para calcular la serie de Fourier, se considera la extensión impar de la función inicial:
Dado que [math] \tilde{g}(x) [/math] es impar, su desarrollo de Fourier se mantiene en [math] x=0 [/math], asegurando la convergencia de la solución.
Por lo tanto, la solución final de la ecuación de calor con las condiciones dadas es: