Análisis de Fourier aplicado a la resolución de la Ecuación del Calor

Leandro G. Ercoli

Estudiante de Ingeniería en Sistemas de Computación

Universidad Nacional del Sur, Avda. Alem 1253, B8000CPB Bahía Blanca, Argentina
leandroercoli@[Link]
Febrero 2014

Resumen: En este trabajo, se presenta la resolución de la ecuación del calor unidimensional mediante la aplicación de las
series de Fourier primero, haciendo uso del método de separación de variables para resolver ecuaciones diferenciales y
luego por medio de transformadas de Fourier. Los resultados obtenidos en el primer caso, se aplican a un ejemplo
concreto del enfriamiento de una barra de sección constante.
Palabras clave: Serie de Fourier, Transformada de Fourier, Ecuación del calor.

I. INTRODUCCIÓN
En física, se define al “flujo de calor” como el proceso de transferencia de energía debido a la diferencia
de temperatura entre dos puntos. El calor es una forma de energía y su unidad en el sistema SI es el Joule (J).
Por su parte, la “temperatura” de una sustancia es una medida de la energía cinética de sus átomos. Según las
definiciones de Isaac Newton “una diferencia de temperatura entre dos sustancias en contacto, causa una
transferencia de energía desde la sustancia con mayor temperatura a la sustancia con menor temperatura, y se
denomina flujo de calor. Si no hay diferencia de temperatura, no hay transferencia de calor”. [1]
A principios del siglo XIX, Joseph Fourier (científico y matemático francés), introdujo la idea de
expandir funciones en series trigonométricas como método para resolver la ecuación de conducción del calor
desarrollada en su “Théorie Analytique de la Chaleur” (Teoría analítica del calor). Su texto clásico se
convirtió desde ese entonces, en la fuente para los métodos modernos de resolución de problemas prácticos
asociados con ecuaciones diferenciales parciales sujetas a condiciones de frontera dadas. [2]

II. FUNDAMENTOS TEÓRICOS

Una serie de Fourier es una expansión de una función periódica f(t) de periodo 2 / en la que el
conjunto base es el conjunto de funciones seno de diferentes amplitudes, fases y periodos, obteniéndose una
representación expandida como:

1
~ cos , (1)
2

donde los coeficientes están dados por las fórmulas de Euler:

2
cos 0, 1, 2, …

2
sen 1, 2, 3, …
Las ideas y técnicas de las series de Fourier se extienden a funciones f(x) no periódicas definidas en la
recta real completa; esto lleva a la integral de Fourier:
1
(2)
2
donde:
(3)
es la transformada de Fourier . Entonces (2) es la transformada inversa de Fourier.
III. SOLUCIÓN DE LA ECUACIÓN DEL CALOR POR SERIES DE FOURIER
Un ejemplo típico de transferencia de calor es encontrar la distribución de temperatura transitoria
unidimensional u(x,t) de una barra de sección constante y longitud L, con una determinada distribución inicial
de temperatura , cuyos extremos (x=0, x=L) se mantienen a temperatura constante. El problema se
resuelve a partir de la ecuación general del calor que en dirección x es:

(4)

donde K: conductividad térmica, : calor específico, : densidad. Suponiendo que los extremos de la barra se
mantienen a 0ºC, se tienen las siguientes condiciones de frontera:
0, 0, , 0 0 (5)
En general, si la temperatura inicial de la barra está representada por una función f(x), la condición
inicial será:
,0 ,0 (6)
Aplicando el método de separación de variables o método del producto a la ecuación de calor (4), se
determinan sus soluciones de la forma , ; que son producto de dos funciones, cada una de
las cuales depende únicamente de una de las variables (x y t). Reemplazando en (4) se obtiene:

" ,

que se puede reescribir así:

"
(7)

Con las variables separadas, el lado izquierdo de (7) es una función sólo de la variable independiente x, y
el lado derecho es una función sólo de la variable independiente t, entonces al cambiar x o t se afectaría
únicamente el primer miembro o el segundo, respectivamente (sin alterar el otro miembro), por lo tanto
ambos miembros deben ser iguales a una constante k. Para k≥0 la solución única u=XT que satisface (5) es
u≡0 (trivial). Para (negativa) se obtienen de (7) dos ecuaciones diferenciales lineales ordinarias:
" 0 (8)
0 (9)

Una solución general de (8) es:

(10)
Las soluciones X y T de (8) y (9) deben satisfacer las condiciones de frontera (5), es decir: 0,
0 0, , 0 para toda 0. Si T≡0, u≡0, entonces 0 0y
por (10) se obtiene 0 0 y entonces:
0, (11)
con 0 (a fin de evitar la solución trivial 0). Las soluciones no triviales se obtienen sólo si la
constante toma ciertos valores. Por (11), y como , las soluciones para son:
, 1, 2, 3, …
Se obtienen así las soluciones no triviales para X:
nπx
sin , n 1, 2, … B una constante arbitraria
L
Se resuelve ahora la otra ecuación diferencial (9), que para es: 0, donde .
Las soluciones a esta ecuación son: , 1, 2, … ; donde es una constante. Por lo
tanto las funciones:
 , 1, 2, … (12)

