See

discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/281998659

DEDUCCIN DE LA ECUACIN DE CALOR POR

COORDENADAS RECTANGULARES Y
CURVILINEAS.

Article June 2013

CITATIONS READS

0 5,684

2 authors, including:

Fernando Zorto
National Autonomous University of Honduras
16 PUBLICATIONS 0 CITATIONS

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

Design a turbo-Jet Engine. View project

All content following this page was uploaded by Fernando Zorto on 22 September 2015.

The user has requested enhancement of the downloaded file.

 1

UNIVERSIDAD NACIONAL AUTONOMA DE HONDURAS

FACULTAD DE INGENIERIA
ESCUELA DE INGENIERIA MECNICA Y ELECTRICA
DEPARTAMENTO DE INGENIERIA MECNICA

Transferencia de Calor IM-414

Deduccin de la Ecuacin General de Calor por Conduccin

Coordenadas Rectangulares y Curvilneas

Ing. Jorge Gallo Navarro

ME., MSc., MBA.

Ing. Fernando Jos Zorto Aguilera

ME.

Resumen

En el Estudio de transferencia de calor, uno del mecanismo ms utilizado en la ingeniera y

como primer tema de cualquier curso de transferencia calor, se estudia a fondo lo que es la
transferencia de calor por conduccin, en este caso, analizamos un cuerpo rectangular para poder
definir la ecuacin general de calor en coordenadas cartesianas, es decir, que el vector del flujo de
calor ser ( ) y de este deduce dicha ecuacin general de calor por conduccin. Se plantea
tambin que en algunos casos especiales un espacio euclidiano es muy complicado para poder
resolver ciertos problemas en la transferencia por conduccin, en este caso, se expone la solucin
general para el cambio de coordenadas rectangulares a cualquier Espacio no Euclidiano en
coordenadas curvilneas.
 2

Introduccin:
En el presente trabajo se deducir la ecuacin general de calor por conduccin para
cualquier tipo de espacio vectorial, recordando, que el flujo de calor es un vector y que la
resolucin de dicha ecuacin ser determinada por la forma geomtrica del solido que transfiere
energa mediante la conduccin en su medio.

Objetivos:
1. Definir una ecuacin que generalice el fenmeno fsico de la transferencia de calor a travs de
la conduccin y los fenmenos adicionales de energa que ocurren en el cuerpo en
coordenadas rectangulares.

2. Definir una solucin general para la transformacin de espacio de uno euclidiano con base
cannica a otro espacio curvilneo ortonormal u ortogonal.

3. Definir la ecuacin que generalice el fenmeno fsico de la transferencia de calor a travs de la

conduccin y los fenmenos adicionales de energa que ocurren en un cuerpo curvilneo.

ANALISIS
Ecuacin General de Calor por Conduccin
La ecuacin de calor general en coordenadas rectangulares por conduccin resulta de un
balance de energa del vector de flujo calor, la energa generada y la energa que se almacena en
dicho sistema. Esto es expresado de la siguiente manera:

(1.1)

( ) ( ) ( ) (1.2.)

Aplicando la serie de Taylor tendremos que los son iguales a lo

indicado en la ecuacin 1.3, donde i representa cualquiera de las
variables independientes coordenadas, tal como se indica en las
ecuaciones 1.4 a la 1.6 respectivamente para x, y y z:
(1.3.)

Por lo tanto:
(1.4.)

(1.5.)

(1.6.)

Sustituyen las ecuaciones 1.3. a la 1.6 en 1.2 Tendremos que

( ) ( ) ( ) (1.7.)
 3

Acomodando los trminos, se tiene que:

( ) ( ) ( ) ( ) (1.8.)

Se cancelan los flujos de calor en la direccin x, y y z, con lo cual obtenemos la nueva ecuacin:

( ) (1.9.)

Al desarrollar un poco ms la igualdad anterior, se observa lo siguiente:

