República Bolivariana de Venezuela

Universidad Nacional Abierta

Área de Matemática

ECUACIÓN GENERAL DE CALOR EN R3

Y UNA APLICACIÓN EN MAPLE

Tutor: Javier Quintero

Gualberto Adolfo Salinas Linares

V-16444149
0414-7582976
gualberto164@[Link]
2

INDICE

Pag.

I. INTRODUCCIÓN ...................................................................................................................... 3

II. PLANTEAMIENTO DEL PROBLEMA ....................................................................................... 4

III. OBEJTIVOS DEL TRABAJO ....................................................................................................... 5

I. OBJETIVO GENERAL ........................................................................................................ 5

II. OBJETIVOS ESPECÍFICOS ................................................................................................ 5

IV. RELACIÓN DE LAS ACTIVIDADES A SEGUIR Y UN CRONOGRAMA TENTATIVO

DEL DESARROLLO DE LAS MISMAS ...................................................................................... 6

V. ENFOQUE METODOLÓGICO .................................................................................................. 7

VI. BIBLIOGRAFÍA O REFERENCIAS ............................................................................................. 8

Capítulo 1

Series de Fourier

1.1. Funciones periódicas, funciones pares y funciones impares

Definición 1 Una función f : A ⊆ R → R es peródica si existe un número T > 0 tal que

f (t + T ) = f (t ) ∀t ∈ A (1,1)

La constante mínima T que satisface la relación (1.1) s llama peíodo de la función.

Definición 2 Una función f : A ⊆ R → R es par sí

f (−t ) = f (t ) ∀ ∈ A.

Una función f : A ⊆ R → R es impar sí

f (−t ) = − f (t ) ∀ ∈ A.

Propiedades de las funciones pares e impares.

a) El gráfico de una función par es simétrico respecto al eje vertical.

b) El gráfico de una función impar es simétrico respecto al origen.

c) Si f [−a, a] → R es par, entonces

Z a Z a
f (t )d t = 2 f (t )d t
−a 0

d) Si f [−a, a] → R es impar, entonces

Z a
f (t )d t = 0
−a

e) Sean f 1 , f 2 : A ⊆ R → R funciones impares. Entonces

1) Suma: f 1 + f 2 es impar.
2) Producto: f 1 · f 2 es par.
3) Derivada: f 1′ es par.

3
4 CAPÍTULO 1. SERIES DE FOURIER

f) Sean f 1 , f 2 : A ⊆ R → R funciones pares. Entonces

1) Suma: f 1 + f 2 es par.
2) Producto: f 1 · f 2 es par.
3) Derivada: f 1′ es impar.
1.2. SERIES DE FOURIER 5

1.2. Series de Fourier

El conjunto de funciones

πt πt
½ µ ¶ µ ¶ µ ¶ µ ¶ ¾
2πt 2πt
1, cos , cos , . . . , sen , sen ,...
p p p p
es ortogonal en el intervalo [−p, p]. Ahora dada una función f : [−p, p] → R, lacual es represen-
table mediante una cierta serie trigonométrica, es decir, para cada t ∈ [−p, p] la serie trigonomé-
trica correspondiente converge a la función f (t ), satisfaciendo la igualdad

a0 X∞ · µ
nπ
¶ µ
nπ
¶¸
f (t ) = + a n cos t + b n sen t (I 1 )
2 n=1 p p
Nos preguntamos: ¿Es posible determinar un conjunto de coeficientes a 0 , a n , b n para n = 1, 2, 3, . . .,
de la serie (I 1 )?

Para ello, integramos a ambos lados de (I 1 ) sobre [−p, p] resulta

Z p a0
Z p ∞
X
· Z p nπ
¶
µ Z p µ
nπ
¶ ¸
f (t )d t = dt + an cos t d t + bn sen t dt
−p 2 −p n=1 −p p −p p

