Contenido

Unidad 5: Ecuaciones diferenciales parciales

Ciclo: 2025 - I
Semana 16 - sesión 32

La ecuación de calor
Universidad Nacional Agraria La Molina Con condiciones de Neumann

Departamento Académico de Matemática

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 2 15

Logro de la sesión Barra con extremos aislados (Condiciones de Neumann)

Ocurre un problema ligeramente diferente si los extremos de la barra están aislados para que no
haya paso de calor a través de ellos.

∂u(0, t) ∂u(L, t)
=0 =0
∂x ∂x
Al finalizar la sesión, el estudiante:
• Resuelve problemas con valores iniciales y condiciones de fronteras para la ecuación de calor
unidimensional.
x
0 L

Usando el principio de conducción del calor: esto implica que el gradiente de temperatura debe
anularse en estos extremos; es decir, la tasa de flujo de calor a través de una sección transversal
es proporcional a la tasa de cambio de la temperatura en la dirección x. Ası́,
∂u ∂u
(0, t) = (L, t) = 0, t > 0.
∂x ∂x
Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 3 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 4 15
 En el siguiente ejemplo obtenemos la solución formal del problema de flujo de calor con estas
condiciones en la frontera.
Ejemplo 1 (Ecuación de calor en una varilla con extremos aislados) Veamos las condiciones de frontera

Determine una solución formal del problema de flujo de calor descrito por el problema con ∂u ∂u
(0, t) = X ′ (0) T (t) = 0 y (L, t) = X ′ (L) T (t) = 0.
valores iniciales y en la frontera (Neumann): ∂x ∂x
Si T (t) = 0 entonces u (x, t) = 0 es la solución trivial, pero estamos buscando soluciones no

 ut = βuxx , 0 < x < L, t > 0
ux (0, t) = ux (L, t) = 0, t > 0 triviales, ası́ que

u (x, 0) = f (x) , 0 < x < L X ′ (0) = X ′ (L) = 0
Luego, tenemos el problema con valores en la frontera:
Solución. Usamos el método de separación de variables para suponer que
 ′′
X (x) − λX (x) = 0
u (x, t) = X (x) T (t) . (2)
X ′ (0) = X ′ (L) = 0
Al sustituir en la ecuación obtenemos las dos ecuaciones la ecuación auxiliar asociada de (2) es r2 − λ = 0.
X ′′ (x) − λX (x) = 0 y T ′ (t) − λβT (t) = 0, (1)

donde λ es una constante.

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 5 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 6 15

Analicemos los siguientes casos para solución de (2):

• Segundo caso: Si λ = 0, entonces la raı́z es r = 0 y la solución es

√ X (x) = C1 + C2 x.
• Primer caso: Si λ > 0, entonces las raı́ces son r = ± λ y la solución es
√ √
Utilizando las condiciones de contorno: X ′ (0) = 0 y X ′ (L) = 0, se tiene C2 = 0. Por lo
− λx λx que, obtenemos una solución constante X (x) = C1 . Asimismo, reemplazando λ = 0 en
X (x) = C1 e+ C2 e
 √  √ √ √ T ′ (t) − λβT (t) = 0
Derivando esta solución, tenemos: X ′ (x) = C1 − λ e− λx + C2 λe λx . Luego,

( ′ √ √ obtenemos que
X (0) = − λC1 + λC2 = 0 T ′ (t) = 0 ⇐⇒ T (t) = C.
√ √ √ √
X ′ (L) = C1 − λ e− λL + C2 λe λL = 0 Por lo que, para λ = 0 tenemos la siguiente solución:

de donde C1 = C2 = 0. En consecuencia, se obtiene únicamente la solución trivial, la cual a0

u0 = u(x, t) = X(x)T (t) = C1 C = .
se descarta, dado que nuestro propósito es hallar soluciones no triviales del problema (2). 2

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 7 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 8 15
 √
• Tercer caso: Si λ < 0, entonces las raı́ces de la ecuación caracterı́stica son r = ±i −λ, y Además, resolviendo la segunda ecuación diferencial ordinaria de (1): T ′ (t) = λβT (t),
la solución general es: obtenemos
√ √ nπ 2
T (t) = eλβt+C5 = Be−β ( L ) t
 
X(x) = C1 cos −λx + C2 sen −λx .
o equivalentemente
Su derivada respecto a x es: nπ 2
√ √ √ √  Tn (t) = Be−β ( L ) t, n = 1, 2, 3, . . .
′
X (x) = −C1 −λ sen −λx + C2 −λ cos −λx .
Por lo tanto, combinando las dos soluciones
Aplicando las condiciones de frontera:
√ ) t C cos nπ x ,
2
nπ
 