( ( ) ( ) ( ) (1.10.)

Ahora se debe recordar cmo se define el flujo de calor dado que en este caso ( ), es decir
que el flujo de calor es un vector en un espacio coordenado Cartesiano. Al encontrarse en un
espacio rectangular se dice que la base que define dicho espacio se encuentra conformada por los

vectores ( ) ( ) ( ), que se denominan base cannica. Planteando lo anterior el flujo de calor

por medio de la ley de Fourier:

( ) (1.11.)
Donde
(1.12.)
(1.13.)

(1.14.)

Retornemos a la ecuacin 1.9 y sustituir los flujos de calor respectivos de las ecuaciones de la 1.12
hasta 1.14 :
( ) ( ) ( ) (1.15.)

( ) ( ) ( ) (1.16.)

Pasamos a multiplicar el signo negativo a fuera de los parntesis:

( ) ( ) ( ) (1.17.)

En los siguientes pasos se encontrar la ecuacin general de calor por conduccin. Ahora se debe
determinar los diferenciales de rea y volumen involucrados en la ecuacin.

(1.18.)
(1.19.)
(1.20.)
(1.21.)
Hay que hacer una mencin de donde salen estos trminos de las diferenciales reas, es el
rea transversal del flujo de calor de un cubo infinitesimal y el dv se relaciona al volumen
 4

infinitesimal del mismo cubo. Sustituyendo las ecuaciones 1.18 al 1.20 en la 1.17 se obtienen la
ecuacin 1.22

( ) ( ) ( ) (1.22.)

Los diferenciales que se encuentran dentro del parntesis son constantes y se pueden pasar a
multiplicar a fuera del mismo y ser divididos a ambos extremos de la igualdad, lo que da como
resultado la eliminacin del elemento diferencial de volumen

( ) ( ) ( ) (1.23.)

( ) ( ) ( ) (1.24.)

Asumiendo una Conductividad Trmica constante, pasamos el trmino a fuera de los parntesis y
luego dividindolo a ambos lados de la igualdad, donde se tiene:

( ) ( ) ( ) (1.25.)

( ) ( ) ( ) (1.26.)
Y, recordando que:
( ) (1.27.)

Al sustituir la relacin 1.27 en 1.26 donde corresponda para la variable independiente respectiva,
se obtiene que:
(1.28.)

El trmino se conoce como , donde es la difusividad trmica que representa la velocidad

de propagacin del calor en el medio slido. Al realizar todas estas operaciones matemticas, la
expresin resultante se conoce como Ecuacin General de Conduccin de Calor, en el caso que
la conductividad trmica sea constante. La ecuacin 1.28 tambin se conoce como la Ecuacin de
Fourier-Biot y, en condiciones especficas, se reduce a estas tres formas conocidas como
Ecuacin de Poisson, Ecuacin de Difusin de Calor y Ecuacin de Laplace respectivamente:

1) Rgimen Estacionario: (Ecuacin de (1.29.)

Poisson )

(1.30.)
2) Rgimen transitorio sin generacin de calor:
(Ecuacin de difusin)

3) Rgimen estacionario y sin Generacin de calor: (1.31.)

(Ecuacin de Laplace)
 5

Ecuacin de Calor en Coordenadas Curvilneas

En los problemas reales de transferencia de calor existirn
diversas maneras o mtodos para poder darles alguna solucin,
esto conlleva que algunas formas geomtricas por donde se
manifiesta el fenmeno sea muy difcil analizarlas por coordenadas
rectangulares, lo cual hace mucho ms fcil su estudio si
transformamos el espacio por donde existe el fenmeno de
transferencia de calor. De esta manera, si se plantearn en otro
espacio coordenado ajeno a un Espacio Euclidiano normal a uno
no Euclidiano se logra generar una transformacin curvilnea entre
el espacio coordenado cartesiano al nuevo sistema curvilneo
ortonormal como son los ejemplos de los espacios esfricos y los
cilndricos (Tou, 2011).

