Ecuaciones Diferenciales Parciales Tema 5.

Ecuación de calor unidimensional

SEMANA 04. CLASE 08. JUEVES 23/10/14

1. Deducción de la ecuación del calor en una varilla unidimensional.
Suponga que una varilla circular delgada de longitud L tiene un área transversal A y coincide
con el eje x en el intervalo 0,L  (ver figura). Suponga lo siguiente:

Figura 1. Varilla circular delgada

a. El flujo de calor dentro de la varilla toma lugar solo en la dirección x.

b. La superficie lateral, o curva, de la varilla está aislada; es decir, el calor no escapa de su
superficie.
c. No se genera calor dentro de la varilla.
d. La varilla es homogénea; es decir, su masa por unidad de volumen ρ es una constante.
e. El calor específico γ y la conductividad térmica K del material de la varilla, son constantes.

A fin de obtener la ecuación diferencial parcial que se cumple para la temperatura u(x,t), son
necesarias dos leyes empíricas de conducción del calor:
Ley empírica 1. La cantidad de calor Q en un elemento de masa m es
Q = γmu , (1)
donde u es la temperatura del elemento.
Ley empírica 2. La rapidez de flujo de calor Qt por la sección transversal indicada en la
figura, es proporcional al área A de la sección transversal y la derivada parcial con respecto a
x de la temperatura:
Qt = −KAux . (2)

Puesto que el calor fluye en la dirección de temperatura decreciente, el signo menos en (2) se
usa para asegurar que Qt sea positiva para ux < 0 (flujo de calor a la derecha) y negativa
para ux > 0 (flujo de calor a la izquierda). Si la parte circular de la varilla mostrada en la
figura entre x y x + ∆x es muy delgada, entonces u(x,t) se toma como la temperatura
aproximada en cada punto del intervalo. Ahora la masa de la porción es m = ρA∆x , y de (1) se
deduce que la cantidad de calor en ésta es
Q = γρA∆xu . (3)

Además, cuando el calor fluye en la dirección positiva de x, se ve de (2) que el calor se

acumula en el corte con razón neta
−KAux (x, t) − −KAux (x + ∆x, t) = KA ux (x + ∆x, t) − ux (x, t) . (4)

Al diferenciar (3) con respecto a t, se ve que esta razón neta se expresa también mediante

José Luis Quintero 1

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Qt = γρ∆xut . (5)
Al igualar (4) y (5), se obtiene
K ux (x + ∆x, t) − ux (x, t)
= ut . (6)
γρ ∆x

