0% encontró este documento útil (0 votos)
43 vistas4 páginas

Entregable 3. Solución de Una Ecuación de Calor 18/08/2025: Brian Izquierdo Cadena Viridiana Jimenez Martinez

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
43 vistas4 páginas

Entregable 3. Solución de Una Ecuación de Calor 18/08/2025: Brian Izquierdo Cadena Viridiana Jimenez Martinez

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

.

BRIAN IZQUIERDO CADENA

VIRIDIANA JIMENEZ MARTINEZ

Entregable 3. Solución de una


ecuación de calor
18/08/2025
Ecuación de calor en una barra 1D — Resolución
paso a paso
1) Datos del problema
Ecuación del calor: ∂u/∂t = α·∂²u/∂x², con α = 0.01.

Dominio: 0 ≤ x ≤ L, L = 10 cm.

Condiciones de frontera: u(0,t) = 0, u(L,t) = 0.

Condición inicial: u(x,0) = 100·sin(πx/L).

2) Separación de variables
Supóngase u(x,t) = X(x)·T(t). Entonces:

u_t = X·T′, u_xx = X″·T.

Sustituyendo en la PDE: X·T′ = α·X″·T ⇒ T′/(αT) = X″/X = −λ (constante de separación).

3) Problema propio espacial


ODE: X″ + λX = 0, con X(0)=0, X(L)=0.

Soluciones no triviales requieren λ>0. Escriba λ = k². Entonces: X(x) = A·sin(kx) + B·cos(kx).

De X(0)=0 ⇒ B=0. De X(L)=0 ⇒ sin(kL)=0 ⇒ kL = nπ, n=1,2,3,…

⇒ k_n = nπ/L, λ_n = (nπ/L)², X_n(x) = sin(nπx/L).

4) Parte temporal
Para cada n: T_n′ + αλ_n T_n = 0 ⇒ T_n(t) = C_n·exp[−αλ_n t] = C_n·exp[−α(nπ/L)² t].

5) Combinación y forma general


u(x,t) = Σ_{n=1}^∞ C_n·exp[−α(nπ/L)² t]·sin(nπx/L).

6) Determinación de coeficientes usando la condición inicial


A t=0: u(x,0) = Σ_{n=1}^∞ C_n·sin(nπx/L) = 100·sin(πx/L).

Use ortogonalidad en [0,L]: ∫₀ᴸ sin(mπx/L)·sin(nπx/L) dx = (L/2)·δ_{mn}.

Multiplique por sin(mπx/L) e integre:

∫₀ᴸ u(x,0)·sin(mπx/L) dx = Σ C_n ∫₀ᴸ sin(nπx/L)·sin(mπx/L) dx = C_m·(L/2).

Por tanto: C_m = (2/L) ∫₀ᴸ u(x,0)·sin(mπx/L) dx.

Con u(x,0)=100·sin(πx/L):
C_m = (2/L) ∫₀ᴸ 100·sin(πx/L)·sin(mπx/L) dx = 100·δ_{m1}.

⇒ C_1 = 100, C_{n≥2}=0. Renombrando C_n ≡ A_n: A_1=100, A_{n≥2}=0.

7) Solución final
u(x,t) = 100·exp[−α(π/L)² t]·sin(πx/L).

Con L=10 y α=0.01: u(x,t) = 100·exp[−(π²/10000) t]·sin(πx/10).

8) Parámetros característicos y gráfica


Tiempo característico: τ = L²/(απ²) ≈ 1013.212. Semivida: t_1/2 = τ·ln 2 ≈ 702.305.

Factores de atenuación exp[−α(π/L)² t]: t=500→0.610, t=1000→0.373, t=2000→0.139.

9) Comentarios
La condición inicial excita únicamente el primer modo propio, por lo que toda la dinámica
es un decaimiento exponencial de ese modo. La forma espacial seno se conserva y la
amplitud disminuye hasta tender a 0 °C.
• Boyce, W. E., & DiPrima, R. C. (2017). Elementary Differential Equations and
Boundary Value Problems (11th ed.). Wiley.

• Zill, D. G. (2018). Ecuaciones diferenciales con problemas con valores en la frontera


(10.ª ed.). Cengage Learning.

• Haberman, R. (2012). Applied Partial Differential Equations with Fourier Series and
Boundary Value Problems (5th ed.). Pearson.

• Strauss, W. A. (2007). Partial Differential Equations: An Introduction (2nd ed.). Wiley.

• Kreyszig, E. (2011). Advanced Engineering Mathematics (10th ed.). Wiley.

También podría gustarte