Como cada una de las funciones cos( nπ nπ

p t ), sen( p t ) para n ≥ 1 es ortogonal a 1 en [−p, p], se
tiene:
Z p a0 1
Z p
f (t )d t = 2p =⇒ a0 = f (t )d t (I 2 )
−p 2 p −p
³ ´
nπ
Por otro lado, multiplicando a ambos lados de (I 1 ) por cos p t e integrando a ambos lados
sobre [−p, p], obtenemos

Z p nπ
¶µ
a0 p
Z µ
nπ
¶
X∞ · Z p µ
nπ
¶ µ
nπ
¶
cos t f (t )d t = cos t dt + an cos t cos t dt+
−p p 2 −p p n=1 −p p p
Z p
nπ nπ
µ ¶ µ ¶
bn cos t sen t dt (I 3 )
−p p p

Como
Z p nπ
µ¶
cos t d t = 0, m>0
−p p
(
Z p nπ
µ¶ µ
nπ
¶
0, si m , n
cos t cos t dt =
−p p p p, si m = n
Z p
nπ nπ
µ ¶ µ ¶
cos t sen t dt = 0
−p p p
En consecuencia, (I 2 ) se reduce a
Z p nπ
µ ¶
nπ
¶
1
Z p µ
f (t )cos t d t = an p
=⇒ f (t )cos t dt an =
(I 4 )
−p p −p p p
³ ´
nπ
Finalmente, multiplicando a ambos lados de (I 1 ) por sen p t e integrando,
6 CAPÍTULO 1. SERIES DE FOURIER

Z p nπ
µ ¶
a0 p
Z µ
nπ
¶
X∞ · Z p µ
nπ
¶ µ
nπ
¶
sen t f (t )d t = sen t dt + an sen t cos t dt+
−p p 2 −p p n=1 −p p p
Z p
nπ nπ
µ ¶ µ ¶
bn sen t sen t dt (I 5 )
−p p p

Como
Z p nπ
µ ¶
sen t dt = 0 ,m > 0
−p p
Z p
nπ nπ
µ ¶ µ ¶
sen t cos t dt = 0
−p p p
Z p (
nπ nπ 0, si m , n
µ ¶ µ ¶
sen t sen t dt =
−p p p p, si m = n
De modo que (I 5 ) se reduce a
Z p p
nπ nπ
µ ¶ µ ¶
1
Z
f (t )cos t d t = bn p =⇒ bn = f (t )sen t dt (I 6 )
−p p p −p p
A la serie trigonométrica (I 1 ), con coeficientes a 0 , a n y b n definidas por (I 2 ), (I 4 ) y (I 6 ), respec-
tivamente, se le llama Serie de Fourier de la función f . Los coeficientes obtenidos de (I 2 ), (I 4 ) y (I 6 )
se conocen como coeficientes de f .

1.3. Ecuaciones diferenciales

Cuando hablamos de una ecuación diferencial, esto nos lleva a pensar en una ecuación con
derivadas continuas y ′ , y ′′ , . . .. Por ello cuando resolvamos ecuaciones diferenciales nos topare-
mos con ecuaciones del estilo y ′′ + 3y ′ + y = 0. Dada la importancia del estudio de las ecuaciones
diferenciales, las cuales se presentan en numerosos modelos naturales, pasaremos a definirla.

Definición 3 Una ecuación que contenga una o más derivadas en las variables dependientes, con
respecto a una o más variables independientes se llamará ecuación diferencial (ED).

Las ecuaciones diferenciales las podemos clasificar según su tipo, orden y linealidad. Si clasifi-
camos las EDs por su orden, estaremos interesados en el mayor orden de derivada que aparezca,
por ejemplo:

d2y dy
+ 3 − 2y = e x
d x2 dx
En este caso diremos que la anterior ecuación diferencial es de orden dos, ya que la derivada
de mayor orden es dos. Si por el contrario, nos interesa clasificar las EDs por su linealidad, en-
tonces tendremos que son lineales o no lineales. Así, una ED de n-orden se dice que es lineal si
F (x, y ′ , y ′′ , . . . , y (n) ) = 0 es lineal en y ′ , y ′′ , . . . , y (n) . Es decir, F es lineal cuando:

a n (x)y (n) + a n−1 (x)y (n−1) + · · · + a 1 (x)y ′ + a 0 (x)y − g (x) = 0

Donde los coeficientes a 0 , a 1 , . . . , a n de y, y ′ , . . . , y (n) dependen de la variable independiente x.
1.4. ECUACIÓN DIFERENCIAL DE CALOR EN UN MEDIO ISOTRÓPICO 7

Ahora, si queremos clasificar las ecuaciones diferenciales por su tipo, tendremos que clasifi-
carlas como ecuaciones diferenciales ordinarias (EDO) o ecuaciones diferenciales parciales (EDP).

Entonces, si una ecuación contiene solamente derivadas ordinarias de una o más variables
dependientes con respecto a una sola variable independiente, se le llama ecuación diferencial
ordinaria (EDO). Como ejemplo tenemos:

dy d2y d y dy dy
− 9y = e x , + − 4y = 0 y + = 3x − y
dx d x2 d x dt dt
Ahora, una ecuación que posea derivadas parciales de una o más variables dependientes de
dos o más variables independientes la llamaremos ecuación diferencial parcial (EDP). Algunas
EDPs son:

∂2 u ∂2 u ∂2 u ∂2 u ∂u
+ =0 , = 2 +3
∂x 2 ∂y 2 ∂x 2 ∂t ∂t
A su vez, las EDPs tambien tienen su clasificación. Sea la siguiente ecuación diferencial:

∂2 u ∂u ∂2 u du du
A 2 + 2B +C 2 + D +E +Fu +E = 0 (a)
∂x ∂x∂y ∂y dx dy
Donde los coeficientes A, B, . . . , E son constantes. Entonces se tiene la siguiente definición

Definición 4 La ecuación diferencial parcial de segundo orden lineal con coeficientes (a) se llama:
de tipo hiperbólico si B 2 − AC > 0,
de tipo parabólico si B 2 − AC = 0 y
de tipo elíptico si B 2 − AC < 0

En nuestro estudio nos centraremos en las ecuaciones diferenciales lineales de segundo orden
parabólicas, ya que la ecuación de calor es de tipo parabólico.

1.4. Ecuación diferencial de calor en un medio isotrópico

La ecuación diferencial fundamental de calor en un medio isotrópico se deduce de la ecuación
(1.1), explicaremos y demostraremos la ecuación de conductividad de calor: Para ello, considere-
mos el caso de un sólido a través del cual fluye calor, por ejemplo, un elemento de volumen en
forma de paralelepípedo rectangular con centro en un punto P (x, y, z), en este trazamos ejes coor-
denados x, y y z que pasan por el punto P de longitudes 2d x, 2d y y 2d z.

Denotemos por T = T (t , x, y, z) la concentración de calor del cuerpo en el punto P (x, y, z) en el

tiempo t .

Como la transferencia de calor en un sólido predomina la conductividad de estos, entonces se

tiene que la fórmula de Fourier es una cantidad unidireccional, es decir, mide le flujo de calor en
una dirección, así tenemos los flujos de calor en tres direcciones dada por:

∂T
F x = −k4d yd z (1,2)
∂x

∂T
F y = −k4d xd z (1,3)
∂y
8 CAPÍTULO 1. SERIES DE FOURIER

∂T
F z = −k4d xd y (1,4)
∂z
Donde

k es la conductividad térmica del sólido, para el flujo unidireccional (según x,y,z).

F x , F y y F z es el fujo de calor que pasa según el área de sección transversal a travé del puntp
P en la dirección del eje x, eje y y eje z.

4d yd z, 4d xd z y 4d xd y es el área de sección transversal de la transferencia de calor.

∂T ∂T
y ∂T
∂x , ∂y ∂z nos dan las variaciones de la temperatura T en las direcciones de los ejes coor-
denados en las proximidades del punto P .

El flujo de calor va de alta a baja temperatura, para que exista un flujo de calor debe haber
una variación de temperatura. Ahora como estamos analizando este paralelepípedo rectangular,
vamos a analizar todos los flujos que pasan a través de él, así tomando ABC D y A ′ B ′C ′ D ′ las caras
perpendiculares al eje x como en la figura 1.2. Entonces el flujo de calor que ingresa al sólido a
travé de la cara ABC D en dirección de x − d x viene dada por

∂F x
F x−d x = F x − d x
∂x
De manera análoga, el flujo de calor que sale del sólido a través de la cara A ′ B ′C ′ D ′ viene dada
por

∂F x
F x+d x = F x + d x
∂x
Así,

∂F y
F y−d y = F y − d y
∂y
Análogamente, para las otras caras

∂F y
F y+d y = F y + d y
∂y
∂F z
F z+d z = F z + d z
∂z
Como estamos analizando los flujos de calor (o flujos de energía) debe existir un balance de
energía en este sólido, esto es lo que se conoce como la conservación de la energía (primera ley de
la termodinámica):

Ė ent r ad a + Ė g ener ad a − Ė sol i d o = Ė acumul ad a

Como hemos considerado que dentro del sólido no se genera calor, es decir, Ė g ener ad a = 0. Así,
se tiene que

Ė ent r ad a − Ė sol i d o = Ė acumul ad a

esto es,

∂F x ∂F y ∂F z ∂F x ∂F y ∂F z
µ ¶ µ ¶
Fx − d x + Fy − d y + Fz − d z − Fx + d x + Fy + d y + Fz + d z = Ė acum (1,5)
∂x ∂y ∂z ∂x ∂y ∂z
1.4. ECUACIÓN DIFERENCIAL DE CALOR EN UN MEDIO ISOTRÓPICO 9

Ahora, la energía almacenada

∂T
Ė al m = d mc (1,6)
∂t
donde d m = ρdV , m = masa del bloque, como el material es muy pequeño la masa viene dada
por un d m, ρ es la densidad y c el calor específico, de modo que ρc es la capacidad calorífica por
unidad de volumen.

De (1.5) y (1.6), resulta

∂F x ∂F y ∂F z ∂T
−2d x − 2d y − 2d z = ρcdV (1,7)
∂x ∂y ∂z ∂t
Sustituyendo (1.2), (1.3) y (1.4) en (1.7) se tiene

∂2 T ∂2 T ∂2 T ∂T
µ ¶ µ ¶ µ ¶
−2d x −k4d yd z 2 − 2d y −k4d xd z 2 − 2d z −k4d xd y 2 = ρcdV
∂x ∂y ∂z ∂t

∂2 T ∂2 T ∂2 T ∂T
k8d xd yd z + k8d xd yd z + k8d xd yd z = ρcdV (1,8)
∂x 2 ∂y 2 ∂z 2 ∂t
Como la diferencia de volumen es dV = 8d xd yd z, sustituyendo en (1.8)

∂2 T ∂2 T ∂2 T ∂T
µ ¶
kdV + 2 + 2 = ρcdV
∂x 2 ∂y ∂z ∂t
resulta

∂T k ∂2 T ∂2 T ∂2 T
µ ¶
= + + (1,9)
∂t ρc ∂x 2 ∂y 2 ∂z 2
Esta es la ecuación de calor en coordenadas cartesianas en R3 , la ecuación (1.9) se puede res-
cribir en forma matricial

∂T ∂ k ∂T ∂ k ∂T ∂ k ∂T
µ ¶ µ ¶ µ ¶
= + +
∂t ∂x ρc ∂x ∂y ρc ∂y ∂z ρc ∂z
k ∂T k ∂T k ∂T
µ ¶
=▽ + +
ρc ∂x ρc ∂y ρc ∂z
 k  
ρc
0 0
k
= ▽  0 ρc 0 ▽T
  
k
0 0 ρc
k
· ¸
= di v ▽T
ρc
Ahora bien, estamos interesados en la ecuación unidimensional en estado no estacionario con
la variación de la temperatura T según la posición de x y el tiempo t en el cubo con centro P (x, y, z)
de longitud 2a en el eje x y de temperatura inicial f (x) que se extiende a lo largo del eje x y cuyos
extremos se mantienen a una temperatura constante T1 y temperatura T0 . Si

El flujo de calor se produce solamente en la dirección del eje x.

No se pierde calor a través de las superficies laterales del cubo.

No se genera calor dentro del cubo.

10 CAPÍTULO 1. SERIES DE FOURIER

El cubo es homogéneo, es decir, su densidad por unidad de longitud es constante.

Su calor específico es constante, entonces la temperatura T (t , x) del cubo en la dirección del

eje x esta dada por la solución del problema con condiciones iniciales y de frontera.

La misma cumple una condición inicial la cual se puede específicar en la forma:

T (x, 0) = f (x) x ∈ [0, a] (1,11)

y las condiciones de frontera no homogenea de temperatura específica se puede dar como:

T (0, t ) = T1 , T (a, t ) t ≥0 (1,12)

donde T1 , T2 son temperaturas específicas constantes en la superficie en x = 0 y x = a respec-
tivamente.

Ahora, cuando la concentración de calor es constante, soluciones generales de la ecuación de

calor se puede obtener para las condiciones inicieles y de frontera, para ello existe el método de
separación de variables o método de Fourier como tambien se conoce, siqgue siendo un método
de gran importancia y de uso frecuente en la actualidad, para este método a la ecuación (1.10),
(1.11) y (1.12) consideramos que una solución de esta ecuación es dada por:

T (x, t ) = Y (x) +C (x, t ) (1,13)

donde Y (x) es la concentración del estado estacionario, sin generación de calor y C (x, t ) es la
concentración del estado no estacionario, las cuales verifican:

∂2 Y (x)
(
∂x 2
=0 (1,14)
Y (0) = T1 , .Y (a) = T2 (1,15)
y

∂C (x,t ) ∂2C (x,t )
 ∂t = k 2 ∂x 2

 (x, t ) ∈ [0, a] × [0, +∞) (1,16)
C (x, 0) = f (x) − Y (x) 0<x <a (1,17)


C (0, t ) = 0,C (a, t ) = 0 t >0 (1,18)
Para hallar la soluciòn de (1.14) y (1.15) integramos a ambos lados de (1.14) con respecto a x,
resulta

∂2 Y ∂Y
Z
d x = 0 =⇒ = k1 (1,19)
∂x 2 ∂x
integrando a ambos lados de (1.19) con respecto a x, se tiene

∂Y
Z Z
d x = k1 d x =⇒ Y (x) = k 2 + k 1 x (1,20)
∂x
Aplicando las condiciones iniciales de (1.15), obtenemos

Y (0) = T1 =⇒ k 2 + K 1 (0) = T1

Y (a) = T2 =⇒ k 2 + K 1 (a) = T2

(
k 2 = T1
=⇒
k 1 = T2 −T
a
1
1.4. ECUACIÓN DIFERENCIAL DE CALOR EN UN MEDIO ISOTRÓPICO 11

Sustituyendo k 1 y k 2 en (1.20), obtenemos

x
Y (x) = T1 + (T2 − T1 ) (1,21)
a
es la solución de estado estacionario (1.14) - (1.15).

Ahora resolveremos el estado no estacionario, la solución de (1.16) con las condiciones (1.17) -
(1.18) se obtiene aplicando el método de separación de variable. Supongamos que la C (x, t ) tiene
la forma:

C (x, t ) = X (x)U (t ) (1,22)

que es el producto de dos funciones, una de las cuales depende sólo de la variable x y la otra
de la variable t . Mediante diferenciación de (1.22), se obtiene

∂2C ∂C
= X ′′ (x)U (t ) y = X (x)U ′ (t ) (1,23)
∂x 2 ∂t
Sustituyendo (1.23) en (1.16),

X ′′ (x) U ′ (t )
k2 =
X (x) U (t )
La expresión del primer miembro depende solamente de x, mientras que la segunda expresión
sólo depende de t ; por consiguiente, se concluye que ambas expresiones deben ser igual a una
constante. En la práctica es conveniente escribir esta constante real como λ2 , o bien como −λ2 se
tiene

X ′ (x) 1 U ′ (t )
= 2 = −λ2 .
X (x) k U (t )
Y de lo anterior obtenemos dos ecuaciones diferenciales ordinarias

X ′′ (x) + λ2 X (x) = 0 (1,24)

U ′ (t ) + k 2 λ2U (t ) = 0 (1,25)
las soluciones generales de (1.24) y (1.25) son

X (x) = Acos(λx) + B sen(λx) (1,26)

2 2
λ t
U (t ) = C e −k (1,27)
Aplicando las condiciones de frontera (1.18) y del hecho que sen(0) = 0, resulta

X (0) = Acos(0) + B sen(0) =⇒ A=0

X (a) = Acos(λa) + B sen(λa) =⇒ B sen(λa) = 0

Luego, B , 0 y sen(λa) = 0. Si sen(λa) = 0, entonces λa = ar c sen(0) donde λa = nπ∀n =
1, 2, 3, . . .. Por lo tanto,
nπ
λ= ∀n ≥ 1. (1,29)
a
Se obtiene así las soluciones de (1.24), que satisfacen las (1.18) como
12 CAPÍTULO 1. SERIES DE FOURIER

³ nπ ´
X n (x) = B n sen x ∀n ≥ 1 (1,30)
a
Las soluciones correspondientes de (1.27) son

2 π2
−k 2 n t
Un (t ) = C n e a2 ∀n ≥ 1 (1,31)

En consecuencia,

2 π2 ³ nπ ´
−k 2 n t
C n (x, t ) = X n (x)Un (t ) = b n e a2 sen x ∀n ≥ 1
a
donde b n = B n C n , son las soluciones de la ecuación diferencial parcial (1.16) que satisfacen
la condición (1.18), entonces por el principio de superposición, la solución de (1.18) resulta de la
siguiente forma:

∞ 2 π2 ³ nπ ´
X −k 2 n t
C (x, t ) = bn e a2 sen x (1,32)
n=1 a

la cual satisface las condiciones de frontera (1.18). Y por la condición (1.17) y (1.32), se obtiene

∞
X ³ nπ ´
C (x, 0) = b n sen x
n=1 a

³ nπ ´ ∞
X
f (x) − Y (x) = x b n sen
(1,33)
n=1 a

Multiplicando a ambos lados de (1.33) por sen nπ

¡ ¢
a x e integrando respecto a x, se tiene
Z a ³ nπ ´ ∞ Z a ³ nπ ´ ³ nπ ´
′ ′ ′ ′
X
[ f (x ) − Y (x )]sen x dx = bn sen x sen x dx (1,34)
0 a n=1 0 a a

Por otro lado, como

Z a ³ nπ ´ ³ nπ ´
sen x sen x dx = 0 , si m,n
0 a a
Z a ³ nπ ´ a
sen 2 x dx = , si m=n (1,35)
0 a 2
resulta que
Para la solución de este problema nos regiremos por el método de separación de variables. Este
método (inventado en s. XVIII por Daniel Bernoulli) no se puede utilizar en todos los casos, sólo
se limita a ciertas condiciones, las cuales veremos.

Al aplicar el método de separación de variables tendremos en cuante que trabajaremos con

EDPs lineales y homogéneas, conceptos que veremos a continuación.

Definición 5 Un operador lineal se define como aquél que satisface

L(c 1 u 1 + c 2 u 2 ) = c 1 L(u 1 ) + c 2 L(u 2 )

para dos funciones cualesquiera u 1 y u 2 , donde c 1 y c 2 son constantes arbitrarias.

1.4. ECUACIÓN DIFERENCIAL DE CALOR EN UN MEDIO ISOTRÓPICO 13

Una ecuación lineal para la incógnita u es de la forma:

L(u) = f (1,10)
donde L es un operador lineal y f es conocida. Ahora si f = 0, entonces (1.10) se convierte
en L(u) = 0, llamada así una ecuación lineal homogénea. Como ejemplo tenemos la ecuación de
calor unidemensional:

∂u ∂2 u
−k 2 = 0
∂t ∂x
La propiedad fundamental de los operadores lineales nos permite sumar soluciones de ecua-
ciones tal que:

Definición 6 Principio de superposicion. Si u 1 y u 2 satisfacen una ecuación lineal homogénea,

entonces cualquier combinación lineal de ellas, c 1 u 1 + c 2 u 2 , también satisface la misma ecuación
lineal homogénea.

En el método de separación de variables se intenta encontrar soluciones que sean de la forma:

u(x, t ) = φ(x)G(t ) (1,11)

Donde φ(x) es una función sólo de x y G(t ) una función sólo de t . La función definida en (1.11)
debe cumplir la ecuación en derivadas parciales lieal homogénea y las condiciones de contorno.

Así, podemos decir que nuestra meta es solucionar la siguiente ecuación lineal homogénea

∂u ∂2 u
=k 2 0 < x < L, t >0
∂t ∂x
con las siguientes condiciones de contorno

u(0, t ) = 0,
u(L, t ) = 0
y la condición inicial

u(x, 0) = f (x)
El primer paso a seguir será sustiruir la solución producto (1.11) en la ecuación tal que

∂u dG
= φ(x) ,
∂t dt
∂2 u d 2 u
= G(t )
∂x 2 d x 2
quedando así

dG ∂2 φ
φ(x)
= k 2 G(t )
dt ∂x
Podemos notar que podemos separar variables si dividimos entre φ(x)G(t ) y la constante k
quedando así:

1 dG 1 ∂2 φ
= (1,12)
kG d t φ ∂x 2
14 CAPÍTULO 1. SERIES DE FOURIER

Ahora vemos que tenemos una ecuación donde una función del tiempo es igual a otra función
del espacio. Si x y t son variables arbitrarias entonces x no puede ser una función de t , por lo tanto
es necesario que ambos miembros de (1.12) sean iguales a la misma constante:

1 dG 1 ∂2 φ
= = −λ (1,13)
kG d t φ ∂x 2
Donde λ es una constante arbitraria, llamada constante de separación. La ecuación (1.13) da
paso a construir dos ecuaciones diferenciales ordinarias, una en función del tiempo y otra en fun-
ción del espacio tal que:

∂2 φ
= −λφ (1,14)
∂x 2
dG
= −λkG (1,15)
dt
Las soluciones producto (1.11) también deben cumplir las dos condiciones de contorno hom-
géneas. Asi, u(0, t ) = 0 implica que φ(0)G(t ) = 0, en este caso existen solo dos posibilidades: o
G(t ) ≡ 0 o bien φ(0) = 0. Si G(t ) ≡ 0, entonces por (1.11) la solución producto es idénticamente
cero (u(x, t ) ≡ 0). Este caso no nos interesa porque u(x, t ) ≡ 0 es la solución trivial ya que cumple
con las condiciones de contorno y la ecuacion deferencial parcial homgéneas.

Para obtener las soluciones no triviales se debe tener que

φ(0) = 0 (1,16)

φ(L) = 0 (1,17)
Donde (1.17) viene de utilizar la otra condición de contorno u(L, t ) = 0.

Ahora pasaremos a resolver la EDO en función del tiempo. La ecuación dependiente del tiempo
no tieene condiciones adicionales, simplemente es

dG
= −λkG
dt
Así, esta ecuación (1.15) es una ecuación diferencial ordinaria de primer orden con coeficien-
tes constantes. Al encontrar su solución nos damos cuenta que es bastante fácil de resolver, ya
que casi todas las EDOs (lineales y homégeneas) con coeficientes constantes se pueden resolver
buscando soluciones exponenciales de la forma G = e r t . Sustituyendo en la ecuación el polinomio
carracterístico obtenemos que r = −λk. Entonces la solución general de (1.15) es

G(t ) = ce −λkt
Se ha tomado como hecho que en las ecuaciones lineales homogéneas si e −λkt es una solu-
ción, entonces ce −λkt es también una solución para cualquier constante arbitraria c.

Ahora veremos la parte dependiente de x de la supuesta solución producto φ(x) una ecuación
diferencial ordinaria de segundo orden con dos condiciones de contorno homogéneas:

d 2φ
= −λφ, (1,18)
d x2
φ(0) = 0,
φ(L) = 0.
1.4. ECUACIÓN DIFERENCIAL DE CALOR EN UN MEDIO ISOTRÓPICO 15

Donde lo llamaremos problema de contorno para una ecuación diferencial ordinaria. Al igual
que en el caso dependiente sólo del tiempo t , aqui podemos suponer que si φ(x) = 0, entonces sa-
tisface la EDO y las condiciones de contorno homogéneas sin importar el valor de la constante de
separación λ. A ésta la llamamos solución trivial del problema de contorno. Existen algunos valo-
res especiales de λ que llamaremos autovalores del problema de contorno para los cuales existen
soluciones no triviales φ(x), la cuales existen solamente para algunos valores de λ; a éstas las lla-
maremos autofunción correspondiente al autovalor λ.

Ahora se determinaran los autovalores λ, es decir, ¿Para que valores de λ existen soluciones no
triviales? Se puede resolver (1.18) directamete, la EDO se segundo orden es lineal, homogénea y
con coeficientes constantes. Se generarám dos soluciones independientes con forma exponencial
φ = e r x . Asi, sustituyendo esta exponencial en la ecuación diferencial se obtiene el polinomio ca-
racterístico r 2 = −λ. Dependiendo del valor de λ las soluciones correspondientes a las dos raíces
tienen propiedades muy diferentes. Por ello veremos los siguientes casos:

1. λ > 0, aqui tenenmos

p que las dos raíces son imaginarias puras y complejas conjugadas una
de otra, r = + − i λ.

2. λ = 0 aqui las dos raíces se confunden en una, r = 0.

p
3. λ < 0 donde las dos raíces son reales y distintas, r = + − −λ.

4. λ es compleja.

Se ignorará el último caso ya que los autovalores deben ser reales para que existan soluciones
no triviales del problema de contorno (1.18). Ahora veremos el primer caso:

Autovalores y autofunciones (λ > 0). Asi se tiene el problema de contorno

d 2φ
= −λφ, (1,19)
d x2
φ(0) = 0,

φ(L) = 0.

Si λ > 0 las soluciones oscilan, ya que las soluciones exponenciales tienen componentes ima-
p
+−i λx
p p
ginariaspr = e . Como cos λx y sen λx son combinaciones lineales independientes de
r = e +−i λx , asi obtenemos la solución general de (1.19):
p p
φ = c 1 cos λx + c 2 sen λx

Al aplicar las condiciones de contorno: Si φ(0) = 0, entonces:

0 = c1 ,
p
es decir, se cancela el término del coseno cuando x = 0. Asi, se tiene que c 2 sen λx. Solo resta
aplicar la otra condición de contorno, φ(L) = 0, lo cual implica que
p
0 = c 2 sen λL.
p
Podemos observar que para que se verifique esta condición debe ser, o bien c 2 = 0 o sen λL =
0. Si c 2 = 0, entonces φ(x) ≡ 0 ya que previamente teníamos que c 1 = 0, la cual es la solución trivial,
pero nosotros estamos interesados en soluciones no triviales. Asi, los autovalores deben cumplir
16 CAPÍTULO 1. SERIES DE FOURIER

p
sen λL = 0
p
Gracias a estudios previos del comportamiento de la función seno tenemos que λL = nπ. Por
lo tanto, los autovalores λ son:
³ nπ ´2
λ= , n = 1, 2, 3, . . . (1,20)
L
Y la autofunción correspondiente sería
p nπx
φ(x) = c 2 sen λx = c 2 sen (1,21)
L