Por último, al tomar el límite de (6) cuando ∆x → 0 , se obtiene la ecuación (K / γρ)uxx = ut (La
definición de la segunda derivada parcial viene dada por la expresión
u (x + ∆x, t) − ux (x, t)
uxx = lím x .
∆x → 0 ∆x
Normalmente se establece que k = K / γρ , y esta constante se llama difusividad térmica.

2. Ejemplo ilustrativo 1 (Caso sencillo). Considere una barra delgada o varilla de largo L con
una distribución longitudinal de temperatura f(x) y cuyos extremos se mantienen a una
temperatura constante de cero grados en todo instante. Si la varilla satisface las suposiciones
enunciadas en la sección anterior, entonces la temperatura u(x,t) en la varilla se determina a
partir del problema de valores en la frontera
∂2u ∂u
k 2 = , k > 0, 0 < x < L, t > 0 , (7)
∂x ∂t
u(0, t) = 0, u(L, t) = 0, t > 0 , (8)
u(x, 0) = f(x), 0 < x < L . (9)

Este problema representa el caso de una ecuación de calor homogénea con condición inicial no
nula y condiciones de frontera nulas.
Solución.
Usando el método de separación de variables con la suposición u(x, t) = X(x)T(t) ≠ 0 se tiene
que
X '' T '
=
X kT
Si se asume −λ como constante de separación, se llega a
X '' T '
= = −λ .
X kT
Separando y considerando las condiciones de frontera nulas se tiene que
X ''+ λX = 0 , X(0) = 0 , X(L) = 0 (10)
T '+ KλT = 0 (11)

Ya se ha obtenido resolviendo (10) la solución de la ecuación (7) sujeta a las

condiciones de frontera (8). Se vió que (10) tenía una solución no trivial sólo si el parámetro
λ adoptaba los valores
2
 nπ 
λn =   , n = 1, 2,... (12)
 L 
Las correspondientes soluciones de (10) eran entonces
Xn (x) = Cnsen(nLπx ), n = 1, 2,... (13)

José Luis Quintero 2

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Se hace notar que para un valor de λ distinto de los dados en (12), la única solución de (10)
es la función X(x) = 0 . Puesto que la solución de (11) es
2 π2 /L2 )t
Tn (t) = Dne−kλnt = Dne−k(n ,
los productos
2 π2 /L2 )t  nπ 
un = Xn (x)Tn(t) = Ane−k(n sen  x (14)
 L 
satisfacen la ecuación diferencial (7) y las condiciones de frontera (8) para cada valor del
entero positivo n. Por conveniencia se han reemplazado las constantes CnDn por An . Para que
las funciones dadas en (14) satisfagan la condición inicial (9), se tendría que elegir
coeficientes constantes An de modo que
 nπ 
un (x, 0) = f(x) = Ansen  x . (15)
 L 
Por el principio de superposición, la función
∞ ∞

u(x, t) =
∑ ∑
n =1
un =
n =1
Ane−k(n
2 π2 /L2 )t  nπ 
sen  x
 L 
(16)

también satisface (7) y (8). Sustituyendo t = 0 en (16) se obtiene

∞

u(x, 0) = f(x) =
∑n =1
 nπ 
Ansen  x .
 L 

Se ve que esta última expresión es el desarrollo en serie de senos en medio intervalo. Por
consiguiente si se hace la identificación An = bn , n = 1, 2,3,... , se tiene que
L

∫
2  nπ 
An = f(x)sen  x  dx .
L
0
 L 
Se concluye que la solución del problema de condición de frontera descrito en (7), (8) y (9)
está dada por la serie infinita
∞
 
∑∫
L
2   nπ   −k(n2π2 /L2 )t  nπ 
u(x, t) = f(x)sen  x  dx e sen  x .
L   L    L 
n =1  0 

3. Ejemplo ilustrativo 2 (Caso general). La temperatura u(x,t) en una varilla no aislada de

longitud L, con condiciones de frontera no nulas se determina a partir del problema de valores
en la frontera no homogéneo
kuxx + f(x, t) = ut , k > 0, 0 < x < L, t > 0 ,
u(0, t) = α1(t), u(L, t) = α2 (t), t > 0 ,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) . En tal sentido
k Fxx + wxx  + f(x, t) = Ft + wt , k > 0, 0 < x < L, t > 0 ,
F(0, t) + w(0, t) = α1(t), F(L, t) + w(L, t) = α2 (t), t > 0,
F(x, 0) + w(0, t) = ϕ(x), 0 < x < L .
Construcción de F(x,t).
Si se exige que Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) , el problema anterior se convierte en

José Luis Quintero 3

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

kwxx + f(x, t) = Ft + wt , k > 0, 0 < x < L, t > 0 ,

w(0, t) = 0, w(L, t) = 0, t > 0,
w(0, t) = ϕ(x) − F(x, 0), 0 < x < L .
el cual equivale a
kwxx + G(x, t) = wt , k > 0, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0,
w(0, t) = H(x), 0 < x < L .
donde G(x, t) = f(x, t) − Ft (x, t) y H(x) = ϕ(x) − F(x, 0)
El problema anterior es un problema de calor no homogéneo con condiciones de frontera nulas
y condición inicial no nula.
Volviendo a la construcción de F(x,t) se tiene el siguiente problema:
Problema.
Determinar F(x,t) si Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) .
Solución.

Fxx (x, t) = 0 ⇒ Fx (x, t) =

∫ Fxx (x, t)dx =
∫ 0dx = g2 (t)

⇒ F(x, t) =
∫ ∫
Fx (x, t)dx = g2 (t) dx = g1(t) + xg2 (t)