As mismo, Esta transformacin al nuevo espacio coordenado nos

ayuda a expandir nuestra visin sobre el problema de transferencia
de calor. Adicionalmente, se dice que un espacio vectorial no
Euclidiano, donde un sistema coordenado curvilneo se encuentra
definido por la base conformada por los vectores , , caso
contrario, lo que ocurre en un sistema Cartesiano, donde, la base
cannica est definida por x, y, z, as como se observa en la figura
1.1, En este caso, los vectores curvilneos se pueden expresar de
la siguiente manera:

( ) (1.32.)

( ) ( ) ( ) (1.33.)
Figura 1.1 Espacio Vectoriales a)
(1.34.)
Rectangular b) Espacio Curvilneo

Donde son las coordenadas de nuestro nuevo espacio coordenado, as mismo, , ,

se conocen como factores escalares, lo cuales se deben encontrar para realizar la transformacin.

La pregunta fundamental es como definir o realizar una Transformacin de un Espacio a otro;

Segn Osizik (1993) se debe definir en primera instancia la distancia de un vector en la base
cannica, sea, en sus correspondientes coordenadas rectangulares. Esto se calcula segn
Grossman (2011) de la siguiente manera:
(1.35.)

Al determinar que el vector tiene caractersticas infinitesimales, dado que se mueve de un punto a
otro, esa distancia sufre una pequea variacin en su posicin, la cual matemticamente puede ser
expresada como una diferencial del vector V, esto tambin, diferenciar cada una de sus
componentes rectangulares y curvilneas:

( ) (1.36.)
Donde la nueva distancia ser:
(1.37.)

Ahora bien, si una componente se puede expresar en trminos de otras componentes del nuevo
espacio curvilneo, podemos decir que cada componente cartesiana puede ser expresada como
una diferencial total correspondiente a sus nuevas coordenadas :
( ) (1.38.)

( ) ( ) ( ) (1.39.)
 6

La diferencial total se determina de la siguiente manera.