X ′ (0) = C2 −λ = 0 ⇒ C2 = 0. un (x, t) = Xn (x) Tn (t) = Be−β ( L
1 n = 1, 2, 3, . . .
L
Evaluando ahora en x = L:
√ √  Aplicando el Principio de Superposición, agregamos todas las soluciones particulares un y
X ′ (L) = −C1 −λ sen −λL = 0. absorbemos las constantes B y C1 en un nuevo coeficiente an para cada n:
Si C1 = 0, se obtiene la solución trivial de (2).
̸ 0, se requiere que:
Si C1 = ∞ ∞
a0 X a0 X nπ 2
 nπ 
p  nπ 2 u (x, t) = + un (x, t) = + an e−β ( L ) t cos x
−λL = nπ =⇒ λn = − , n = 1, 2, 3, . . . 2 n=1
2 n=1
L
L
y las funciones propias asociadas son:
donde esta función, satisface la ecuación de calor con las condiciones iniciales y de contorno.
p   nπ 
X(x) = C1 cos −λx =⇒ Xn (x) = C1 cos x , n ≥ 1.
L
Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 9 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 10 15

Observación fı́sica sobre el término a0 :

Determinamos los coeficientes {an } utilizando la condición inicial La solución puede interpretarse como la suma de:
a
∞ • una distribución de temperatura en estado estacionario: 0 , constante e independiente del
a0 X  nπ 
2
u (x, 0) = + an cos x = f (x) (3) tiempo, y
2 n=1
L
∞  nπ  nπ 2
x e−β ( L ) t , que se atenúa progresivamente
X
• una distribución transitoria: an cos
Pero (3) es la serie de Fourier de cosenos de f (x) sobre el intervalo [0, T ], la cual nos permite n=1
L
calcular los coeficientes de Fourier mediante: conforme t → ∞.
a0
Desde el punto de vista fı́sico, el término representa el valor promedio de la distribución de
Z L  nπ  2
2 temperatura inicial:
an = f (x) cos x dx, n = 0, 1, 2, · · ·
1 L
Z
L 0 L a0
= f (x) dx,
2 L 0
lo que concuerda con el hecho de que, en ausencia de flujo de calor hacia o desde el exterior, la
conducción tiende a igualar la temperatura en toda la barra.

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 11 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 12 15
 Ejemplo 2
Encuentre la solución del problema de flujo de calor Para n ≥ 1:

2 π
Z π
 ut = 3uxx , 0 < x < π, t > 0
Z  
2 x sen nx π sen nx
−

∂u ∂u an = x cos (nx) dx = dx
(0, t) = (π, t) = 0, t ≥ 0 π 0 π n 0 0 n
 ∂x ∂x  n 

u (x, 0) = x, 0 ≤ x ≤ π 2h cos nx π i 2 (−1) − 1
= 0+ = .
π n2 0 π n2
Solución. Tenemos que β = 3, L = π y f (x) = x, y que la solución está dada por:
Finalmente, la solución es
∞
a0 X 2
u (x, t) = + an e−3n t cos (nx) π
∞  n 
2 X (−1) − 1 −3n2 t
2 n=1 u (x, t) = + e cos (nx)
2 π n=1 n 2

donde ∞
2
Z L  nπ  π 4X 1 2

an = f (x) cos x dx, n = 0, 1, 2, · · · = − e−3(2k−1) t cos ((2k − 1)x) .

L L 2 π (2k − 1)2
0 k=1

Para n = 0,
π π
x2
Z  
2 2
a0 = xdx = = π.
π 0 π 2 0

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 13 15 Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 14 15

Cierre

• ¿Qué aprendimos el dı́a de hoy?

Departamento Académico de Matemática (UNALM) Unidad 5: Ecuaciones diferenciales parciales Ciclo: 2025 - I 15 15