Al imponer la primera condición de borde se tiene que F(0, t) = α1(t) = g1(t) . Al imponer la
segunda condición de borde se tiene que
α2 (t) −α1(t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
.
Por lo tanto se obtiene la función definida por
F(x, t) = α1(t) + x
L α2 (t) − α1(t) .
En tal sentido, G(x, t) = f(x, t) − α1' (t) − x  α2' (t) − α1' (t) y H(x) = ϕ(x) − α1(0) − x
α2 (0) − α1(0) .
L   L

Construcción de w(x,t).
Ensayos y expresiones de interés:
Sean
∞ ∞ ∞

w(x, t) =
∑
n =1
wn(t)sen(nLπx ) , wt (x, t) =
∑
n =1
wn' (t)sen(nLπx ) , wxx (x, t) = −
∑n =1
wn(t)(nLπ )2 sen(nLπx ) ,

∑
L

∫
2
G(x, t) = Gn (t)sen(nLπx ) , con Gn (t) = G(x, t)sen(nLπx )dx, n = 1,...
L
0
n =1

Resolución de la ecuación diferencial parcial con los ensayos y expresiones de

interés:
Sustituyendo en la ecuación kwxx + G(x, t) = wt se tiene que
∞ ∞ ∞

−k
∑
n =1
wn (t)(nLπ )2 sen(nLπx ) +
∑n =1
Gn (t)sen(nLπx ) =
∑
n =1
wn' (t)sen(nLπx ) .

Reacomodando términos:
∞

∑ n =1
wn' (t) + k(nπ )2 wn (t) − Gn(t) sen(nπx ) = 0 ⇒ wn' (t) + k(nπ )2 wn(t) = Gn (t) .
 L  L L

 
Ecuación lineal de primer orden

José Luis Quintero 4

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Como
∞

w(x, t) =
∑n =1
wn(t)sen(nLπx ) ,

entonces de la condición inicial del problema se tiene

∞

w(x, 0) = H(x) =
∑
n =1
wn (0)sen(nLπx )

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
L

∫
2
wn(0) = H(x)sen(nLπx )dx, n = 1,...
L
0

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de primer orden
sujeta a una condición:
L

∫
2
wn' (t) + k(nLπ )2 wn (t) = Gn(t) , wn(0) = H(x)sen(nLπx )dx, n = 1,...
L
0

Solución general:
 
∫
−k( nπ )2 t k(nπ )2 t −k(nπ )2 t
wn(t) = e L  Gn (t)e L dt + C = e L Mn(t) + C  .
 
Valor de C: C = wn (0) − Mn (0) .
Solución particular:
−k( nπ )2 t
wn(t) = e L Mn (t) + wn (0) − Mn(0) .

Resumen de pasos a seguir:

Paso 1. Construcción de G(x,t).

G(x, t) = f(x, t) − Ft (x, t) = f(x, t) − α1' (t) − x α2' (t) − α1' (t)
L  
Paso 2. Cálculo de los Gn (t) .
L

∫
2
Gn (t) = G(x, t)sen(nLπx )dx, n = 1,...
L
0

Paso 3. Construcción de H(x).

H(x) = ϕ(x) − F(x, 0) = ϕ(x) − α1(0) − x
L
α2 (0) − α1(0)
Paso 4. Cálculo de wn (0) .
L

∫
2
wn(0) = H(x)sen(nLπx )dx, n = 1,...
L
0
Paso 5. Cálculo de Mn (t) y Mn (0) .

∫
k(nπ )2 t
Mn (t) = Gn (t)e L dt

Paso 6. Cálculo de los wn(t) .

−k( nπ )2 t
wn(t) = e L Mn (t) + wn (0) − Mn(0)

José Luis Quintero 5

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Paso 7. Construcción de w(x,t).

∞

w(x, t) =
∑
n =1
wn(t)sen(nLπx )

Paso 8. Construcción de u(x,t).

∞

u(x, t) = F(x, t) + w(x, t) = α1(t) + x

L
α2 (t) − α1(t) +
∑
n =1
wn (t)sen(nLπx ) .

4. Ejemplo ilustrativo 3. Resolver

uxx + 1 = ut , 0 < x < π, t > 0
u(0, t) = 2, u(π, t) = 3, t > 0
u(x, 0) = 4, 0 < x < π
Solución.
Paso 1. Construcción de G(x,t).
G(x, t) = f(x, t) − Ft (x, t) = 1 − 0 = 1
Paso 2. Cálculo de los Gn (t) .
π

∫
2 2 2 
Gn (t) = sen(nx)dx = 1 − cos(nπ) = 1 − (−1)n 
π πn  πn  
0

Paso 3. Construcción de H(x).

x x
H(x) = ϕ(x) − F(x, 0) = 4 − 2 − =2−
π π
Paso 4. Cálculo de wn (0) .
π π π

∫ ∫ ∫
2  x 4 2
wn(0) =  2 − π  sen(nx)dx = π sen(nx)dx − xsen(nx)dx
π   π2
0 0 0
4  2 4 2 2 
= 1 − (−1)n  + (−1)n = − (−1)n = 2 − (−1)n 
πn   πn πn πn πn  
Paso 5. Cálculo de Mn (t) y Mn (0) .

∫
2  2 2  2 2 
Mn(t) = 1 − (−1)n  en t dt = 1 − (−1)n  en t ⇒ Mn (0) = 1 − (−1)n 
πn   πn3   πn3  
Paso 6. Cálculo de los wn(t) .
2t  2  n  n2t 2  2 
wn(t) = e−n  3 1 − (−1)  e + 2 − (−1)n  − 3 1 − (−1)n  
 πn πn   πn   
Paso 7. Construcción de w(x,t).
∞

w(x, t) =
∑
n =1
e−n
2t  2 
 πn
n  n2t 2  n 2  n 
 3 1 − (−1)  e + πn 2 − (−1)  − 3 1 − (−1)   sen(nx)
πn 

Paso 8. Construcción de u(x,t).

∞
x
u(x, t) = 2 + +
π ∑ n =1
e−n
2t  2 
 πn
n  n2t 2  n 2  n 
 3 1 − (−1)  e + πn 2 − (−1)  − 3 1 − (−1)   sen(nx)
πn 

José Luis Quintero 6

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

5. Ejemplo ilustrativo 4 (Caso general). La temperatura u(x,t) en una varilla no aislada de

longitud L, con condiciones de frontera no nulas donde hay flujo de calor en la superficie de la
varilla en estudio hacia el medio que la rodea cuya temperatura se mantiene constante,
digamos um se determina a partir del problema de valores en la frontera no homogéneo
kuxx + β u + δ = ut , k > 0, β, δ ∈ R, 0 < x < L, t > 0
u(0, t) = α1(t), u(L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) . En tal sentido
k Fxx + wxx  + β F(x, t) + w(x, t) + δ = Ft + wt , k > 0, β, δ ∈ R, 0 < x < L, t > 0 ,
F(0, t) + w(0, t) = α1(t), F(L, t) + w(L, t) = α2 (t), t > 0,
F(x, 0) + w(0, t) = ϕ(x), 0 < x < L .
Construcción de F(x,t).
Si se exige que Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) , el problema anterior se convierte en
kwxx + β F(x, t) + w(x, t) + δ = Ft + wt , k > 0, β, δ ∈ R, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0,
w(0, t) = ϕ(x) − F(x, 0), 0 < x < L .
el cual equivale a
kwxx + βw(x, t) + G(x, t) = wt , k > 0, β ∈ R, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0,
w(0, t) = H(x), 0 < x < L .
donde G(x, t) = βF(x, t) + δ − Ft (x, t) y H(x) = ϕ(x) − F(x, 0)
El problema anterior es un problema de calor no homogéneo con condiciones de frontera nulas
y condición inicial no nula.
Volviendo a la construcción de F(x,t) se tiene el siguiente problema:
Problema.
Determinar F(x,t) si Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) .
Solución.