(( ) ( ) ( )) (( ) (
) ( )) (( ))
) ( ) (
(1.40.)
Desarrollando la serie de Taylor para los primeros trminos, en forma genrica se tiene la igualdad
dada en la ecuacin 1.41:
(1.41.)
Y sustituyendo sta en la ecuacin 1.40, el vector queda:
(( ) ( ) ( ) ( ))
(( ) ( ) ( ) ( )) ((
) ( ) ( ) ( )) (1.42.)
Eliminando los primeros tres pares de trminos semejantes del respectivo vector, se obtiene:
( ) ( )
( ) (1.43.)

Reagrupando la expresin anterior como sumatorias tendremos lo siguiente:

( ) ( ) ( ) (1.44.)
Dado que:
( ) (1.45.)

( ) (1.46.)

( ) (1.47.)
Los diferenciales de x, y, z son vectores con coordenadas curvilneas as como lo muestra la
ecuacin 1.32 a la 1.34, por lo tanto, podemos sustituir estas expresiones en 1.36, as que,

(1.48.)

( ) ( ) ( ) (1.49.)

(1.50.)

( ) ( ) ( ) (1.51.)

(1.52.)

( ) ( ) ( ) (1.53.)
 7

Sustituyendo en la ecuacin 1.37

=( ) ( ) ( ) ( ) ( ) ( )

( ) ( ) ( ) (1.54.)

Agrupando trminos semejantes:

(( ) ( ) ( ) ) (( ) ( ) ( ) )

(( ) ( ) ( ) ) (1.55.)
Entonces:

(1.56.)

Dnde:

(( ) ( ) ( ) ) (1.57.)

(( ) ( ) ( ) ) (1.58.)

(( ) ( ) ( )) (1.59.)
De esta manera se encuentran los factores escalares para hacer la transformacin, ahora el nuevo
vector transformado tendr las siguientes coordenadas:
( ) (1.60.)
Coordenadas Cilndricas
Esta es la manera general para cambiar el espacio coordenado para un sistema ortonormal
curvilneo, ahora se analizara que es lo que ocurre en una figura geomtrica igual a un cilindro, las
coordenadas internamente de un vector en cualquier cilindro se definen de la siguiente manera:
( ) ( ) (1.61.)

Por lo tanto se debe definir como se expresan sus

coordenadas rectangulares a coordenadas cilndricas:
(1.62.)

Adicionalmente, se debe de recordar que analizamos este

caso para determinar el vector de flujo de calor en
coordenadas cilndricas en ese caso tenemos entonces lo
siguiente:
( ) (1.63.)
Si analizamos el flujo de calor en una seccin infinitesimal
de un cilindro, en el cual, existe flujo de calores entrantes,
como, salientes donde se gener calor y se almacena parte
de la energa, se podra deducir una Ecuacin que nos
ayude a definir la transferencia de calor en coordenadas
cilndricas. En este caso determinemos con un balance de
Figura 2.2 Espacio Cilndrico energa cual sera dicha ecuacin:

(1.64.)
 8

( ) ( ) ( ) (1.65.)

Aplicando la serie de Taylor tendremos que los son iguales a:

(1.66.)
Por lo tanto:
(1.67.)
(1.68.)
(1.69.)

Sustituyen las ecuaciones 1.66 hasta la 1.68 en 1.64 se tendr que.

( ) ( ) ( ) (1.70.)

Acomodando los trminos se tiene:

( ) ( ) ( ) ( ) (1.71.)

Se cancelan los flujos de calor de los primeros tres pares de trminos y se obtiene la nueva
ecuacin diferencial parcial:

(1.72.)

Al desarrollar un poco ms la igual anterior tendramos la expresin 1.73:

( ) ( ) ( ) (1.73.)

Ahora bien debemos determinar en este caso el vector del flujo de calor como funcin de las
coordenadas cilndricas:

( ) (1.74.)

Donde los diferenciales de rea y de volumen son los siguientes:

(1.75.)
(1.76.)
(1.77.)
(1.78.)
En este caso determinemos el cambio de variables, para ello se necesita deducir los Factores
escalares obtenidos anteriormente en ecuaciones 1.57 a la 1.58 respectivamente; de igual manera,
sustituyendo sus equivalentes en coordenadas cilndricas dadas en ecuacin 1.62, haciendo las
derivadas parciales respectivas y considerando la unidad trigonomtrica definida por ( )
( ) , se obtienen las ecuaciones 1.79 a la 1.81

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) ) ( )
( ) (1.79.)

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) ) ( )
( ) (1.80.)
 9

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) )
(1.81.)

Por tanto tenemos que:

(1.82.)
(1.83.)
(1.84.)
Sustituyendo en la ecuacin 1.11, obtenemos que:

( ) (1.85.)

Dnde:
(1.86.)

(1.87.)

(1.88.)

(1.89.)
Por tanto:
( ) (1.90.)

Dnde:

; ; (1.91.)

Este sera el vector del flujo de calor segn la ley de Fourier en coordenadas Cilndricas ahora
sustituyendo las ecuaciones 1.90 en la ecuacin 1.72:

( ) ( ) ( ) (1.92.)

Sacando los trminos constantes de las respectivas derivadas parciales:

( ) ( ) ( ) (1.93.)

Sustituyendo la ecuacin 1.89 en la 1.92:

( ) ( ) ( ) ) ( ) (
(1.94.)
Dividiendo por el trmino que corresponde al diferencial de volumen y eliminando el
producto de los signos menos, se obtiene:

( ) ( ) ( ) (1.95.)

Aunque la conductividad trmica no es constante, para el caso ideal y deduccin de la ecuacin de

transporte de calor, se asumir constante, con lo que se obtiene la ecuacin 1.96:

( ) ( ) ( ) (1.96.)
 10

Esta expresin es la Ecuacin General de Calor por conduccin en coordenadas cilndricas.

As como las coordenadas Rectangulares vamos otras variaciones para este determinado espacio
curvilneo:

1. Rgimen Estacionario:
(Ecuacin de Poisson en ( ) ( ) ( ) (1.97.)
coordenadas cilndricas)
2. Rgimen transitorio sin
generacin de calor: ( ) ( ) ( ) (1.98.)
(Ecuacin de difusin
coordenadas cilndricas)
3. Rgimen estacionario y sin
Generacin de calor: (Ecuacin ( ) ( ) ( ) (1.99.)
de Laplace en coordenadas
cilndricas)

Coordenadas Esfricas
Ya analizado el mtodo de transformacin de espacio podemos analizar un nuevo conjunto de
coordenadas curvilneas, ahora se realizar el en una figura geomtrica igual a una esfrica, las
coordenadas internamente de un vector en cualquier cilindro se definen de la siguiente manera:

( ) ( ) (1.100.)

Por lo tanto se debe definir como se expresan sus

coordenadas rectangulares a coordenadas cilndricas:

(1.101.)

Adicionalmente, se debe de recordar que analizamos este

caso para determinar el vector de flujo de calor en
coordenadas cilndricas en ese caso tenemos entonces lo
siguiente:
( ) (1.102.)

Si analizamos el flujo de calor en una seccin infinitesimal de

un cilindro, en el cual, existe flujo de calores entrantes,
como, salientes donde se gener calor y se almacena parte
de la energa, se podra deducir una Ecuacin que nos ayude
Figura 3.3 Espacio Esfrico a definir la transferencia de calor en coordenadas cilndricas.

En este caso determinemos con un balance de energa cual sera dicha ecuacin:

(1.103.)

( ) ( ) ( ) (1.104.)

Aplicando la serie de Taylor tendremos que los son iguales a:

(1.105.)
 11

Por lo tanto:
(1.106.)

(1.107.)

(1.108.)

Sustituyen las ecuaciones 1.106 a la 1.108 en 1.104 y agrupando loe trminos, se obtiene:

( ) ( ) ( ) ( ) (1.109.)

Cancelando los primeros tres pares de trminos correspondientes al flujo de calor, se obtiene la
nueva ecuacin:

( ) (1.110.)

Al desarrollar un poco ms la igual anterior tendramos lo siguiente:

* ( ) ( ) ( ) (1.111.)

Ahora bien debemos de determinar en este caso el vector del flujo de calor como funcin de las
coordenadas cilndricas
( ) (1.112.)
Donde los diferenciales de rea y de volumen son los siguientes:
(1.113.)
(1.114.)
(1.115.)
(1.116.)

De igual manera al caso cilndrica, se necesita determinar el cambio de variables, para esto
necesitamos deducir los Factores escalares obtenidos anteriormente en ecuaciones 1.57 a la 1.58
respectivamente. Sustituyendo sus equivalentes en coordenadas esfricas dadas en ecuacin
1.101, haciendo las derivadas parciales respectivas y considerando la unidad trigonomtrica
definida por ( ) ( ) ( ) , se obtienen las ecuaciones 1.117 a la
1.119

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) )
( ) ( ) ( ) (1.117.)

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) )

( ) ( ) (1.118.)

(( ) ( ) ( ) ) (( ( )) ( ( )) ( ( )) )
( ) ( ) ( ) (1.119.)
 12

Considerando los respectivos diferenciales se tiene:

(1.120.)
(1.121.)
(1.122.)

Sustituyendo en la 1.112 las ecuaciones 1.120 a la 1.122:

( ) (1.123.)
Dnde:
(1.124.)
(1.125.)
(1.126.)
(1.127.)

Por tanto sustituyendo en 1.123 las ecuaciones 1.124 hasta la 1.127:

( ) (1.128.)
Dnde:

; ; (1.129.)