son soluciones de la ecuación de calor (4) que satisfacen las condiciones de frontera (5). Estas funciones son
las llamadas eigenfunciones o funciones características del problema, que corresponden a los eigenvalores o
valores característicos √ / . Para obtener una solución que también satisfaga la condición inicial
(6), se considera una serie de estas eigenfunciones (12):

, , (13)

A partir de esta expresión y de (6) se tiene:

,0

Por lo tanto los términos deben elegirse de modo que , 0 quede como la serie de senos de Fourier
de ; así:
2
1, 2, … (14)

Suponiendo que es continua por secciones en el intervalo 0 y que tiene derivadas laterales
en todos los puntos interiores de ese intervalo, la serie de Fourier (13) con coeficientes (14) es la solución del
problema físico planteado. Debido al factor exponencial, todos los términos de (13) tienden a cero cuando t
tiende a infinito; la rapidez del decremento varía con n.

IV. EJEMPLO USANDO LA SOLUCIÓN DE LA ECUACIÓN DEL CALOR POR SERIES DE FOURIER
Si en el caso particular de la barra bajo estudio la distribución de temperatura inicial está dada por:
0 /2
/2
se obtienen los coeficientes:
/
2
.
/

Si n es par, integrando se obtiene 0. Si n es impar, se tiene:

4 4
1, 5, 9, … 3, 7, 11, …

Por lo tanto, la solución es:

4 1 3 3
, . (15)
9
En la figura 1 se presentan los resultados de u(x,t) para distintos tiempos. Aquí puede verse que la
distribución de temperatura de la barra va cambiando hacia una forma más “aplanada” a medida que
transcurre el tiempo y se enfría como consecuencia de la pérdida de calor hacia los extremos. Finalmente,
cuando se alcance el estado estacionario (t>>2 s), toda la barra estará a la temperatura de los extremos (0 ºC).

V. SOLUCIÓN DE LA ECUACIÓN DEL CALOR POR TRANSFORMADA DE FOURIER

Ahora se resolverá la ecuación del calor, con las mismas condiciones iniciales y de frontera, haciendo uso
de las relaciones (2) y (3). Se buscan soluciones a (4) que tengan la forma , , tal que X(x) se
mantenga finito para | | ∞. Las soluciones generales que se ajustan a estos requerimientos son múltiplos
 Figura 1: Solución del ejemplo para , 1 y varios valores de t.

de exp , donde k es cualquier número real. Se construye por el método de superposición la

siguiente solución general:

, exp . (16)

En la solución por series de Fourier, las condiciones iniciales en x=0 y x=L requerían seleccionar sólo
una familia finita de números para k, pero ahora k puede tomar cualquier número real. Además se resolvió la
ecuación tomando una suma de eigenfunciones o funciones características; ahora esta suma se convierte en
una integral. Se puede satisfacer la condición inicial (6) por (3) y (2) si /2 . Introduciendo (3)
en (16), la solución se puede expresar en términos de la distribución inicial f(x), como
1
, exp .
2
Invirtiendo el orden de integración y haciendo la integración con respecto a k, se encuentra la solución:
1 /
, . (17)
√4

VI. CONCLUSIÓN
El análisis de Fourier resulta ser una importante herramienta en la resolución de problemas en muchas
ramas de la ingeniería (transferencia de calor y masa, análisis de ondas y señales, etc.). En particular se
estudió en este informe la aplicación de series y transformada de Fourier para encontrar la distribución de
temperatura no estacionaria unidimensional en una barra sometida a condiciones iniciales y de frontera
especificadas. Queda demostrado que las herramientas matemáticas estudiadas en esta asignatura son de
aplicación directa para la resolución de problemas ingenieriles reales.

REFERENCIAS
[1] B. Osgood, “Lecture Notes for EE 261, The Fourier Transform and its Applications”, Standford
University, disponible en [Link] [consultada el
3 de Febrero de 2014].
[2] G. James, Matemáticas Avanzadas para Ingeniería, Pearson Educación, 2002.
[3] F. P. Incropera, Fundamentals of Heat and Mass Transfer, J. Wiley & Sons, 2007.
[4] E. Kreyszig, Matemáticas Avanzadas para Ingeniería, vol II, Limusa Wiley, 2004.
[5] C. C. Lin & L. A. Segel, Mathematics Applied to Deterministic Problems in the Natural Sciences,
Macmillan Publishing, 1974.