Fxx (x, t) = 0 ⇒ Fx (x, t) =

∫ Fxx (x, t)dx =
∫ 0dx = g2 (t)

⇒ F(x, t) =
∫ ∫
Fx (x, t)dx = g2 (t) dx = g1(t) + xg2 (t)

Al imponer la primera condición de borde se tiene que F(0, t) = α1(t) = g1(t) . Al imponer la
segunda condición de borde se tiene que
α2 (t) −α1(t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
.
Por lo tanto se obtiene la función definida por
F(x, t) = α1(t) + x
L
α2 (t) − α1(t) .
En tal sentido, G(x, t) = β α1(t) + x
L
(α2 (t) − α1(t)) + δ − α1' (t) − Lx (α'2 (t) − α1' (t)) y
H(x) = ϕ(x) − α1(0) − x
L α2 (0) − α1(0) .
Construcción de w(x,t).
Ensayos y expresiones de interés:
Sean

José Luis Quintero 7

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

∞ ∞ ∞

w(x, t) =
∑
n =1
wn(t)sen(nLπx ) , wt (x, t) =
∑n =1
wn' (t)sen(nLπx ) , wxx (x, t) = −
∑ n =1
wn(t)(nLπ )2 sen(nLπx ) ,

∑
L

∫
2
G(x, t) = Gn (t)sen(nLπx ) , con Gn (t) = G(x, t)sen(nLπx )dx, n = 1,...
L
0
n =1

Resolución de la ecuación diferencial parcial con los ensayos y expresiones de

interés:
Sustituyendo en la ecuación kwxx + β w(x, t) + G(x, t) = wt se tiene que
∞ ∞ ∞ ∞

−k
∑
n =1
wn (t)(nLπ )2 sen(nLπx ) +β
∑ n =1
wn (t)sen(nLπx ) +
∑
n =1
Gn (t)sen(nLπx ) =
∑n =1
wn' (t)sen(nLπx ) .

Reacomodando términos:
∞

∑n =1
w' (t) + k(nπ )2 − β  w (t) − G (t) sen(nπx ) = 0 ⇒ w' (t) + k(nπ )2 − β  w (t) = G (t) .
 n  L  n n  L n  L  n


n

Ecuación lineal de primer orden

Como
∞

w(x, t) =
∑n =1
wn(t)sen(nLπx ) ,

entonces de la condición inicial del problema se tiene

∞

w(x, 0) = H(x) =
∑
n =1
wn (0)sen(nLπx )

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
L

∫
2
wn(0) = H(x)sen(nLπx )dx, n = 1,...
L
0

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de primer orden
sujeta a una condición:
L

∫
2
wn' (t) + k(nLπ )2 − β  wn (t) = Gn(t) , wn(0) = H(x)sen(nLπx )dx, n = 1,...
  L
0

Solución general:
β− k( nπ )2  t  k( nπ )2 −β  t  β − k(nπ )2  t
wn(t) = e  L 

 ∫ Gn (t)e  L 
dt + C  = e 

L 
Mn (t) + C  .

Valor de C: C = wn (0) − Mn (0) .

Solución particular:
β −k( nπ )2  t
wn(t) = e  L 
Mn (t) + wn (0) − Mn (0) .
Resumen de pasos a seguir:

Paso 1. Construcción de G(x,t).

G(x, t) = β α1(t) + x
L
(α2 (t) − α1(t)) + δ − α1' (t) − Lx (α'2 (t) − α1' (t))

José Luis Quintero 8

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Paso 2. Cálculo de los Gn (t) .

L

∫
2
Gn (t) = G(x, t)sen(nLπx )dx, n = 1,...
L
0

Paso 3. Construcción de H(x).

H(x) = ϕ(x) − F(x, 0) = ϕ(x) − α1(0) − x
L
α2 (0) − α1(0)
Paso 4. Cálculo de wn (0) .
L

∫
2
wn(0) = H(x)sen(nLπx )dx, n = 1,...
L
0
Paso 5. Cálculo de Mn (t) y Mn (0) .
k( nπ )2 −β  t
Mn(t) =
∫ Gn (t)e 
 L  dt

Paso 6. Cálculo de los wn(t) .

β− k(nπ )2  t
wn(t) = e  L 
Mn (t) + wn (0) − Mn (0)
Paso 7. Construcción de w(x,t).
∞

w(x, t) =
∑ n =1
wn(t)sen(nLπx )

Paso 8. Construcción de u(x,t).

∞

u(x, t) = F(x, t) + w(x, t) = α1(t) + x

L
α2 (t) − α1(t) +
∑n =1
wn (t)sen(nLπx ) .

6. Ejemplo ilustrativo 5 (Caso particular del ejemplo ilustrativo 4). Si hay flujo de calor
en la superficie de la varilla en estudio hacia el medio que la rodea cuya temperatura se
mantiene constante, digamos um la ecuación de calor se convierte en kuxx − h(u − um ) = ut . La
temperatura u(x,t) de la varilla es u(0, t) = u(L, t) = 0 y u(x, 0) = f(x) . Demuestre que la
solución de la ecuación kuxx − hu = ut sujeta a las condiciones anteriores viene dada por la
expresión
∞

u(x, t) = e−ht
∑ n =1
cne−kn
2 π2 t /L2  nπx 
sen 
 L 
,

donde
L

∫
2  nπx 
cn = f(x)sen   dx .
L  L 
0

Solución.
Paso 1. Construcción de G(x,t). G(x, t) = 0 .
Paso 2. Cálculo de los Gn (t) . Gn (t) = 0, n = 1,...
Paso 3. Construcción de H(x). H(x) = f(x) .
Paso 4. Cálculo de wn (0) .
L

∫
2
wn(0) = f(x)sen(nLπx )dx, n = 1,...
L
0

José Luis Quintero 9

Ecuaciones Diferenciales Parciales Tema 5. Ecuación de calor unidimensional

Paso 5. Cálculo de Mn (t) y Mn (0) . Mn(t) = 0 ⇒ Mn(0) = 0 .

Paso 6. Cálculo de los wn(t) .
L
−  −h − k(nπ )2  t − h + k( nπ )2  t

∫
    2
wn(t) = wn (0)e L
=e L
f(x)sen(nLπx )dx
L
0

Paso 8. Construcción de u(x,t).

∞

∑
L

∫
−ht −k( nπ )2 t 2
u(x, t) = F(x, t) + w(x, t) = e cne L sen(nLπx ) , cn = f(x)sen(nLπx )dx, n = 1,...
L
0
n =1

José Luis Quintero 10