Este sera el vector del flujo de calor segn la ley de Fourier en coordenadas Esfricas ahora
sustituymoslo en la ecuacin 1.110:
( ) ( ) ( )
(1.130.)
( ) ( ) ( )
(1.131.)

Sustituyendo la ecuacin 1.127 en la 1.131

( ) ( ) ( )
(1.132.)

Dividiendo la ecuacin 1.132 por el trmino , se obtiene:

( ) ( ) ( ) (1.133.)

De igual manera que en los dos casos anteriores, se asume una conductividad trmica constante y
eliminando el signo menos por el producto de los dos signos menos, se tiene que:

( ) ( )
( ) ( ) (1.134.)
 13

Esta expresin es la Ecuacin General de Calor por conduccin en coordenadas esfricas.

As como las coordenadas Rectangulares vamos otras variaciones para este determinado espacio
curvilneo:
1. Rgimen Estacionario:
(Ecuacin de Poisson
en coordenadas ( ) ( )
( ) ( )
esfricas)
(1.135.)
2. Rgimen transitorio sin
generacin de calor:
(Ecuacin de difusin ( ) ( )
( ) ( )
coordenadas esfricas)
(1.136.)
3. Rgimen estacionario y sin
Generacin de calor:
(Ecuacin de Laplace en
coordenadas esfricas) ( ) ( )
( ) ( )
(1.137.)

Conclusiones

1. En el balance de Energa que existe en un cuerpo multidimensional, es afectado por un flujo de

calor incidiendo en todas sus dimensiones; as mismo, se le aplica una generacin calor y
ocurre una almacenamiento de energa, se logra determinar que la ecuacin de calor por
conduccin es una ecuacin diferencial parcial; cuya solucin ser una funcin matemtica de
sus coordenadas en el espacio vectorial que sea estudiado. Adicionalmente esta solucin
determina la distribucin de Temperatura en cualquier punto del espacio vectorial en el estudio.
En este caso, se deben plantear las condiciones de frontera o las condiciones fsicas para
poder determinar la solucin ms adecuada para dicha ecuacin.

2. Se puede realizar un cambio de espacio vectorial, con el simple hecho de realizar un cambio
de base entre 2 espacios ortonormales u ortogonales, en este caso, trasladar un Espacio
Euclidiano rectangular con base cannica a otro espacio vectorial no Euclidiano con base
curvilnea. Esto se logra hacer simplemente observando la variacin existente en
desplazamiento de posicin de un vector, esto se define como una diferencial total a las
coordenadas curvilneas que definen a las coordenadas rectangulares en el espacio
rectangular. Por tanto, se puede concluir que la trasformacin de un espacio a otro es una
linealizacin del mismo con un factor de escala determinado, es decir que los trminos son los
factores escalares a1,a2,a3, los cuales se multiplican a los vectores de la base cannica, lo que
conlleva que la base nueva es una base linealmente independiente

3. La transferencia de calor para un espacio coordenado curvilneo puede tambin encontrarse

una expresin matemtica que modele dicho mecanismo de intercambio de calor, haciendo las
respectivas transformaciones de las coordenadas del espacio, se puede determinar cualquier
ecuacin general de calor para un sinfn de espacios curvilneos, en nuestro caso los
determinados en coordenadas Cilndricas y esfricas.

Bibliografa
1. Cengel, Y. A. (2007). Transferencia de Calor y Masa, un enfoque Prctico. Mexico D.F.:
Mcgraw Hill.
2. Grossman, S. (2011). Algebra lineal. Mexico DF: McGraw-Hill.
3. Holman, J. P. (2010). Heat Transfer. New York: Mcgraw-Hill.
4. Osizik, N. (1993). Heat conduction . New York: John Wiley and Sons.
5. Prieve, D. C. (2000). A course in Fluid Mechanics with Vector Field Theory. Pittburgh:
Carnegie Mellon University.
6. Tou, S. (2011). Visualization of Fields and Applications in Engineering. John Wiley & Sons.

View publication stats