UNIVERSIDAD AUTÓNOMA METROPOLITANA

Unidad Iztapalapa
Notas de curso

Transferencia de Calor

Francisco J. Valdés Parada

Transferencia de Calor
Notas de clase
Francisco J. Valdés Parada
Transferencia de Calor
Última actualización: 31 de julio de 2021

c Francisco J. Valdés Parada

Universidad Autónoma Metropilitana-Iztapalapa. T-261 Av. San Rafael Atlixco No. 186, Col.
Vicentina, Del. Iztapalapa, C. P. 09340, CDMX, México

e-mail: [email protected]
 Índice general

I Parte Uno: Fundamentos

1 Introducción: Repaso de mecánica de fluidos . . . . . . . . . . . . . . . . . 11

1.1 Generalidades sobre el modelado 12
1.2 Un primer vistazo a la transferencia de calor 13
1.3 Tres tipos de regiones y de derivadas temporales 14
1.4 Teoremas integrales 15
1.4.1 Teoremas del transporte y de la divergencia . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4.2 Teorema de la divergencia y fórmula de Green . . . . . . . . . . . . . . . . . . . . . . . . 15
1.5 Conservación de masa 16
1.6 Conservación de cantidad de movimiento lineal 18
1.7 Ecuación de energía mecánica 20
1.8 Resumen 22

2 Deducción de la ecuación de energía . . . . . . . . . . . . . . . . . . . . . . . . 23

2.1 Introducción 23
2.2 Primera ley de la termodinámica 23
2.3 Ecuación de energía total 26
2.4 Ecuación de energía térmica 27
2.5 Ley de Fourier 29
2.6 Mecanismos de transferencia de calor 32
2.6.1 Conducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
2.6.2 Convección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.6.3 Radiación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.7 Condiciones de frontera 34
2.7.1 Condición de salto del flux de calor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34
2.7.2 Condición de salto de la temperatura . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.7.3 Clasificación de condiciones de frontera desde el punto de vista matemático 37
2.8 Números adimensionales asociados a la transferencia de calor 39
2.9 Resumen 41

II Parte Dos: Aplicaciones

3 Conducción de calor en estado estacionario . . . . . . . . . . . . . . . . . . 45

3.1 Introducción 45
3.2 Dos alternativas para simplificar modelos 45
3.2.1 Simplificaciones basadas en estimados de orden de magnitud . . . . . . . . . . . 45
3.2.2 Simplificaciones basadas en un proceso de promediado . . . . . . . . . . . . . . . . 48
3.3 Conducción de calor en una placa 49
3.3.1 Temperaturas superficiales especificadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
3.3.2 Placa con generación de calor uniforme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.3.3 Transferencia de calor considerando convección externa . . . . . . . . . . . . . . . 53
3.4 Conducción unidimensional en estado estacionario en coordenadas cilín-
dricas 55
3.5 Conducción de calor en coordenadas esféricas 57
3.6 Una introducción al escalamiento: La aleta de enfriamiento 59
3.6.1 Análisis basado en suponer que la temperatura en la superficie es igual a la
temperatura promedio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3.6.2 Análisis basado en en la suposición de desviaciones periódicas . . . . . . . . . . 63
3.7 Aletas de enfriamiento cilíndricas 65
3.8 Aletas de enfriamiento anulares 66

4 Transferencia de calor en más de una dirección . . . . . . . . . . . . . . . 71

4.1 El problema de Sturm-Liouville 71
4.2 Coordenadas rectangulares 72
4.2.1 Transporte bidimensional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72
4.2.2 Transporte por conducción en tres dimensiones . . . . . . . . . . . . . . . . . . . . . . . 78
4.3 Conducción de calor transitoria en una barra cilíndrica 82
4.4 Transporte en medios semi-infinitos 87
4.4.1 Solución exacta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88
4.4.2 Flux de calor en la frontera . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4.4.3 Ley de enfriamiento en la frontera . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
4.4.4 Solución aproximada par el método de von Kármán-Pohlhausen . . . . . . . . . . 92
5 Soluciones numéricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5.1 Coordenadas cilíndricas 95
5.1.1 Transporte axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
5.1.2 Transporte no axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
5.2 Problema de aplicación 99
5.3 Aleta de enfriamiento revisada 99
5.4 Conducción de calor transitoria 100
5.5 Conducción y convección en estado transitorio 102

6 Capa límite . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107

6.1 Ecuaciones de capa límite de energía y cantidad de movimiento para
flujo laminar 107
6.1.1 Planteamiento y simplificación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
6.1.2 Solución aproximada de las ecuaciones de capa límite . . . . . . . . . . . . . . . . 110
6.2 Capa límite cuando se especifica el flux de calor en la superficie 114

7 Modelos macroscópicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 117

7.1 Introducción 117
7.2 Transferencia de calor en un panel fotovoltaico 117
7.3 Tanque agitado homogéneo 120
7.3.1 Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 120
7.4 Tanque agitado no homogéneo 122
7.5 Transferencia de calor en un lecho fijo 123
7.6 Transferencia de calor en intercambiadores de calor de doble tubo 128
7.7 Transferencia de calor en un pozo geotérmico 133

Bibliografía . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
I
Parte Uno: Fundamentos

1 Introducción: Repaso de mecánica de

fluidos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.1 Generalidades sobre el modelado
1.2 Un primer vistazo a la transferencia de calor
1.3 Tres tipos de regiones y de derivadas temporales
1.4 Teoremas integrales
1.5 Conservación de masa
1.6 Conservación de cantidad de movimiento lineal
1.7 Ecuación de energía mecánica
1.8 Resumen
2 Deducción de la ecuación de energía 23
2.1 Introducción
2.2 Primera ley de la termodinámica
2.3 Ecuación de energía total
2.4 Ecuación de energía térmica
2.5 Ley de Fourier
2.6 Mecanismos de transferencia de calor
2.7 Condiciones de frontera
2.8 Números adimensionales asociados a la transferen-
cia de calor
2.9 Resumen
 1. Introducción: Repaso de mecánica de fluidos

El transporte de energía es un tema crucial para cualquier alumno interesado en comprender

las ciencias básicas e ingeniería. Aunque el enfoque que tiene este texto está orientado hacia
la ingeniería, sus fundamentos recaen en la mecánica del medio continuo y las aplicaciones de
los conceptos vistos aquí trascienden los casos meramente industriales clásicos. De los textos de
termodinámica se sabe que la energía de un sistema se transfiere a través de sus interacciones
con los alrededores en forma de calor y trabajo. Sin embrago, en los textos introductorios de
termodinámica el interés está en los estados del sistema sin prestar atención al tipo y naturaleza de
las interacciones. El propósito de este texto es precisamente comprender estas interacciones, en
particular en la conversión de energía térmica (o calor) a energía interna de un sistema.
Los alumnos que toman este curso normalmente ya han tenido una primera experiencia con
los fenómenos de transporte mediante el curso de Mecánica de Fluidos. En el presente texto se
seguirá un esquema similar al utilizado en Mecánica de Fluidos, esto es: Primeramente se deducirán
cuidadosamente las ecuaciones que gobiernan el transporte de energía y las correspondientes condi-
ciones de frontera. Posteriormente, se estudiarán problemas con diferentes grados de complejidad,
empezando con problemas meramente conductivos en estado estacionario (ecuaciones ordinarias),
seguido de problemas en estado transitorio (ecuaciones parciales). Posteriormente, se acoplará la
transferencia de calor con el transporte de cantidad de movimiento en problemas de capa límite,
para después discutir acerca de la transferencia de calor en régimen turbulento. Finalmente, se
presentará un esquema de promediado de las ecuaciones a la escala puntual (o microscópica)
para así desarrollar las ecuaciones válidas a un nivel de escala macroscópico. Las herramientas
matemáticas para cumplir este objetivo son las mismas que las que se requieren para deducir las
ecuaciones a escala microscópica. Esto es, los teoremas del transporte y de la divergencia.
Para un mejor aprovechamiento del texto es recomendable que los alumnos cuenten con bases
sólidas en la solución (analítica y numérica) de ecuaciones diferenciales ordinarias y parciales,
análisis vectorial y el uso de software como Comsol Multiphysics. Si los alumnos no cuentan con
dichas bases, se recomienda que el profesor dedique algunas sesiones a remediar esta situación.
El resto de este capítulo está dedicado a repasar la deducción de las ecuaciones gobernantes
de la transferencia total de masa y de cantidad de movimiento en una fase y en una interfase.
 12 Capítulo 1. Introducción: Repaso de mecánica de fluidos

a) Matemático
b) Determinista
d) c) Estocástico
d) Fundamental
b) e) Observacional
∮𝛂d𝝈
a) c) e)

Modelado
ii)
i) Fı́sico
iii) ii) Operación
i) iii) Diseño
iv) Rediseño
iv)

Figura 1.1: Esquema de la clasificación del modelado. Las flechas amarillas llevan al tipo de modelo
a utilizar en este texto.

El Propósito es revisar los pasos involucrados en la deducción de estas ecuaciones así como el
herramental matemático necesario. Estos elementos son de especial importancia para la deducción
de la ecuación de energía, como se discute en el Capítulo 2.

1.1 Generalidades sobre el modelado

Un modelo puede entenderse, en general, como una representación de la realidad, la cual puede
llevarse a cabo de manera experimental (en cuyo caso se habla de modelos físicos) o teórica (que
para los propósitos de este texto corresponde a los modelos matemáticos) como se muestra en la
figura 1.1 (Herrera y Pinder, 2012). Un modelo físico puede definirse como una representación de
la realidad por medios experimentales. El trabajo experimental puede, a grandes rasgos, clasificarse
en tres actividades: operación (usar un equipo existente), diseño (construir un equipo o técnica
experimental) y rediseño (mejoramiento de una técnica o equipo ya existente).
Por su parte, los modelos matemáticos pueden definirse como una representación matemática de
aspectos no matemáticos (Aris, 1994). Como primera clasificación, se pueden distinguir dos tipos
de modelos matemáticos: los deterministas, en los cuales las mismas entradas producen siempre
las mismas salidas y los estocásticos donde lo anterior no siempre se cumple y el resultado está
sujeto a probabilidades (Himmelblau y Bischoff, 1968). De acuerdo al origen de su formulación,
los modelos matemáticos deterministas se clasifican en: los modelos basados en la observación y
los modelos basados en principios fundamentales. Los primeros (también conocidos como modelos
empíricos) necesitan de modelos experimentales que aporten datos confiables para posteriormente
encontrar expresiones matemáticas que los representen (Gershenfeld, 1998). En este caso, un
modelo será más aceptable que otro en función de qué tan cerca se ajusten las predicciones teóricas
a los datos experimentales. La formulación de los modelos basados en principios fundamentales,
 1.2 Un primer vistazo a la transferencia de calor 13

para los casos que se abordan en este texto, tienen su origen en la termodinámica de los sistemas
fuera de equilibrio (Groot y Mazur, 2011). Este tipo de modelos son predictivos, por lo que tras
ser resueltos, las predicciones que aportan deben validarse al compararse con datos obtenidos
experimentalmente. En este caso un modelo será mejor que otro no solo en la medida en que pueda
reproducir los resultados experimentales sino en cuanto a los fundamentos teóricos propuestos que
describan la situación física de interés. En síntesis, los modelos matemáticos se complementan
con los modelos físicos, de manera que los modelos basados en la observación, son modelos a
posteriori y los modelos basados en principios fundamentales son modelos a priori de la actividad
experimental.

1.2 Un primer vistazo a la transferencia de calor

La transferencia de calor corresponde a una porción de la primera ley de la termodinámica que
describe la conversión de calor a energía interna. A esta porción de la primera ley se le conoce
como la ecuación de energía térmica. La otra parte de la primera ley corresponde a la ecuación de
energía mecánica, la cual se estudia en los cursos de mecánica de fluidos y describe la conversión
de trabajo en energía cinética y potencial. De esta forma, la suma de las ecuaciones de energía
térmica y energía mecánica dan como resultado la ecuación de energía total, que no es otra que la
primera ley de la termodinámica y enuncia que la suma de las energías interna, cinética y potencial
es igual a la suma del calor y el trabajo de un sistema.
A grandes rasgos, la ecuación de energía térmica se compone de cuatro partes: la acumulación
de calor, el transporte por convección, el transporte por conducción y la fuentes (o sumideros) de
calor. La acumulación de calor es proporcional al cambio de la temperatura con el tiempo respecto
a un determinado punto en el espacio. El transporte por convección, se refiere a la transferencia
de calor debida al movimiento global de un fluido. Evidentemente, este modo de transferencia
de calor no está presente en sólidos ni en fluidos en condiciones hidrostáticas. En contraste, el
transporte conducción se da en cualquier fase y está asociado a las vibraciones de las moléculas.
Este tipo de transporte es proporcional a la segunda derivada espacial de la temperatura. Por último,
las fuentes o sumideros en la ecuación de energía térmica incluyen contribuciones debidas a la
radiación térmica, disipación viscosa, reacciones químicas, entre otras. Como puede notarse, los
tres modos de transferencia de calor (convección, conducción y radiación) están presentes en la
ecuación de energía térmica.
La parte introductoria de este texto trata sobre un repaso de los fundamentos de la mecánica
de fluidos y el herramental matemático indispensable. Posteriormente, se enfoca la atención en la
deducción de la ecuación de energía térmica en una versión microscópica que aplica a cada punto de
una fase. Una vez deducida esta ecuación, se obtendrá la condición de salto de la densidad de flujo de
calor entre fases y se postulará a la ley de enfriamiento de Newton como su complemento. Con estos
elementos a la mano, se pueden estudiar una amplia variedad de problemas de transferencia de calor
en estado estacionario y en estado transitorio. Para lo cual se usarán tanto soluciones analíticas como
numéricas de ecuaciones diferenciales ordinarias y parciales. A grandes rasgos, el procedimiento a
seguir en la solución de problemas de transferencia de calor es el siguiente: 1) Comprensión de la
situación física, que permita identificar potenciales simplificaciones y suposiciones; 2) Usar, si es el
caso, suposiciones y restricciones de escala (espacial y/o temporal) para simplificar la ecuación de
energía térmica y las condiciones de frontera correspondientes; 3) Resolver el modelo simplificado
de manera analítica o numérica y 4) Analizar los resultados y discutir acerca de los alcances y
limitaciones del modelo.
Al igual que en los cursos de mecánica de fluidos, este texto abarca temas como el transporte
de capa límite, el transporte turbulento de calor y los modelos macoscópicos de transferencia de
calor. El tema de capa límite servirá para deducir la ley de enfriamiento de Newton, mientras que
los temas de transporte turbulento y modelado macroscópico sirven de vínculo con una amplia
 14 Capítulo 1. Introducción: Repaso de mecánica de fluidos

variedad de aplicaciones prácticas de interés. Además de estos temas, se incluyen algunos tópicos
relevantes como el transporte de calor en medios porosos y el transporte en sistemas energéticos
tales como pozos geotérmicos y paneles fotovoltaicos.

1.3 Tres tipos de regiones y de derivadas temporales

Para los propósitos de este texto, es conveniente adoptar la hipótesis del continuo. Bajo esta
perspectiva, la materia se considera que está constituida por fases, las cuales están separadas entre
sí por interfases. Esta postura es una simplificación de la realidad, donde se sabe que dichas fases
se componen de átomos y moléculas, los cuales están separados por espacio vacío. Para describir lo
que ocurre en una o varias fases se pueden definir tres tipos de regiones espaciales de estudio:
Definición 1.3.1 — Regiones espaciales.
1. Región fija (Euleriana): Su forma y tamaño son fijos en el espacio y en el tiempo y
permite entradas y salidas de materia a través de su superficie.
2. Región material (Lagrangiana): A pesar de que su forma y tamaño pueden variar en
la posición y el tiempo, tiene la restricción de que siempre debe contener a la misma
cantidad de materia. Su movimiento corresponde a la velocidad del fluido, v.
3. Región arbitraria (general): No tiene restricciones respecto a su forma y tamaño, su
movimiento es independiente y su velocidad se denota como w. Permite entrada y salida
de materia a través de sus fronteras.

Estas regiones corresponden a tres tipos de observadores imaginarios que podrían considerarse
para estudiar un determinado proceso. El uso de estos tres t Asociadas a estas tres regiones, pueden
definirse los siguientes tres tipos de derivadas temporales:
Definición 1.3.2 — Derivadas temporales.
1. Derivada parcial (Euleriana): Este tipo de derivada consiste en medir los cambios
temporales de una propiedad, f , desde un punto fijo en el espacio. Usando al vector r
para representar la posición de dicho punto, se tiene entonces que

df ∂f
= (1.3.1)
dt r ∂t

2. Derivada material (Lagrangiana): Este tipo de derivada está asociada a la región ma-
terial, por lo que la posición desde la que se mide a la propiedad f se mueve con la
velocidad del fluido. Bajo esta perspectiva, la propiedad es función de la posición y del
tiempo, por lo que usando la regla de la cadena se tiene que
i=3
Df df dri df ∂f
= +∑ = +v·∇f (1.3.2)
Dt dt r i=1 dt t dri t ∂t

Aquí se usó la definición v = dr

dt t .
3. Derivada total (general): Esta derivada está asociada a una región arbitraria, la cual
se mueve a una velocidad w. El desarrollo es similar al de la derivada material, con
la diferencia que la derivada temporal del vector de posición corresponde ahora a la
velocidad arbitraria del observador, w, esto es:
i=3
df df dri df ∂f
= +∑ = +w·∇f (1.3.3)
dt dt r i=1 dt t dri t ∂t
 1.4 Teoremas integrales 15

1.4 Teoremas integrales

1.4.1 Teoremas del transporte y de la divergencia
En los párrafos anteriores se definieron derivadas asociadas a propiedades definidas en cada
punto de un sistema (es decir, propiedades puntuales); sin embargo, a menudo el interés está en
calcular derivadas temporales de integrales volumétricas de propiedades puntuales. El resultado de
tomar la derivada total de una integral en una región arbitraria constituye el teorema general del
transporte (la deducción de este teorema puede consultarse en los textos como los de Aris, 1990;
Gray y Gray, 2017)

Teorema 1.4.1 — Teorema general del transporte. Sea f una función continua en el espacio
y en el tiempo, la derivada total de la integral volumétrica de f en una región arbitraria satisface
la siguiente igualdad

d ∂f
Z Z Z
f (r,t) dV = dV + n · ( f w) dA (1.4.1)
dt ∂t
VA VA AA

Ahora bien, cuando se desea tomar la derivada material de una integral volumétrica correspon-
diente a una región material, el resultado es el teorema del transporte de Reynolds, el cual puede
expresarse como:

Teorema 1.4.2 — Teorema del transporte de Reynolds. Sea f una función continua en el
espacio y en el tiempo, la derivada material de la integral volumétrica de f en una región material
satisface la siguiente igualdad

D ∂f
Z Z Z
f (r,t) dV = dV + n · ( f v) dA (1.4.2)
Dt ∂t
VM VM AM

Por último, en el caso de la derivada parcial de una integral volumétrica correspondiente a una
región fija, el resultado es simplemente la siguiente identidad

∂f
Z Z
∂
f (r,t) dV = dV (1.4.3)
∂t ∂t
Vf Vf

Como puede notarse, las ecuaciones (1.4.2) y (1.4.3) pueden deducirse a partir de la ecuación
(1.4.1) fijando w = v y w = 0, respectivamente. De la misma forma, Las definiciones de la derivada
material (ecuación 1.3.2) y de la derivada parcial (ecuación 1.3.1) pueden deducirse a partir de la
definición de la derivada total dada en la ecuación (1.3.3) al fijar w = v y w = 0, respectivamente.

1.4.2 Teorema de la divergencia y fórmula de Green

Al comparar los teoremas general del transporte y del transporte de Reynolds con la identidad
dada en la ecuación (1.4.3) pareciera que las integrales de superficies en ambos teoremas corrigen,
por así decirlo, a dicha identidad al tomar en cuenta el desplazamiento de las regiones arbitraria
y material, respectivamente. Para el desarrollo de las ecuaciones de transporte, es conveniente
expresar dichas integrales de superficie como integrales volumétricas, para ello se recurre al teorema
de la divergencia o teorema de Gauss (1813), el cual se enuncia como sigue:

Teorema 1.4.3 — Teorema de la divergencia. Sea a una función vectorial diferenciable en

 16 Capítulo 1. Introducción: Repaso de mecánica de fluidos

el espacio, la integral volumétrica de su divergencia satisface la siguiente igualdad

Z Z
∇ · a dV = n · a dA (1.4.4a)
V A

Aquí la región V (de superficie A ) se refiere a una región arbitraria, material o fija. El uso de
este teorema puede extenderse al operador gradiente y rotacional como sigue
Z Z
∇a dV = na dA (1.4.4b)
V A
Z Z
∇ × a dV = n × a dA (1.4.4c)
V A

La ecuación (1.4.4b) puede aplicarse también para el caso en que la función sea un escalar.

Para cerrar esta sección, es conveniente analizar el caso en el que el teorema de la divergencia
se aplica al producto φ ∇ϕ, siendo φ y ϕ funciones escalares continuas en la posición. Sustituyendo
a = φ ∇ϕ en ambos lados de la ecuación (1.4.4a) se obtiene
Z Z Z
2


∇ · (φ ∇ϕ) dV = φ ∇ ϕ + ∇φ · ∇ϕ dV = n · (φ ∇ϕ) dA (1.4.5)
V V A

O bien, el resultado de tomar a = ϕ∇φ es

Z Z
2


ϕ∇ φ + ∇φ · ∇ϕ dV = n · (ϕ∇φ ) dA (1.4.6)
V A

Al restar la ecuación (1.4.6) a la ecuación (1.4.5) se obtiene la fórmula de Green:

Teorema 1.4.4 — Fórmula de Green. Sean φ y ϕ funciones escalares continuas con segundas
derivadas espaciales definidas, entonces se cumple que
Z Z
2 2


φ ∇ ϕ − ϕ∇ φ dV = n · (φ ∇ϕ − ϕ∇φ ) dA (1.4.7)
V A

Esta fórmula será de especial importancia en la solución de problemas de valor a la frontera

(Haberman, 2012) como se verá más adelante en el texto. De hecho, Valdés-Parada y col. (2007)
mostraron que con esta fórmula se pueden obtener soluciones numéricas de problemas no lineales.
Además, esta fórmula es el fundamento del método numérico del elemento a la frontera con el cual
se pueden resolver una amplia variedad de problemas de valor inicial y a la frontera (Katsikadelis,
2016). En los párrafos siguientes se revisa el uso del herramental matemático presentado arriba a la
deducción de las ecuaciones de conservación de masa, de cantidad de movimiento y de energía
mecánica en una fase.

1.5 Conservación de masa

El principio de conservación de masa es bien conocido y puede enunciarse en cualquiera de
las tres regiones mencionadas arriba. Como la región arbitraria es la más general, es conveniente
expresar en palabras este principio de conservación como sigue
1.5 Conservación de masa 17

Ecuación de continuidad (forma macroscópica con palabras)

( ) ( )
La razón de cambio en el tiempo Tasa neta de entrada de masa
= (1.5.1)
de la masa en la región por las fronteras de la región

El lado izquierdo de esta ecuación puede expresarse en lenguaje matemático como sigue
( )
La razón de cambio en el tiempo d
Z
= ρ dV (1.5.2)
de la masa en la región dt
VA

Note que la masa de la región se calculó sumando los productos de la densidad por el diferencial de
volumen a lo largo de VA . Por su parte, la tasa de entrada de masa por las fronteras se traduce en
lenguaje matemático mediante la siguiente integral de superficie
( )
Tasa neta de entrada de masa Z
= − n · ρ(v − w) dA (1.5.3)
por las fronteras de la región
AA

donde el signo negativo se debe a que el vector unitario apunta hacia afuera de la superficie AA .
Sustituyendo las ecuaciones (1.5.2) y (1.5.3) en la ecuación (1.5.1) da como resultado

Ecuación de continuidad (forma macroscópica)

d
Z Z
ρ dV + n · ρ(v − w) dA = 0 (1.5.4)
dt
VA AA

Esta forma del principio de conservación de masa es especialmente útil para estudiar problemas
macroscópicos. Sin embargo, por el momento el interés está en deducir su forma microscópica,
es decir una ecuación que sea válida en cada punto de una fase. Aplicando el teorema general del
transporte (ecuación 1.4.1) al primer término de la ecuación anterior lleva a la siguiente expresión
Z Z
∂ρ
dV + n · ρv dA = 0 (1.5.5)
∂t
VA AA

O bien, usando el teorema de la divergencia en el segundo término y agrupando los dos términos
resultantes en una sola integral
Z  
∂ρ
+ ∇ · (ρv) dV = 0 (1.5.6)
∂t
VA

Esta integral es satisfecha ya que el integrando es cero, por lo tanto se deduce de esta expresión la
ecuación de continuidad
Ecuación de continuidad (forma microscópica)

∂ρ
+ ∇ · (ρv) = 0 (1.5.7)
∂t
Esta ecuación es válida en cada punto de una fase y requiere además de una ecuación de estado
que relacione a la densidad con la presión y temperatura. Para casos en los que la densidad sea
constante en el espacio y en el tiempo, la ecuación anterior se deduce a
∇·v = 0 (1.5.8)
la cual expresa la condición de flujo incompresible.
 18 Capítulo 1. Introducción: Repaso de mecánica de fluidos

Ejercicio 1.1 Muestre que si se comienza por el principio de conservación de masa en una
región material o arbitraria, se puede llegar a la ecuación (1.5.5). 

Ejercicio 1.2 — Formulación alternativa de la ecuación de continuidad. Utilice la defi-

nición de derivada material y lleve a cabo las operaciones algebraicas necesarias para demostrar
que la ecuación de continuidad dada en la ecuación (1.5.7) puede expresarse como
Dρ
+ ρ (∇ · v) = 0 (1.5.9)
Dt


Ejercicio 1.3 — Teorema modificado de Reynolds. Utilice la ecuación de continuidad para

deducir el teorema modificado de Reynolds
D Dψ
Z Z
ρψ dV = ρ dV (1.5.10)
Dt Dt
VM VM

donde ψ es una función continua en el espacio y en el tiempo.



1.6 Conservación de cantidad de movimiento lineal

El principio de conservación de cantidad de movimiento lineal engloba a las tres leyes de

Newton y corresponde a la primera ley de Euler como lo explica (Whitaker, 2009). Para un dominio
arbitrario, este principio se puede expresar en palabras como (Gray y Gray, 2017)

Ecuación de movimiento lineal (forma macroscópica en palabras)

   
La razón de cambio en el tiempo de 
    Flujo neto de entrada de  
la cantidad de movimiento en = cantidad de movimiento por
 
  
un dominio arbitrario las fronteras del dominio
  
( ) ( )
Fuerzas que actúan Fuerzas que actúan en
+ + (1.6.1)
dentro del dominio las fronteras del dominio

Traduciendo al lenguaje matemático cada término de esta ecuación se obtiene

Ecuación de movimiento lineal (forma macroscópica)

d
Z Z Z Z
ρv dV + n · ρv(v − w) dA = ρb dV + tn dA (1.6.2)
dt
VA AA VA AA

donde b representa las fuerzas volumétricas por unidad de masa, mientras que tn es el vector de
fuerzas superficiales. De acuerdo al teorema fundamental de Cauchy, dicho vector se relaciona con
el tensor de esfuerzos totales como (Whitaker, 1992a):

tn = n · T (1.6.3)
1.6 Conservación de cantidad de movimiento lineal 19

Sustituyendo esta definición en la ecuación (1.6.2) y usando el teorema general del transporte en el
primer término del lado izquierdo se obtiene que
Z Z Z Z
∂ ρv
dV + n · ρvv dA = ρb dV + n · T dA (1.6.4)
∂t
VA AA VA AA

O bien, tras aplicar el teorema de la divergencia y agrupar todos los términos dentro de una sola
integral volumétrica:
Z  
∂ ρv
+ ∇ · (ρvv) − ρb − ∇ · T dV = 0 (1.6.5)
∂t
VA

De donde se extrae la primera ecuación de Cauchy

Primera ecuación de Cauchy

Dv
ρ = ρb + ∇ · T (1.6.6)
Dt
Para obtener este resultado se utilizó la forma microscópica de la ecuación de continuidad (ecuación
1.5.7) así como la definición de la derivada material (ecuación 1.3.2). El tensor de esfuerzos
totales puede descomponerse en los esfuerzos perpendiculares (debidos a la presión) y tangenciales
(debidos a los esfuerzos viscosos) al fluido como se muestra a continuación

T = −pI + τ (1.6.7)

donde el signo de menos en el primer término se debe a que la presión se aplica hacia el fluido y el
sentido del esfuerzo es hacia afuera del fluido. El tensor de esfuerzos cortantes en elsegundo término
de la ecuación anterior se puede relacionar con el tensor de deformación D = 12 ∇v + (∇v)T de


manera lineal de acuerdo a la ley de Newton de la viscosidad (Batchelor, 2000)

 
2
τ = 2µD + I κ − µ (∇ · v) (1.6.8)
3
donde µ es el coeficiente de viscosidad dinámica, κ es la viscosidad de dilatación (o viscosidad
de volumétrica). Los coeficientes µ y κ pueden cambiar con la presión y la temperatura y para
el caso de un flujo Newtoniano se consideran independientes de D. Se sabe que κ = 0 para gases
monoatómicos y es generalmente aceptado que κ  µ para otros fluidos. Una excepción de esta
suposición es un líquido que contenga burbujas, como un refresco. Para casos en los que el flujo se
incompresible, la ecuación de continuidad adquiere la forma dada en la ecuación (1.5.8) y el tensor
de esfuerzos cortantes se reduce a

τ = µ ∇v + (∇v)T
 
(1.6.9)

Sustituyendo esta expresión en la primera ecuación de Cauchy (ecuación 1.6.6) da lugar a la

conocida ecuación de Navier (1822)-Stokes (1880) (ver ejercicio 1.5):

Ecuación de Navier-Stokes

∂v
ρ + ρv · ∇v =− ∇p + ρg + µ∇2 v (1.6.10)
∂t | {z } |{z} |{z} | {z }
esfuerzos superficiales
|{z} esfuerzos inerciales esfuerzos esfuerzos superficiales
acumulación de normales volumétricos tangenciales
cantidad de movimiento
 20 Capítulo 1. Introducción: Repaso de mecánica de fluidos

donde el vector de fuerza volumétrica fue reemplazado por el vector de gravedad. Note que se
incluyó el significado físico de cada término. Esta ecuación, junto con la ecuación de continuidad
dada en la ecuación (1.5.8) forman un sistema cerrado de ecuaciones diferenciales, cuyas incógnitas
son la presión y la velocidad y con las cuales se puede estudiar una amplia variedad de problemas
como se muestra en los textos de mecánica de fluidos (Whitaker, 1992a; Batchelor, 2000).
Para concluir, esta sección es interesante notar que al dividir la ecuación (1.6.10) entre la
densidad y al reacomodar los términos resultantes se obtiene la siguiente expresión
∂v 2
∂t | ·{z
+ v ∇v} = ν∇
| {z v} + |{z}
f (1.6.11)
|{z} convección difusión fuente
acumulación

donde ν = µ/ρ es la viscosidad cinemática. La ecuación anterior está escrita en la forma clásica que
incorpora términos de acumulación, convección y difusión así como un término fuente f = −∇w + g,
donde w = p/ρ es el trabajo termodinámico específico (o por unidad de masa). Note que −∇w es
una fuente interna, mientras que g es una fuente externa. Esta forma matemática de la ecuación de
conservación de cantidad de movimiento es consistente con la que se usa en transferencia de calor
y en transferencia de masa (Bird y col., 2006).
Ejercicio 1.4 Utilice la ecuación de continuidad y la definición de derivada material para
deducir la ecuación (1.6.6) a partir de la ecuación (1.6.5). 

Ejercicio 1.5 Tome la divergencia de ambos lados de la ecuación (1.6.9) y tome en cuenta la
condición de flujo incompresible, para demostrar que ∇ · τ = µ∇2 v para casos en los que el
coeficiente de viscosidad sea constante. 

1.7 Ecuación de energía mecánica

Para los desarrollos que siguen es necesario deducir la ecuación de energía mecánica en una
fase. Esta ecuación no parte en sí de un principio fundamental, sino que se deduce a partir de
las ecuaciones de movimiento. Para ilustrar lo anterior se toma el producto interno de la primera
ecuación de Cauchy con el vector de velocidad del fluido v,
Dv
ρ · v = ρb · v + (∇ · T) · v (1.7.1)
Dt
Para avanzar, se toman en cuenta las siguientes identidades

Dv 1 Dv2
·v = (1.7.2a)
Dt 2 Dt
(∇ · T) · v = ∇ · (T · v) − T : ∇v (1.7.2b)

De esta forma, la ecuación (1.7.1) toma la siguiente forma

ρ Dv2
= ρb · v + ∇ · (T · v) − T : ∇v (1.7.3)
2 Dt
Esta forma en particular de la ecuación de energía mecánica se usará en el siguiente capítulo para
la deducción de la ecuación de energía térmica. Si se sustituye ahora la descomposición del tensor
de esfuerzos totales dada en la ecuación (1.6.7), el resultado es

ρ Dv2
= −∇ · (pv) + ρb · v + p(∇ · v) + ∇ · (ττ · v) − τ : ∇v (1.7.4)
2 Dt
1.7 Ecuación de energía mecánica 21

Más aún, para un flujo Newtoniano e incompresible se tiene que

τ = µ ∇v + (∇v)T
 
(1.7.5)

Por lo tanto, el penúltimo término del lado derecho de la ecuación (1.7.4) toma la siguiente forma

∇ · (ττ · v) = µ∇ · (∇v · v) + µ∇ · (∇v)T · v

(1.7.6)

O bien, tomando en cuenta la identidad

1
∇v · v = ∇v2 (1.7.7)
2
se tiene entonces que
µ 2 2
∇ v + µ∇ · (∇v)T · v


∇ · (ττ · v) = (1.7.8)
2
El último término del resultado anterior puede modificarse al notar que

∇ · (∇v)T · v = ∇ · (v · ∇v) = ∇v : ∇v + v · ∇ (∇ · v) = ∇v : ∇v


(1.7.9)

En la última igualdad se tomó en cuenta la suposición de flujo incompresible. De esta forma, la

ecuación (1.7.4) toma la siguiente forma:

ρ Dv2 1
= −∇ · (pv) + ρb · v + µ∇2 v2 + µ∇v : ∇v − τ : ∇v (1.7.10)
2 Dt 2
Dirigiendo la atención al último término de la ecuación anterior, el resultado de sustituir la ley de
Newton de la viscosidad es:

τ : ∇v = µ∇v : ∇v + µ(∇v)T : ∇v (1.7.11)

Por lo que la forma final de la ecuación de energía mecánica es:

Ecuación de energía mecánica

ρ ∂ v2
 
1
+ v · ∇v = −∇ · (pv) + ρb · v + µ∇2 v2 − µ∇vT : ∇v
2
(1.7.12)
2 ∂t 2

Los términos del lado izquierdo de la ecuación anterior representan la acumulación y el transporte
convectivo de energía cinética, respectivamente. Por su parte, los términos del lado derecho
representan: −∇ · (pv) : la tasa de trabajo de presión ejercida por los alrededores hacia el fluido;
ρb · v la tasa de trabajo debida a fuerzas volumétricas; 12 µ∇2 v2 transporte difusivo de energía
cinética; −µ∇vT : ∇v la tasa de conversión irreversible de energía mecánica a energía interna o
disipación viscosa.
Ejercicio 1.6 Utilice notación indicial para demostrar las identidades dadas en las ecuaciones
(1.7.2), (1.7.7) y (1.7.9). 

Ejercicio 1.7 Deduzca el resultado dado en la ecuación (1.7.12) a partir de tomar el producto
punto de la ecuación de Navier-Stokes con el vector de velocidad. 
 22 Capítulo 1. Introducción: Repaso de mecánica de fluidos

Ejercicio 1.8 Exprese la ecuación (1.7.12) en una forma similar a la dad en la ecuación (1.6.11)
y discuta el significado físico de las fuentes y sumideros de energía resultantes. 

1.8 Resumen
En este capítulo, se revisó la definición y clasificación de un modelo y se ubicó la atención en
los modelos matemáticos deterministas basados en principios fundamentales. Tras dar un panorama
general de la transferencia de calor y revisar el herramental matemático (en concreto, los tres
tipos de derivadas temporales y los teoremas integrales), se dedujeron las formas macroscópica y
microscópica de la ecuación de conservación de masa:

d
Z Z
ρ dV + n · ρ(v − w) dA = 0 (1.8.1a)
dt
VA AA
∂ρ
+ ∇ · (ρv) = 0 (1.8.1b)
∂t
La deducción consistió en partir de un axioma (o principio de conservación en palabras) y, a partir
del uso de los teoremas del transporte y de la divergencia, se extrajo la versión microscópica al
conseguir formular la versión macroscópica en la forma de una integral volumétrica igualada a
cero.
Más adelante, se presentó de forma axiomática el principio de conservación de cantidad de
movimiento lineal, el cual en lenguaje matemático se expresa como sigue

d
Z Z Z Z
ρv dV + n · ρv(v − w) dA = ρb dV + n · T dA (1.8.2a)
dt
VA AA VA AA

note que se tomó en cuenta el teorema fundamental de Cauchy en el último término. Siguiendo
el mismo procedimiento utilizado en la deducción de la ecuación de conservación de masa, se
recuperó la primera ecuación de Cauchy, la cual es la versión microscópica de la ecuación anterior.
Ya que esta ecuación no está cerrada, se utilizó una descomposición del tensor total de esfuerzos
en sus partes normal y tangencial como se muestra en la ecuación (1.6.7), junto con la ley de
Newton de la viscosidad para flujo incompresible (ecuación 1.6.9) para así obtener la ecuación de
Navier-Stokes
∂v
ρ + ρv · ∇v = −∇p + ρg + µ∇2 v (1.8.2b)
∂t
Por último, se dedujeron diversas formas de la ecuación de energía mecánica que resultan de tomar
el producto punto de la primera ecuación de Cauchy (o bien de la ecuación de Navier Stokes) con el
vector de velocidad. En específico, la forma de la ecuación de energía mecánica que será de interés
para las deducciones que siguen es la dada en la ecuación (1.7.3).
 2. Deducción de la ecuación de energía

2.1 Introducción

En el capítulo anterior se revisaron los fundamentos de mecánica de fluidos, así como el

herramental matemático necesario para su deducción. Con estos antecedentes, puede centrarse la
atención ahora en la deducción de la ecuación de energía térmica, la cual es el objeto de estudio del
resto del texto. Para ello, se partirá de la forma macroscópica e integrada en el tiempo de la primera
ley de la termodinámica tal y como se presenta en los cursos introductorios de termodinámica.
A partir de dicha ecuación se deducirá la forma local e instantánea de la ecuación de energía
total. Al restarle la ecuación de energía mecánica a dicha ecuación, el resultado es la ecuación de
energía térmica. El modelo resultante no está cerrado, pues tiene más incógnitas que ecuaciones. En
específico, es necesario contar con una relación entre la densidad de flujo de calor y la temperatura.
Para resolver esta situación se utilizará la ley de Fourier y se discuten además sus variantes.
Por último, se deduce la condición de salto para la densidad del flujo de calor entre dos fases
y se postulará la ley de enfriamiento de Newton. Con el material cubierto en este capítulo, se
contarán con los elementos necesarios para abordar una gran cantidad de problemas asociados a la
transferencia de calor. Los cuales se discutirán en la Parte II de este texto.

2.2 Primera ley de la termodinámica

El principio fundamental de conservación de energía es la primera ley de la termodinámica y

puede enunciarse como sigue

Primera ley de la termodinámica (expresada en palabras)

( ) ( )
El cambio de energía interna Suma del calor y trabajo que actúan
=
cinética y potencial de un sistema sobre la superficie del sistema
(2.2.1)
24 Capítulo 2. Deducción de la ecuación de energía

En lenguaje matemático, la ecuación anterior toma la siguiente forma algebraica:

∆U + ∆Ec + ∆E p = Q +W (2.2.2)

Donde los cambios de energía interna, cinética y potencial se refieren a la diferencia entre los
estados final e inicial de dichas propiedades. Estos estados inicial y final, pueden asociarse a
instantes de tiempo específicos, los cuales se denotan como t1 y t2 , respectivamente. De esta forma,
se tiene que
t=t
Z 2
∆ψ = ψ(t2 ) − ψ(t1 ) = dψ; ψ = U, Ec , E p (2.2.3a)
t=t1

Por su parte Q y W pueden expresarse también como integrales temporales de las tasas a las cuales
se suministra calor y trabajo al sistema en cada instante, esto es,
t=t
Z 2
Q= Q̇dt (2.2.3b)
t=t1
t=t
Z 2
W= Ẇ dt (2.2.3c)
t=t1

Por lo tanto, la primera ley de la termodinámica también puede escribirse como una sola integral
temporal:
t=t
Z 2 

dU + dEc + dE p − Q̇ + Ẇ dt = 0 (2.2.4)
t=t1

De donde se puede extraer la siguiente ecuación diferencial ordinaria:

DU DEc DE p
+ + = Q̇ + Ẇ (2.2.5)
Dt Dt Dt
Esta expresión es la forma instantánea de la primera ley de la termodinámica. Debe quedar claro
que todas las cantidades involucradas en la ecuación anterior son macroscópicas pues se refieren al
sistema, que en este contexto se refiere a un cuerpo o bien a una región material (ver definición
1.3.1). Por ello, las derivadas que aparecen en la ecuación anterior son derivada materiales.
Ya que U, Ec y E p son cantidades macroscópicas, pueden expresarse en términos de sus
contrapartes microscópicas (y por unidad de masa) mediante las siguientes integrales volumétricas
Z
U= ρe dV (2.2.6a)
VM
v2
Z
Ec = ρ dV (2.2.6b)
2
VM
Z
Ep = ρφ dV (2.2.6c)
VM

Donde e y φ denotan la energía interna por unidad de masa y la energía potencial por unidad de
masa, respectivamente.
2.2 Primera ley de la termodinámica 25

Dirigiendo la atención al primer término en el lado derecho de la ecuación (2.2.5), se sabe que
la energía suministrada en forma de calor o por medio de radiación es recibida por la superficie del
cuerpo. Por lo que si se denota por q y qR a la densidad de flujo (o flux) de calor (es decir, el calor
por unidad de área y de tiempo) y de radiación hacia el cuerpo, se tiene entonces que
( )
La tasa a la cual se suministra energía Z
Q̇ = = − n · (q + qR ) dA (2.2.7)
al cuerpo por calor y radiación
AM

donde el signo negativo se debe a que el vector unitario, n, está dirigido hacia afuera del cuerpo,
mientras que los flujos de calor van dirigidos hacia dentro del cuerpo.
Por último, la tasa a la cual se suministra trabajo (es decir, la potencia) en las superficies del
sistema puede expresarse en términos del producto interno del vector de esfuerzos y la velocidad
del fluido, o bien usando el teorema fundamental de Cauchy (ecuación 1.6.3) como sigue:
Z Z
Ẇ = tn · v dA = n · T · v dA (2.2.8)
AM AM

Usando estas definiciones en la ecuación (2.2.5), la forma instantánea de la primera ley de la

termodinámica es ahora:
Z  
D 1
Z
ρe + ρv2 + ρφ dV = n · (−q − qR + T · v) dA (2.2.9)
Dt 2
VM AM

Dirigiendo la atención al último término del lado izquierdo de la ecuación anterior se tiene, con la
ayuda del teorema modificado de Reynolds (ver ecuación 1.5.10), que:
Z  
D Dφ
Z Z
∂φ
ρφ dV = ρ dV = ρ + v · ∇φ dV (2.2.10)
Dt Dt ∂t
VM VM VM

Dado que la energía potencial no cambia con el tiempo, la ecuación anterior se reduce a:
D
Z Z Z
ρφ dV = ρv · ∇φ dV = − ρb · v dV (2.2.11)
Dt
VM VM VM

En la última igualdad se utilizó la definición del vector de fuerzas volumétricas por unidad de masa:

b = −∇φ (2.2.12)
De esta forma, la ecuación (2.2.9) puede escribirse como:
Z  
D 1 2
Z Z
ρe + ρv dV = n · (−q − qR + T · v) dA + ρb · v dV (2.2.13)
Dt 2
VM AM VM

Este resultado es consistente con la siguiente ecuación con palabras que puede encontrarse en textos
de fenómenos de transporte (Bird y col., 2006; Whitaker, 1983)

( ) ( )
La razón de cambio en el tiempo de la La tasa a la cual se suministra energía
=
energía interna y cinética de un cuerpo al cuerpo por calor y trabajo
( )
La tasa a la cual se suministra energía
+
al cuerpo por radiación
(2.2.14)
 26 Capítulo 2. Deducción de la ecuación de energía

Note que en la expresión anterior se denota por calor a la energía térmica transmitida por contacto
con una superficie, mientras que la radiación se refiere a la energía térmica transferida por fotones
que llegan al cuerpo desde una cierta distancia.
Como se mostró en el capítulo anterior, las fuerzas ejercidas sobre un cuerpo pueden ser
volumétricas y superficiales, por lo que la tasa a la cual se ejerce trabajo sobre el cuerpo puede
expresarse a su vez como
( ) Z
La tasa a la cual se suministra energía Z
= n · T · v dA + ρb · v dV (2.2.15)
al cuerpo por trabajo
AM V
| {z } |M {z }
trabajo superficial trabajo volumétrico

Como es evidente, la ecuación (2.2.13) es válida para una región material. Gray y Gray (2017)
proporcionan una versión más general de esta ecuación, la cual es aplicable a una región arbitraria
Ecuación de energía total (forma macroscópica)
Z   Z  
d 1 1 2
Z
ρe + ρv2 dV + ρe + ρv (v−w)·n dV = n·(−q − qR + T · v) dA
dt 2 2
VA AA AA
Z
+ ρb · v dV (2.2.16)
VA
Como puede notarse, para el caso de una región material w = v y se recupera la ecuación
(2.2.13). A pesar de que la ecuación de arriba es más general, los desarrollos que siguen partirán de
la ecuación (2.2.13) con el ánimo de continuar el análisis basado en una región material, aunque el
resultado final es independiente del tipo de región de donde se comience el análisis.

Ejercicio 2.1 Exprese en forma de ecuación con palabras la ecuación (2.2.16). 

Ejercicio 2.2 Partiendo de la ecuación (2.2.16), deduzca la expresión correspondiente para una
región fija. Exprese el resultado tanto en forma matemática como en forma de ecuación con
palabras. 

Ejercicio 2.3 Utilice la descomposición del tensor de esfuerzos totales dada en la ecuación
(1.6.7) para expresar a la ecuación (2.2.13) en términos de la presión y el tensor de esfuerzos
viscosos. 

2.3 Ecuación de energía total

Para deducir la forma local e instantánea de la ecuación de conservación de energía, es con-
veniente expresar la ecuación (2.2.13) como una sola integral de volumen igualada a cero tal y
como se hizo en el caso de la ecuación de continuidad (sección 1.5) y de cantidad de movimiento
(sección 1.6). Como primer paso, se utiliza el teorema modificado de Reynolds (ecuación 1.5.10)
en el término del lado izquierdo de la ecuación (2.2.13), lo que da como resultado
Z     Z     
D 1 2 D 1 2 1 2 1 2
Z
∂
ρe + ρv dV = ρ e + v dV = ρ e+ v +v·∇ e+ v dV
Dt 2 Dt 2 ∂t 2 2
VM VM VM
(2.3.1a)
donde en la última igualdad se utilizó la definición de la derivada material (ecuación 1.3.2). Con
este paso ya no se tiene una derivada de una integral, sino una integral volumétrica de términos con
 2.4 Ecuación de energía térmica 27

derivadas. El siguiente paso es entonces aplicar el teorema de la divergencia en el primer término

del lado derecho de la ecuación (2.2.13), para transformar dicha integral de superficie a la siguiente
integral volumétrica
Z Z
n · (−q − qR + T · v) dA = ∇ · (−q − qR + T · v) dV (2.3.1b)
AM VM

El resultado de sustituir las ecuaciones (2.3.1) en la ecuación (2.2.13) es la siguiente expresión,

Z    
D 1 2
ρ e + v + ∇ · (q + qR − T · v) − ρb · v dV = 0 (2.3.2)
Dt 2
VM

De donde se extrae la ecuación de energía total (Bird y col., 2006)

Ecuación de energía total (versión microscópica)

 
D 1 2
ρ e + v = −∇ · (q + qR ) + ∇ · (T · v) + ρb · v (2.3.3)
Dt 2

A pesar de la estructura que posee la ecuación anterior, su significado físico permanece intacto
respecto a su contraparte macroscópica (ecuación 2.2.13). De hecho, el lado izquierdo de la ecuación
anterior representa la razón de cambio de la energía interna y cinética de un cuerpo (lo cual incluye
tanto la acumulación como el transporte por convección), mientras que el primer término del
lado derecho representa la tasa a la cual se suministra calor (por conducción y radiación) y los
últimos dos términos se refieren al trabajo (superficial y volumétrico) realizado sobre el sistema.
En otras palabras, la ecuación (2.3.3) representa la forma local e instantánea de la primera ley de
la termodinámica. En los siguientes párrafos se muestra cómo se pude obtener, a partir de esta
ecuación, la ecuación de energía térmica.

2.4 Ecuación de energía térmica

La ecuación de energía total (ecuación 2.3.3) combina tanto las contribuciones térmicas como
mecánicas de transporte de energía que actúan en cada punto de una fase. En la sección 1.7, se
presentó la deducción de la ecuación de energía mecánica, la cual, en términos del tensor de
esfuerzos totales es (ver ecuación 1.7.3):
ρ Dv2
= ρb · v + ∇ · (T · v) − T : ∇v (2.4.1)
2 Dt
Para aislar las contribuciones exclusivas por transporte de energía térmica de la ecuación de energía
total, se substrae la ecuación (2.4.1) a la ecuación (2.3.3) para obtener la siguiente expresión
De
ρ = −∇ · (q + qR ) + T : ∇v (2.4.2)
Dt
la cual es una primera forma microscópica de la ecuación de energía térmica. A diferencia de
la ecuación de energía mecánica, la cual requiere conocimiento de la velocidad y el tensor de
esfuerzos totales, los cuales provienen de resolver las ecuaciones de continuidad y de cantidad de
movimiento; la ecuación anterior tiene como incógnitas a la energía interna por unidad de masa,
así como las densidades de flujo por conducción y radiación. Estas variables no están disponibles
en este punto del análisis, por lo que el modelo no está cerrado. Para atender esta situación, se
comienza por sustituir la descomposición del tensor de esfuerzos totales dada en la ecuación (1.6.7)
en el último término de la ecuación anterior para obtener
De
ρ = −∇ · q − ∇ · qR − p∇ · v + τ : ∇v (2.4.3)
Dt
28 Capítulo 2. Deducción de la ecuación de energía

Con el fin de expresar al lado izquierdo de la ecuación (2.4.3) en función de la temperatura del
cuerpo, se recurre a la definición de entalpía por unidad de masa (Callen, 1985):
p
h = e+ (2.4.4)
ρ
Aplicando el operador de derivada material en ambos lados de la ecuación anterior da como
resultado
Dh De 1 Dp p Dρ
= + − (2.4.5)
Dt Dt ρ Dt ρ 2 Dt
De
Despejando de la ecuación anterior el término Dt y, tomando en cuenta la forma de la ecuación de
continuidad dada por (ver ecuación 1.5.9)
Dρ
= −ρ∇ · v (2.4.6)
Dt
se obtiene el siguiente resultado
De Dh 1 Dp p
= − − ∇·v (2.4.7)
Dt Dt ρ Dt ρ
Sustituyendo la expresión anterior en el término del lado izquierdo de la ecuación (2.4.3), da lugar
a
Dh Dp
ρ = −∇ · q − ∇ · qR + τ : ∇v + (2.4.8)
Dt Dt
Lo que se ha conseguido es expresar a la energía interna por unidad de masa en términos de la
entalpía por unidad de masa y la presión. Sin embargo, sin una expresión para la entalpía, no hay
un progreso real en el análisis. Para atender esta situación, es conveniente recordar que la entalpía
es función de la temperatura y la presión (ver, por ejemplo Callen, 1985) , y dicha dependencia
puede expresarse como sigue

Dh Dh DT Dh Dp DT 1 Dp
= + = Cp + (1 − αT ) (2.4.9)
Dt DT p Dt Dp T Dt Dt ρ Dt

donde Cp y α son la capacidad calorífica a presión constante y el coeficiente de expansión térmica,

respectivamente (ver ejercicio 2.4). Este último coeficiente se define como

1 Dρ
α =− (2.4.10)
ρ DT p

Sustituyendo la ecuación (2.4.9) en la ecuación (2.4.8) y usando la definición de la derivada

material (ecuación 1.3.2) permite llegar a la siguiente expresión
 
∂T
ρCp + v · ∇T = −∇ · q + Φ (2.4.11)
∂t
Para hacer más breve esta expresión se utilizó la siguiente definición
Dp
Φ = τ : ∇v + T α − ∇ · qR (2.4.12)
Dt
Como en el caso de la conservación de masa para una especie química (ver, por ejemplo el texto
de Bird y col., 2006), note que si Φ > 0, existe generación de calor en el sistema y si Φ < 0 el
 2.5 Ley de Fourier 29

sistema tiene pérdidas de calor. Aunque se conservará esta forma compacta, no debe olvidarse que
Φ contiene contribuciones debidas a la disipación viscosa, los cambios materiales de la presión y la
radiación.
La ecuación (2.4.11) significa un considerable progreso hacia un modelo cerrado en el que
la incógnita sea la temperatura, ya que la presión y la velocidad se obtienen de las ecuación de
conservación de masa y cantidad de movimiento. Sin embargo, aún permanecen dos incógnitas
más que deben determinarse, estas son la densidad de flujo de calor por conducción y radiación (q
y qR , respectivamente). La segunda puede determinarse a partir del conocimiento de la potencia
de la fuente radiativa y se aborda en otra parte del texto. En los siguientes párrafos se presenta la
relación de cerradura entre el flux conductivo y el gradiente de la temperatura.

Ejercicio 2.4 — Entalpía y expansión térmica. Partiendo de la siguiente relación entre la

entalpía por unidad de masa, la entropía por unidad de masa (s) y la presión para sistemas
cerrados (es decir, regiones materiales),
Dp
Dh = T Ds + (2.4.13)
ρ

y tomando en cuenta que

Ds Dρ −1
=− (2.4.14)
DP T DT P

demuestre que

Dh 1
= (1 − αT ) (2.4.15)
Dp T ρ


Ejercicio 2.5 — Disipación viscosa. En el caso en el que exista flujo de un fluido, ¿qué signo
debe tener τ : ∇v? 

Ejercicio 2.6 — Forma macroscópica de la ecuación de energía térmica. Integre la

ecuación (2.4.11) en una región arbitraria y explique, bajo qué condiciones, puede obtenerse la
siguiente forma macroscópica de la ecuación de energía térmica

d
Z Z Z
ρCp T dV + n · [ρCp T (v − w) + q] dA = Φ dV (2.4.16)
dt
VA AA VA

Explique el significado físico de cada término para así poder expresar este resultado en forma
de ecuación con palabras. 

2.5 Ley de Fourier

Como elemento final para cerrar la ecuación de energía térmica, es necesario determinar la
relación que existe entre el flux conductivo de calor q y la temperatura. Con este fin, considere el
sistema experimental mostrado en la figura 2.1a) (Whitaker, 1983) en donde se colocan dos placas
paralelas de acero entre una cama de aire. Inicialmente el sistema estaba a 293 K, posteriormente la
lámina inferior se mantuvo a esa temperatura, mientras que la placa superior se somete a una fuente
de calentamiento que va de Q̇/A = 100 a 1500 W/m2 . Se midió la temperatura en la superficie de
30 Capítulo 2. Deducción de la ecuación de energía

material a T1

aislante

aire

material a T0

a)

40

Q̇ 30

A
20

10 1

400 600 800 1,000 1,200 1,400

∆T /L
b)
Figura 2.1: a) Esquema del dispositivo experimental para determinar la conductividad térmica; b)
Resultados experimentales del flux de calor vs ∆T /L.

la lámina superior para cada valor de Q̇/A manteniendo constate la separación de las placas en
L = 0.1 m. Otra opción es mantener el flujo de calor constante y cambiar la separación entre las
placas. Los resultados se muestran en la figura 2.1b), donde ∆T = T1 − T0 . Como puede notarse,
los datos obedecen una tendencia lineal, esto es:
Q̇ ∆T
q= =k (2.5.1)
A L
donde k es la pendiente de la recta y se conoce como coeficiente de conductividad térmica del aire.
Para este experimento, el valor encontrado fue de 0.032333968 W/(m K). Dicho valor se encuentra
en el rango esperado de la conductividad térmica del aire.
La conductividad térmica depende débilmente de la temperatura y muy poco de la presión para
líquidos y gases. La conductividad térmica de los sólidos es, en general mayor que para los líquidos
y la de éstos es a su vez mayor que la de los gases. El matemático francés, Jean Baptiste Fourier
(1878) fue el primero en notar esta relación lineal entre la densidad de flujo de calor y el cambio de
2.5 Ley de Fourier 31

temperatura a finales del siglo XIX (al lector interesado en la historia del trabajo de Fourier se le
recomienda el texto de Narasimhan, 1999). Tomando el límite cuando la separación de las placas
tiende a cero, lleva a la siguiente forma diferencial de la ley de Fourier:

∂T
qz = −k (2.5.2)
∂z
El signo negativo se agrega para que k siempre sea positiva ya que ∂ T /∂ z es positivo pero el flux de
calor es negativo. Lo anterior se debe a que, de acuerdo a la termodinámica, la transferencia de calor
siempre se da desde la región de mayor temperatura hacia la zona de menor temperatura. Como
puede notarse, la ley de Fourier no fue deducida formalmente, sino que vino de la observación
experimental. Esta es una muestra de la ayuda que pueden aportar los modelos matemáticos basados
en la observación en los modelos matemáticos basados en principios fundamentales (ver sección
1.1).
En forma vectorial, la ley de Fourier puede escribirse como:

Ley de Fourier (forma vectorial)

q = −k∇T (2.5.3)

Para el experimento mencionado arriba, el valor de la conductividad térmica es constante, sin

embargo, puede darse el caso en el que este coeficiente sea función de la temperatura o de la
posición. Para materiales anisótropos, la conductividad térmica es un tensor de segundo orden y la
ley de Fourier se expresa en una forma más general como:

q = −K · ∇T (2.5.4)

Los valores de la conductividad térmica de los materiales se obtienen de experimentos, correlaciones

o bien a partir de teoría.
Utilizando la forma de la ley de Fourier dada en la ecuación (2.5.3), la ecuación de energía
térmica se expresa como:

Ecuación de energía térmica

 
∂T
ρCp + v · ∇T = ∇ · (k∇T ) + Φ (2.5.5)
∂t

de acuerdo a este resultado, la transferencia de calor se refiere a la energía térmica transferida debido
a cambios espacio temporales de la temperatura (Bergman y col., 2017). El primer término del lado
izquierdo de esta ecuación representa la acumulación de calor, el segundo término representa el
transporte por convección de calor, el primer término del lado derecho representa la conducción
de calor y el último término se refiere a la fuente o sumidero de calor en la fase. Siempre que
este último término sea conocido, la ecuación anterior está cerrada pues la única incógnita es la
temperatura. En los siguientes párrafos se discute acerca de los tres mecanismos de transferencia
de calor incluidos en la ecuación anterior, para finalmente discutir acerca de las condiciones de
frontera correspondientes.

Ejercicio 2.7 — Difusividad térmica. El coeficiente de difusividad térmica se define como

sigue

k
αT = (2.5.6)
ρCp
 32 Capítulo 2. Deducción de la ecuación de energía

Con base en las unidades de k, ρ y Cp deduzca las unidades de αT y exprese a la ecuación de

energía térmica en función de la difusividad térmica.


Ejercicio 2.8 — Ecuación de Cattaneo-Vernotte. La ley de Fourier supone que el gradiente

de temperatura se propaga a una velocidad infinita, de tal manera que instantáneamente se produ-
ce un flux conductivo de calor. Para atender esta situación Maxwell (1867), y más recientemente
Cattaneo (1958) y Vernotte (1958) propusieron la siguiente modificación a la ley de Fourier

∂q
tr + q = −k∇T (2.5.7)
∂t
donde tr es una constante de tiempo que se conoce como el tiempo de relajación. Evidentemente,
cuando esta constante es cero, se recupera la ley de Fourier. Demuestre que si se usa la ecuación
de Cattaneo-Vernotte como ecuación de cerradura, entonces la forma cerrada de la ecuación de
energía térmica puede escribirse como sigue
 2
∂ T ∂ (v · ∇T )
  
∂T ∂Φ
ρCptr + + ρC p + v · ∇T = ∇ · (k∇T ) + Φ + tr (2.5.8)
∂t 2 ∂t ∂t ∂t


2.6 Mecanismos de transferencia de calor

A pesar de que son muchas las formas en las que puede transferirse el calor de una región a otra,
desde el punto de vista del continuo, todos los mecanismos pueden entrar en solo tres categorías
que son: conducción, convección y radiación. Estos mecanismos o modos de transferencia de calor
no sólo están presentes en fases (como se puede notar en la ecuación de energía térmica) sino
también en interfases. Aunque ya se han mencionado en párrafos anteriores estos mecanismos, es
conveniente dedicar un espacio a definirlos con precisión.

2.6.1 Conducción
Este modo de transferencia de calor consiste en la transferencia de energía de partículas con
más energía hacia aquellas con menor energía; es decir, se basa en el contacto directo entre las
partículas. Por lo tanto, la conducción de calor es imposible en el vacío y puede tener lugar en
cualquier fase o interfase. De esta forma puede afirmarse que siempre que exista un gradiente de
temperatura, es decir cuando los cambios espaciales de la temperatura sigan una trayectoria definida
de una región más energética a una menos energética, tendrá lugar la transferencia de calor por
conducción. En otras palabras, siempre que se obedezca la ley de Fourier se pude afirmar que se da
este modo de transferencia de calor. Dado que este modo de transferencia de calor se da por los
choques entre partículas y dichos choques son aleatorios, puede hablarse de un transporte difusivo
de calor .
 Ejemplo 2.1 — Transferencia de calor en una hielera. Una hielera tiene las siguientes
dimensiones internas: ancho 51 cm, profundidad: 28.6 cm, alto: 32 cm y externas: ancho: 63 cm,
profundidad: 36.8 cm, alto: 38.1 cm. Si la temperatura en la superficie externa de la hielera es
constante y vale 30 o C y en el interior se tienen paletas que están a -10o C, calcule la tasa de
transferencia de calor en estado estacionario del exterior hacia la hielera en Watts. La conductividad
térmica del plástico es de 0.12 W/(m K).
Solución:
En este problema, el interés está en conocer la transferencia de calor que recibe el prisma
 2.6 Mecanismos de transferencia de calor 33

interno de la hielera. Para ello, se procede de la siguiente manera, primero se utilizará la forma
integrada de la ley de Fourier para calcular el flux de calor sobre cada superficie del prisma y
posteriormente cada flux de multiplicará por el área de la sección transversal a través de la cual se
transporta el calor hacia cada superficie. Por último se suman las contribuciones de cada superficie.
Por conveniencia, se asigna el eje x como aquel que corre a lo largo de la profundidad del
compartimento interno, mientras que el eje y atraviesa la anchura y el eje z la altura. Ya que el
problema especifica que las temperaturas son constantes tanto en el exterior como en el interior de
la hielera, se tiene entonces que
−10 − 30 W
qx = −2k = 80 2
0.12 m
−10 − 30 W
qy = −2k = 117.1 2
0.082 m
−10 − 30 W
qz = −2k = 157.4 2
0.061 m
Multiplicando por la áreas correspondientes de las superficies de la caja interna se obtiene como
resultado la tasa de transferencia de calor que recibe dicha caja
Q̇ = (80)(0.51)(0.32) + (117.1)(0.286)(0.32) + (157.4)(0.51)(0.286) = 46.7W


2.6.2 Convección
Este modo de transporte de energía es el resultado del movimiento de un fluido. Aunque algunos
autores (Gray y Gray, 2017; Bergman y col., 2017) prefieren llamar a este mecanismo advección
y utilizan al término de convección para referirse a la suma de la advección y la conducción (o
difusión) de calor. Como es de esperarse, la convección sólo tiene lugar en fluidos. La convección
puede tener lugar tanto en una fase como en una interfase; para ser más precisos, en una capa límite
como se verá más adelante en el texto. De hecho, es en la capa límite donde cobra más sentido la
definición de convección mencionada arriba como un mecanismo acumulativo de la advección y la
conducción de calor.
De acuerdo a la naturaleza del flujo, la convección puede ser forzada cuando se origina por
mecanismos externos como un agitador o un ventilador, o bien convección natural o libre cuando
el movimiento se da por diferencias de densidades inducidas por diferencias de temperatura en el
fluido, como es el caso de un plato de sopa caliente o una taza de café. A pesar de esta distinción,
no es raro que ambos tipos de convección tengan lugar en diversas aplicaciones, por ejemplo en
las computadoras se tienen ventiladores para que no se sobrecaliente el procesador, lo cual es una
forma de convección forzada. Sin embargo, el aire dentro de la computadora puede experimentar
diferencias de densidad debido a la temperatura y ocasionar un movimiento global, dando lugar
también a convección natural.

2.6.3 Radiación
Es el transporte de energía que se da en ausencia de materia y consiste en el transporte a
través de ondas electromagnéticas. Esta característica distingue a la radiación de la conducción y
convección, los cuales requieren de materia para poder tener lugar. El transporte de energía por
radiación se da a la velocidad de la luz por medio de fotones y se da más eficientemente en el vacío.
Para fines ingenieriles, el interés está en la radiación térmica, la cual se encuentra en el rango de
longitudes de onda entre el infrarrojo (10−2 cm) y el ultravioleta (10−5 cm). La tasa de emisión de
energía radiante desde una superficie sólida está dada por la ley de Stefan-Boltzmann:
qR · n|emisión = qR |emisión = εσ T 4 (2.6.1)
 34 Capítulo 2. Deducción de la ecuación de energía

donde ε ∈ [0, 1] es la emisividad del cuerpo, σ = 5.67 × 10−8 W/m2 K4 es la constante de Steffan-
Boltzmann y T es la temperatura absoluta del cuerpo. Para un cuerpo negro ε = 1 y se conoce
como un radiador ideal, cualquier superficie real emite menos calor que un cuerpo negro y su valor
depende de la superficie del material, rugosidad, etc. Cabe recordar que qR es el flux de calor, es
decir, es la cantidad de calor por unidad de tiempo y de área que se transfiere. Si se desea conocer
la tasa de transferencia de calor de un cuerpo de superficie A1 se tiene entonces que:

Q̇ = A1 qR [=]W (2.6.2)

Hasta este punto, el interés ha estado en la emisión de un cuerpo, sin embargo, es de especial
interés también calcular la tasa de calor recibido por radiación sobre una superficie. En forma
matemática, el flux de calor absorbido sobre una superficie se calcula como

qR |absorción = ασ T 4 (2.6.3)

donde α ∈ [0, 1] es la fracción de energía incidente que absorbe la superficie del cuerpo emisor y se
conoce como absortancia. En un cuerpo negro ε = α = 1; siempre que α < 1 y la superficie sea
opaca, habrán porciones de la radiación incidente que serán reflejadas. De hecho, los valores de α
dependen tanto de la naturaleza de la fuente radiante como de la superficie receptora.
En la realidad todo cuerpo que emite energía también absorbe energía radiante, por lo que la
tasa neta a la cual un cuerpo de área A1 y temperatura superficial Ts transfiere energía térmica por
radiación desde su superficie hacia los alrededores que se encuentran a una temperatura Ta es

= A1 ( qR |emisión − qR |absorción ) = A1 σ εTs4 − αTa4

Q̇ netoemitido
(2.6.4)

Para un cuerpo gris, ε = α 6= 1 y entonces:

= εA1 σ Ts4 − Ta4

Q̇ netoemitido
(2.6.5)

2.7 Condiciones de frontera

2.7.1 Condición de salto del flux de calor
Como último tema de este capítulo, se deduce la condición de frontera para el flux de calor
entre dos fases. Para ello, considere un sistema de dos fases fluidas como el mostrado en la figura
2.2. Es conveniente comenzar el análisis con la forma macroscópica de la ecuación de energía
térmica dada en la ecuación (2.4.16), la cual es válida tanto en las fases como en la zona de cambios

d
Z Z Z
ρCp T dV + n · [ρCp T (v − w) + q] dA = Φ dV (2.7.1)
dt
VA AA VA

Si el centroide de esta región arbitraria de integración se ubica sobre la superficie divisoria entre las
fases, entonces es posible hacer las siguientes descomposiciones
Z Z Z Z
ρCp T dV = ρβ Cpβ Tβ dV + ργ Cpγ Tγ dV + ρsCps Ts dA (2.7.2a)
VA Vβ Vγ Aβ γ
Z Z Z Z
Φ dV = Φβ dV + Φγ dV + Φs dA (2.7.2b)
VA Vβ Vγ Aβ γ
2.7 Condiciones de frontera 35

Sistema aceite-agua

nγ
fase-γ

superficie
divisoria nβ
γ

fase-β

nβ

Figura 2.2: Esquema de un sistema de dos fases con detalle en su superficie divisoria

donde el subíndice s se utiliza para denotar propiedades superficiales. Ahora bien, dado que
AA = Aβ +Aγ , el segundo término en el lado izquierdo de la ecuación (2.7.1) puede descomponerse
como sigue

Z Z   Z  
n·[ρCp T (v − w) + q] dA = nβ · ρβ Cpβ Tβ (vβ − w) + qβ dA+ nγ · ργ Cpγ Tγ (vγ − w) + qγ dA
AA Aβ Aγ
(2.7.2c)

Sustituyendo los resultados dados en las ecuaciones (2.7.2) en la ecuación (2.7.1) da lugar a la
siguiente expresión

d d d
Z Z Z
ρβ Cpβ Tβ dV + ργ Cpγ Tγ dV + ρsCps Ts dA
dt dt dt
Vβ Vγ Aβ γ
Z   Z  
+ nβ · ρβ Cpβ Tβ (vβ − w) + qβ dA + nγ · ργ Cpγ Tγ (vγ − w) + qγ dA
Aβ Aγ
Z Z Z
= Φβ dV + Φγ dV + Φs dA (2.7.3)
Vβ Vγ Aβ γ
36 Capítulo 2. Deducción de la ecuación de energía

Por su parte, las versiones macroscópicas de la ecuación de energía térmica en cada fase pueden
expresarse como sigue

d
Z Z  
ρβ Cpβ Tβ dV + nβ · ρβ Cpβ Tβ (vβ − w) + qβ dA
dt
Vβ Aβ
Z   Z
+ nβ γ · ρβ Cpβ Tβ (vβ − w) + qβ dA = Φβ dV (2.7.4a)
Aβ γ Vβ

d
Z Z  
ργ Cpγ Tγ dV + nγ · ργ Cpγ Tγ (vγ − w) + qγ dA
dt
Vγ Aγ
Z   Z
+ nγβ · ργ Cpγ Tγ (vγ − w) + qγ dA = Φγ dV (2.7.4b)
Aβ γ Vγ

Restando estas ecuaciones a la ecuación (2.7.3) da como resultado

d
Z Z  

ρsCps Ts dA = nβ γ · ρβ Cpβ Tβ (vβ − w) − ργ Cpγ Tγ (vγ − w) + qβ − qγ + Φs dA
dt
Aβ γ Aβ γ
(2.7.5)

Como puede notarse, casi toda la ecuación puede agruparse dentro de una integral de superficie,
excepto por el término de acumulación en el lado izquierdo. Para hacer frente a esta situación, se re-
curre al teorema del transporte superficial (ver Sección 3 en el trabajo de Whitaker, 1992b):

Teorema 2.7.1 — Teorema del transporte superficial. Sea ψs una función continua y definida
en la superficie divisoria, la cual es a su vez un dominio arbitrario, la derivada total de la integral
en dicha superficie obedece la siguiente expresión
Z  
d
Z
∂ ψs
ψs dA = + ∇s · (ψs vs ) + ψs (∇s · nβ γ )(w · nβ γ ) dA (2.7.6)
dt ∂t
Aβ γ Aβ γ

donde:

∇s = ∇ · P; P = I − nβ γ nβ γ (2.7.7)

es el operador gradiente superficial y P es el tensor de proyección tangencial a la superficie

divisoria. Además, el producto ∇s · nβ γ = 2H, donde H es la curvatura media de la superficie
divisoria (Slattery y col., 2006).

Usando este teorema, la ecuación (2.7.5) toma la siguiente forma

Z 
∂ Ts
ρsCps + ∇s · (ρsCps Ts vs ) + 2HρsCps Ts (w · nβ γ )
∂t
Aβ γ



−nβ γ · ρβ Cpβ Tβ (vβ − w) − ργ Cpγ Tγ (vγ − w) + qβ − qγ − Φs dA = 0 (2.7.8)

De donde es posible extraer la siguiente condición de frontera interfacial

 2.7 Condiciones de frontera 37

Condición de salto del flux calor

∂ Ts
ρsCps + ∇s · (ρsCps Ts vs ) + 2HρsCps Ts (w · nβ γ )
| {z ∂t} | {z } | {z }
Acumulación
Convección superficial Efecto del cambio de la superficie
superficial


= nβ γ · ρβ Cpβ Tβ (vβ − w) − ργ Cpγ Tγ (vγ − w) − kβ ∇Tβ − kγ ∇Tγ + Φs
| {z } |{z}
Fuente
Fluxes (convectivo y conductivo) desde las fases hacia la superficie superficial

(2.7.9)
Aquí se usó la ley de Fourier (ecuación 2.5.3) en el flux conductivo de calor de cada fase y
se incluyó el significado físico de cada término. Como puede notarse, la ecuación anterior es una
condición de salto (es decir, de discontinuidad) para el flux de calor.

2.7.2 Condición de salto de la temperatura

El complemento a la condición de frontera del flux de calor es la ley de enfriamiento de Newton,
la cual establece una discontinuidad entre la temperatura de una superfice y la de una fase fluida
circundante. Al lector interesado en la historia del desarrollo de esta condición de frontera se le
recomienda revisar el trabajo de Besson (2012). La ley de enfriamiento de Newton puede expresarse
en forma matemática como sigue:

Ley de enfriamiento de Newton

nβ γ · qβ = hβ γ Tβ − Tγ (2.7.10)

donde hβ γ [=] W/m2 K es el coeficiente interfacial de transferencia de calor. Esta condición de

frontera puede deducirse a partir del estudio de las capas límite hidrodinámica y térmica como
se muestra en capítulos posteriores. De hecho, es posible deducir expresiones para predecir al
coeficiente de salto hβ γ o bien, es posible usar formulaciones basadas en la observación (ver
Capítulo 14 de Bird y col., 2006). Este coeficiente es inversamente proporcional a las resistencias
interfaciales a la transferencia de calor. De esta forma, cuando este coeficiente tiende a cero, es
porque existen considerables resistencias a la transferencia de calor y la ecuación anterior se reduce
a que el flux interfacial es nulo. Por el contrario, cuando este coeficiente es suficientemente grande,
la ecuación anterior se reduce a una condición de continuidad de la temperatura entre las fases. De
acuerdo con Bergman y col. (2017), es posible combinar la ley de enfriamiento de Newton con la
ecuación de Stefan-Boltzmann para el caso en que la superficie emita calor, el resultado es
nβ γ · qβ = hβ γ Tβ − Tγ + εσ (Tβ4 − Tγ4 )


(2.7.11)

2.7.3 Clasificación de condiciones de frontera desde el punto de vista matemático

Las condiciones de frontera presentadas arriba se pueden reducir a menudo a alguna de las
siguientes tres versiones (Arfken y col., 2012):
1. Condiciones tipo Dirichlet: Se especifica el valor de la variable dependiente en alguna
posición. Por ejemplo: Considere una placa de metal de largo L que se localiza en medio de
otras dos placas que se encuentran a 25o C en un extremo (x = 0) y a 80o C en el otro (x = L)
(ver figura 2.3). Esta situación se expresa de la siguiente forma:
En x = 0, T = 25 ◦ C
En x = L, T = 80 ◦ C
donde se especifica el valor da la variable dependiente, en este caso, la temperatura. Note
que el transporte de calor en las interfaces se da por conducción.
 2. Condiciones tipo Neumann: Se especifica el valor de la derivada de la variable
Breviario sobre modelado Francisco J. Valdés Parada
Septiembre,
dependiente en alguna posición. Por ejemplo, si la placa del ejemplo 2012 se
anterior

encuentra
localizaaislada (esdedecir,
en medio otras su
dosdensidad deseflujo
placas que de calora 25°C
encuentran es cero) en extremo
en un un extremo
(x
38 Capítulo 2. Deducción de la ecuación de energía
= 0) y a 80°C en el otro (x = L) (Figura 1.6-1). Esta situación se expresa de la

siguiente forma:

En x = 0, T = 25 °C

En x = L, T = 80 °C

2. Condiciones tipo Neumann: Se especifica el valor de la derivada de la variable

dependiente en alguna posición. Por ejemplo, si la placa del ejemplo anterior se

Figura 1.6-1: Esquema de una placa de largo L en contacto con otras dos placas
encuentra de
Figura 2.3: Esquema aislada
una (es decir,
placa desu densidad
largo L en de flujo de con
contacto calorotras
es cero)
dosenplacas
un extremo
que se encuentran a
que se
temperaturas 25o C y 80oaC,
deencuentran temperaturas de 25°C y 80°C, respectivamente.
respectivamente.

2. Condiciones tipo Neumann: Se especifica el valor de la derivada de la variable dependiente en

alguna posición. Por ejemplo, suponga que ahora la placa del ejemplo anterior se encuentra
aislada (es decir, su densidad de flujo de calor es cero) en un extremo (x = 0) y en el otro
está expuesta a una densidad de flujo de calor constante qL (ver figura 2.4). En este caso, las
condiciones de frontera son:
dT
En x = 0, k =0
dx
Figura 1.6-1: Esquema de una placa de largo L en contacto con otras dos placas
dT
que xse=encuentran
En −k a temperaturas
L, Esquema = qL de 25°C y 80°C, respectivamente.
Figura 1.6-2: dx de una placa de largo L en contacto en un extremo con
una placa aislante y en el otro con una fuente de calor constante.

42
Figura 1.6-2: Esquema de una placa de largo L en contacto en un extremo con
Figura 2.4: Esquema de una placa de largo L en contacto en un extremo con una placa aislante y en
unafuente
el otro con una placa aislante
de calory en el otro con una fuente de calor constante.
constante

3. Condiciones de frontera mixtas (o tipo Robin o de Cauchy): En una condición de frontera de

tipo Cauchy, se especifica una condición tipo Dirichlet en una frontera y una tipo Neumann
en la otra frontera (o viceversa), de ahí el nombre de condiciones mixtas (Lapidus y Pinder,
1999). Por otro lado, una condición de frontera es tipo Robin cuando se especifica el valor
de una combinación lineal de los valores de una función y su derivada. En su forma general
42
unidimensional, una condición tipo Robin se expresa como sigue:
df
En x = a, α f (a) + β =γ (2.7.12)
dx x=a

donde α, β y γ son coeficientes conocidos que podrían (o no) ser funciones de la variable
dependiente f . Note que, si β = 0, la ecuación (2.7.12) se reduce a una condición de frontera
tipo Dirichlet; mientras que, si α = 0, se reduce a una condición tipo Neumann.
Un ejemplo típico de condición de frontera mixta es la ley de enfriamiento de Newton
discutida arriba (ecuación 2.7.10). Por ejemplo, suponga que la placa de los ejemplos
anteriores está expuesta en un extremo (x = L) a una corriente de aire cuya temperatura es
 Figura
2.8 Números 1.6-2: Esquema de
adimensionales una placa de
asociados a largo L en contacto en
la transferencia deuncalor
extremo con 39
una placa aislante y en el otro con una fuente de calor constante.
constante, T∞ y en el otro está en contacto con una placa a temperatura constante (ver figura
2.5). En este caso, el calor se transporta de la superficie de la placa hacia dentro de ésta
por difusión, mientras que de la superficie hacia afuera se transporta por convección1 . Esta
(x = se
situación 0) expresa
y en el otro
como está expuesta a una densidad de flujo de calor constante Q,
sigue:

las condiciones de fronteradTson:

En x = L, n · q = −k = h ( T |x=L − T∞ ) (2.7.13)
dx x=L | {z }
dT
convección
x = 0, −k = 0
| {z }
En
conducción
dx
Por supuesto, la ecuación (2.7.10) puede expresarse en la misma forma que la ecuación
dT
(2.7.12) con una simple manipulación
En x = algebraica.
L, −k = Q
dx

Figura 1.6-3: Esquema de una placa de largo L en contacto en un extremo con

Figura 2.5: Esquema de una placa de largo L en contacto en un extremo con una placa a 80o C y en
el otro con un
unaflujo
placadea aire
80°Ca yT∞
en. el otro con un flujo de aire a T . ∞

Ejercicio 2.9 — Continuidad del flux. Explique bajo qué condiciones, la condición de salto
dada en la ecuación (2.7.9) se reduce a la condición de continuidad del flux conductivo:

nβ γ · kβ ∇Tβ = nβ γ · kσ ∇Tσ (2.7.14)

De preferencia proporcione argumentos pasados en situaciones físicas y en análisis de órdenes

de magnitud. 42 

Ejercicio 2.10 — Forma alternativa de la ley de enfriamiento de Newton. Sin imponer

ninguna suposición, demuestre que la ecuación (2.7.11) puede escribirse como sigue


nβ γ · qβ = hT Tβ − Tγ (2.7.15)

Explique la forma matemática y el significado físico del coeficiente hT . 

Ejercicio 2.11 — Evaporación y condensación. Utilice los modelos matemáticos vistos en

este capítulo y utilícelos para describir los procesos de evaporación y condensación. 

2.8 Números adimensionales asociados a la transferencia de calor

Retornando la atención a la ecuación de energía térmica escrita en la ecuación (2.5.5); en esta
sección se busca reformularla de manera adimensional. Para ello, se definen las siguientes variables
1 Enrealidad el calor se transporta por convección y difusión, pero en muchas situaciones puede suponerse que la
convección domina sobre la difusión.
40 Capítulo 2. Deducción de la ecuación de energía

adimensionales
T − Tmı́n r t v
T∗ = ; r∗ = ; t∗ = ; v∗ = (2.8.1)
Tmáx − Tmı́n Lre f tre f vre f

donde Tmı́n y Tmáx representan los valores mínimo y máximo de las temperaturas de referencia,
mientras que Lre f y tre f son la longitud y el tiempo de referencia del proceso bajo estudio; por
último, vre f es la magnitud de la velocidad de referencia. Usando estas definiciones, la forma
adimensional de la ecuación de energía térmica es

1 ∂T∗
+ Pev∗ · ∇∗ T ∗ = ∇∗2 T ∗ + Φ∗ (2.8.2)
Fo ∂t ∗
Para obtener esta ecuación, se supuso que la conductividad térmica es constante y además se definió

2
ΦLre f
Φ∗ = (2.8.3)
(Tmáx − Tmı́n )k

junto con los números adimensionales de Fourier y de Péclet

kt ∗ vre f LρCp
Fo = 2
; Pe = (2.8.4)
ρCp Lre f k

El número de Fourier representa la razón del transporte difusivo de calor respecto a la tasa de
almacenamiento de calor; mientras que el número de Péclet representa la tasa de transferencia de
calor por convección, respecto al transporte difusivo de calor. Note que las definiciones anteriores
se expresan en términos del transporte difusivo de calor en lugar del transporte conductivo. Esto
se debe a que las expresiones de Fo y Pe están escritas en términos de la difusividad térmica
αT = k/ρCp .
De esta forma, para casos en los que Fo  1, no puede descartarse el término de acumulación,
mientras que para Fo  1 puede suponerse que el transporte de calor se da en estado estacionario.
Por su parte, Pe  1 implica que puede descartarse el transporte por convección respecto al
transporte difusivo, mientras que lo opuesto ocurre para Pe  1. En este momento es conveniente
definir al número de Prandtl como
µ/ρ
Pr = (2.8.5)
k/(ρCp )

Este número depende completamente de las propiedades del fluido y representa la capacidad de un
fluido para disipar cantidad de movimiento de manera viscosa respecto a su capacidad para difundir
energía térmica por conducción. Recordando que el número de Reynolds se define como
ρvre f Lre f
Re = (2.8.6)
µ

Entonces, el número de Péclet puede expresarse como

Pe = RePr (2.8.7)

Dirigiendo la atención a la ley de enfriamiento de Newton entre un sólido (fase-σ ) y un fluido

(fase-β ), la ecuación (2.7.10) se escribe como


−nσ β · kσ ∇Tσ = hσ β Tσ − Tβ (2.8.8)
 2.9 Resumen 41

donde además se sustituyó la ley de Fourier. Reescribiendo la ecuación anterior en términos de

variables adimensionales se obtiene
 
−nσ β · kσ ∇∗ Tσ∗ = hσ β Lre f Tσ∗ − Tβ∗ (2.8.9)

Existen dos maneras de completar la adimensionalización de la ecuación anterior, la primera

consiste en dividir ambos lados entre kσ
 
−nσ β · ∇∗ Tσ∗ = Bi Tσ∗ − Tβ∗ (2.8.10)

donde se definió el número de Biot como

hσ β Lre f
Bi = (2.8.11)
kσ
La segunda forma de adimensionar la ecuación (2.8.9) consiste en dividir ambos lados entre la
conductividad térmica del fluido, kβ , el resultado es
 
−nσ β · kσ β ∇∗ Tσ∗ = Nu Tσ∗ − Tβ∗ (2.8.12)

donde kσ β = kσ /kβ es la razón de conductividades térmicas y el número de Nusselt se define como

hσ β Lre f
Nu = (2.8.13)
kβ
Comparando las definiciones de los números de Biot y de Nusselt, resulta la siguiente expresión

Nu = kσ β Bi (2.8.14)

Como puede notarse tanto el número de Biot como el número de Nusselt relacionan la transferencia
de calor por convección interfacial respecto al transporte de calor por conducción, ya sea del lado
del sólido (Biot) o del fluido (Nusselt). Gracias a estos números adimensionales es fácil identificar
una condición en la que existen considerables resistencias interfaciales a la transferencia de calor
como Bi  1 (o bien Nu  1), mientras que Bi  1 (o bien Nu  1) corresponde a una situación en
la que no hay resistencias interfaciales y por lo tanto puede suponerse continuidad de la temperatura,
como se explicó anteriormente. En la práctica, se utiliza el número de Nusselt para predecir al
coeficiente interfacial de transferencia de calor a partir de los números de Reyndolds y Prandtl. Un
ejemplo de ello es la ecuación de Ranz y Marshall (1952)
√ √ 3
Nu = 2 + 0.6 Re Pr, 0 ≤ Re ≤ 200; 0 ≤ Pr < 250 (2.8.15)

Aunque existen muchas otras correlaciones como lo discuten Aissa y col. (2015) y Bird y col.
(2006). De manera general Bird y col. (2006, ver tabla 14.1-1) muestran que el coeficiente h es
típicamente mayor en líquidos que en gases tanto para condiciones de convección libre como para
convección forzada.
Ejercicio 2.12 ¿ Para qué condición del número de Biot puede suponerse que la temperatura
del sólido será uniforme ? 

2.9 Resumen
En este capítulo se dedujeron tres ecuaciones de transporte de energía. La más completa es la
ecuación de energía total:
 
D 1
ρ e + v2 = −∇ · (q + qR ) + ∇ · (T · v) + ρg · v (2.9.1)
Dt 2
42 Capítulo 2. Deducción de la ecuación de energía

La cual es el resultado de sumar la ecuación de energía mecánica:

ρ Dv2
= ρg · v + ∇ · (T · v) − T : ∇v (2.9.2)
2 Dt
y de energía térmica
De
ρ = −∇ · (q + qR ) + T : ∇v (2.9.3)
Dt
Esta última será la ecuación de trabajo durante este curso y puede expresarse en términos de la
temperatura mediante la ley de Fourier como:

∂T
ρCp + ρCp v · ∇T = ∇ · (k∇T ) + Φ (2.9.4)
| {z∂t} | {z } |Conducción {z } |{z}
Convección Fuente o sumidero
Acumulación

donde en el término fuente puede estar involucrado el transporte por radiación.

En este momento es conveniente hacer una distinción entre transferencia y transporte de calor.
En el segundo caso se involucra tanto la conducción como la convección, mientras que en el primer
caso, sólo se involucran términos conductivos, por lo que la ecuación de energía térmica se reduce
a:
∂T
ρCp = ∇ · (k∇T ) + Φ (2.9.5)
| {z∂t} | {z } |{z}
Conducción Fuente o sumidero
Acumulación

La cual en estado estacionario se reduce a:

0 = ∇ · (k∇T ) + Φ
|{z} (2.9.6)
| {z }
Conducción Fuente o sumidero

En el próximo capítulo se abordan problemas de aplicación de esta forma simplificada de la

ecuación de energía térmica.
Junto con la ecuación de energía térmica, se presentaron las correspondientes condiciones de
salto tanto para el flux de calor como para la temperatura y se explicó que en muchos casos pueden
reducirse a formas sencillas dadas por las condiciones de frontera tipo Dirichlet, Neumann o Robin.
De hecho, en general las condiciones de salto son de tipo Robin.
II
Parte Dos: Aplicaciones

3 Conducción de calor en estado estacio-

nario . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.1 Introducción
3.2 Dos alternativas para simplificar modelos
3.3 Conducción de calor en una placa
3.4 Conducción unidimensional en estado estaciona-
rio en coordenadas cilíndricas
3.5 Conducción de calor en coordenadas esféricas
3.6 Una introducción al escalamiento: La aleta de en-
friamiento
3.7 Aletas de enfriamiento cilíndricas
3.8 Aletas de enfriamiento anulares
4 Transferencia de calor en más de una di-
rección . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
4.1 El problema de Sturm-Liouville
4.2 Coordenadas rectangulares
4.3 Conducción de calor transitoria en una barra cilín-
drica
4.4 Transporte en medios semi-infinitos
5 Soluciones numéricas . . . . . . . . . . . . . . . 95
5.1 Coordenadas cilíndricas
5.2 Problema de aplicación
5.3 Aleta de enfriamiento revisada
5.4 Conducción de calor transitoria
5.5 Conducción y convección en estado transitorio
6 Capa límite . . . . . . . . . . . . . . . . . . . . . . . . 107
6.1 Ecuaciones de capa límite de energía y cantidad
de movimiento para flujo laminar
6.2 Capa límite cuando se especifica el flux de calor
en la superficie
7 Modelos macroscópicos . . . . . . . . . . . . 117
7.1 Introducción
7.2 Transferencia de calor en un panel fotovoltaico
7.3 Tanque agitado homogéneo
7.4 Tanque agitado no homogéneo
7.5 Transferencia de calor en un lecho fijo
7.6 Transferencia de calor en intercambiadores de ca-
lor de doble tubo
7.7 Transferencia de calor en un pozo geotérmico
Bibliografía . . . . . . . . . . . . . . . . . . . . . . . . . 134
 3. Conducción de calor en estado estacionario

3.1 Introducción
En este capítulo se estudian algunas aplicaciones de la ecuación de energía térmica y de las
condiciones de frontera para condiciones en las que el transporte se da en estado estacionario y
los efectos del transporte convectivo son despreciables respecto al transporte por conducción. En
términos de números adimensionales, lo anterior equivale a restringir el análisis para condiciones
en las que 1  Fo y Pe  1. De esta forma, el único modo de transferencia de calor remanente
es la conducción; si además se supone constante la conductividad térmica, la ecuación de energía
térmica se reduce a la ecuación de Poisson

k∇2 T = −Φ (3.1.1)

Para condiciones en las que Φ = 0, la ecuación anterior se reduce a la ecuación de Laplace

∇2 T = 0 (3.1.2)

La versión transitoria de esta última ecuación se conoce, en física y matemáticas, como la ecuación
de calor y se estudia en capítulos posteriores de este texto. Físicamente estas condiciones de
transporte son aplicables a sólidos y a fluidos bajo condiciones hidrostáticas. De lo expuesto arriba
cabe notar que incluso bajo la forma de la ecuación de Laplace, la forma simplificada de la ecuación
de energía térmica es una ecuación diferencial parcial. Con la intención de reducirla a una ecuación
diferencial ordinaria se pueden usar simplificaciones basadas en análisis de orden de magnitud
o bien basadas en un proceso de promediado. Ambas propuestas son de utilidad y es importante
saber elegir cuándo elegir una u otra. Por ello, antes de comenzar con aplicaciones particulares se
explicará en qué consisten estas simplificaciones.

3.2 Dos alternativas para simplificar modelos

3.2.1 Simplificaciones basadas en estimados de orden de magnitud
El orden de magnitud de una función se define como sigue:
46 Capítulo 3. Conducción de calor en estado estacionario

Definición 3.2.1 — Orden de magnitud. Sea f una función definida en R3 , su orden de

magnitud es
1
Z
O( f ) = k f k dV (3.2.1)
V
V
R
En donde V = dV .
V

Esta definición también es aplicable a funciones definidas en R2 y R1 y puede usarse para escalares,
vectores y tensores. Sin embargo, no siempre se cuenta con los valores de la función para poder
calcular su orden de magnitud. En estos casos se recurre a estimados de orden de magnitud, los
cuales se definen como
Definición 3.2.2 — Estimado de orden de magnitud. Sea f una función definida en R3 , su
estimado de orden de magnitud corresponde a su valor promedio y está dado por
1
Z
O( f ) = f dV = h f i (3.2.2)
V
V

De esta forma, el estimado de orden de magnitud de una primera derivada espacial en x definida
en un dominio x ∈ [a, b], es:
x=b
f (b) − f (a)
 
∂f 1 ∂f
Z
O = dx = (3.2.3)
∂x b−a ∂x b−a
x=a

Whitaker (sección 2.9 1983) sugiere adoptar las siguientes suposiciones:

f (b) = O( f ); f (a) = O( f ) (3.2.4)
Dado que cuando se estiman ordenes de magnitud, los signos positivo y negativo no se toman en
cuenta, la ecuación (3.2.3) puede expresarse como
 
∂f f
=O (3.2.5)
∂x b−a
Este resultado puede extenderse a una derivada de cualquier orden para dar
∂n f
 
f
=O (3.2.6)
∂ xn (b − a)n
De esta forma, si se tiene la ecuación de Laplace para la temperatura en tres direcciones en
coordenadas cartesianas:
∂ 2T ∂ 2T ∂ 2T
0= + 2 + 2 (3.2.7)
∂ x2 ∂y ∂z
no se pueden calcular los ordenes de magnitud de las derivadas pues no se conoce la forma de la
función T , pero en su lugar se pueden proponer los siguientes estimados de orden de magnitud
∂ 2T
 
T
=O (3.2.8a)
∂ x2 L2
∂ 2T
 
T
= O (3.2.8b)
∂ y2 W2
∂ 2T
 
T
2
=O (3.2.8c)
∂z H2
3.2 Dos alternativas para simplificar modelos 47

Esta situación corresponde a la transferencia de calor por conducción en estado estacionario en

una placa de largo L, ancho W y alto H. Note que para proponer estos estimados se requiere que el
sistema coordenado se ubique en una de las esquinas de la placa. De esta forma si, por ejemplo, el
largo es mucho mayor que el ancho y el largo, se tiene que

W 2  L2 , H 2  L2 (3.2.9)

Con base en esta separación de longitudes características de la placa y tomando en cuenta los
estimados de orden de magnitud dados en las ecuaciones (3.2.8), es razonable concluir que

∂ 2T ∂ 2T ∂ 2T ∂ 2T
 ;  (3.2.10)
∂ x2 ∂ y2 ∂ x2 ∂ z2
Note que para propósitos prácticos sólo es necesaria una de las dos anteriores desigualdades para
deducir que la ecuación (3.2.7) puede reducirse a

∂ 2T ∂ 2T
0= + 2 (3.2.11)
∂ y2 ∂z
En este momento vale la pena recapitular el procedimiento para llevar a cabo una simplificación
por medio de estimados de orden de magnitud

Pasos a seguir para simplificar una ecuación mediante estimados de orden de magnitud
1. Utilice las fórmulas dadas en las ecuaciones (3.2.2) y (3.2.6) para estimar los órdenes de
magnitud de las distintas derivadas.
2. Con base en información de las dimensiones físicas del sistema, establezca una restricción
de escala.
3. Con el sustento de la restricción de escala y de los estimados del paso 1, proponga, si es
el caso, una suposición que lleve a descartar un término.
4. Elimine al término que se estimó como el más pequeño de la ecuación original.

En este punto, vale la pena cuestionar ¿cuál es la ganancia de llevar a cabo estos pasos y por qué
no simplemente descartar el término con base en la intuición? La respuesta esta pregunta está en
que al conocer las restricciones de escala, se tiene una idea mucho más clara sobre el rango de
aplicabilidad del modelo resultante que con base en la mera intuición. Whitaker (1988) es más
sistemático en este tema y propone tres niveles de simplificación: En el primer nivel simplemente
se elimina un término con base en la intuición. En el segundo nivel se hace una suposición como
la presentada en las desigualdades (3.2.9). En el tercer nivel se encuentran las restricciones de
escala, las cuales sustentan las suposiciones y llevan a la eliminación del término en cuestión como
consecuencia lógica.

Ejercicio 3.1 — Orden de magnitud vs. estimado de orden de magnitud pt. 1. Compare
el orden de magnitud y el estimado de orden de magnitud para las siguientes funciones definidas
en el intervalo x = 0 a x = L:
1. f (x) = c
2. f (x) = cx
3. f (x) = cx 2

2cx, 0 6 x 6 L/2
4. f (x) =
cL, L/2 6 x 6 L

 48 Capítulo 3. Conducción de calor en estado estacionario

Ejercicio 3.2 — Orden de magnitud vs. estimado de orden de magnitud pt. 2. Repita
el ejercicio anterior para d f /dx. 

Ejercicio 3.3 Discuta sobre la diferencia, si es que la hay, entre estimar el orden de magnitud
del gradiente de un vector y su divergencia. 

3.2.2 Simplificaciones basadas en un proceso de promediado

Como puede notarse las simplificaciones mediante estimados de orden de magnitud se funda-
mentan en restricciones de escala y no necesariamente consideran la información de las condiciones
de frontera. Cuando la región de promediado es tal que abarca a alguno(s) de los límites del dominio,
entonces se puede reducir la dimensión espacial del problema sin recurrir a ninguna aproximación.
Esto es posible ya que la operación matemática de promediado es una transformación integral
que lleva a un reacomodo de información, más no a una pérdida de la misma (Wood, 2009; Wood
y Valdés-Parada, 2013). Por ejemplo, si se desea eliminar la dependencia en la dirección x de la
temperatura en la ecuación (3.2.7), se puede definir el siguiente operador de promediado

x=L
1
Z
x
hT i = T dx (3.2.12)
L
x=0

donde el superíndice x se usa para indicar la dirección en la cual se está llevando la operación de
promediado. Note que hT ix es función de y y z, pues se eliminó la dependencia con x. Aplicando
este operador a la ecuación (3.2.7) da como resultado

∂ 2 hT ix ∂ 2 hT ix
 
1 ∂T ∂T
0= − + + (3.2.13)
L ∂x x=L ∂x x=0 ∂ y2 ∂ z2

Si se conoce el valor de la derivada de la temperatura en x = 0 y x = L, pueden sustituirse

directamente estos valores en la ecuación diferencial y tratar ahora con una ecuación diferencial
cuya variable dependiente es hT ix . Como puede notarse, esto implica que la información que se
encontraba en las fronteras ahora forma parte de la ecuación diferencial. Esta es una muestra del
reacomodo de información que se mencionó arriba. Sin embargo, si por ejemplo, se contara con
condiciones de tipo Dirichlet en las fronteras x = 0 y x = L, entonces este procedimiento no ayuda
mucho y se recomendaría más usar el método basado en estimados de ordenes de magnitud. Estos
dos tipos de formas de simplificación se usan a lo largo del capítulo para reducir la dimensión de
las ecuaciones diferenciales.
Para cerrar esta sección vale la pena comentar que hay casos en los que es conveniente que
la región de promediado tenga una longitud característica, r0 que sea mas grande que la longitud
característica del continuo, `, y a la vez mucho menor que la longitud característica de todo el
sistema, L; esto es

`  r0  L (3.2.14)

Cuando se aplica un operador de promediado con estas características a sistemas multifásicos,

en ocasiones, es posible definir regiones representativas (ver sección 1.1 en Bear, 2018). Cuando
este es el caso, se pueden ver a este tipo de sistemas como nuevos medios continuos y este es el
fundamento de técnicas de promediado como homogenización (Auriault y col., 2009), promedio
volumétrico (Whitaker, 1999), la teoría de promediado restringida por la termodinámica (Gray y
Miller, 2014) y el método de homogenización adjunta (Bottaro, 2019).
 3.3 Conducción de calor en una placa 49

3.3 Conducción de calor en una placa

3.3.1 Temperaturas superficiales especificadas
Considere la transferencia de calor por conducción en estado estacionario en una placa de largo
L, profundidad W y altura H como la mostrada en la figura 2.3. Esta placa no tiene fuentes de
calentamiento pero mantiene una cara (en x = 0) a temperatura constante (T0 ) y la cara opuesta (en
x = L) a una temperatura TL < T0 . Lo anterior se traduce en que los extremos en la dirección x no
ofrecen resistencia a la transferencia de calor. Además, las otras caras de la placa se mantienen
aisladas.
Para las condiciones descritas arriba, la ecuación de energía térmica se reduce a la ecuación de
Laplace, la cual en coordenadas cartesianas se escribe como sigue:
∂ 2T ∂ 2T ∂ 2T
0= + 2 + 2 (3.3.1a)
∂ x2 ∂y ∂z
Más aún, de acuerdo a la descripción de lo que ocurre en las superficies de la placa, se tienen las
siguientes condiciones de frontera, las cuales son de tipo Dirichlet y Neumann (ver sección 2.7.3):

En x = 0, T = T0 (3.3.1b)
En x = L, T = TL (3.3.1c)
∂T
En y = 0, H =0 (3.3.1d)
∂y
∂T
En z = 0,W =0 (3.3.1e)
∂z
De acuerdo a la descripción del problema, la temperatura máxima se encuentra en x = 0, y su valor
mínimo está en x = L. Esto sugiere que los cambios más importantes de la temperatura ocurren en
la dirección x. Ya que se conocen las condiciones de frontera para el flux de calor en las direcciones
y y z, la opción más razonable para reducir la dimensión del problema es utilizar un proceso de
promediado. Con este fin, se define el siguiente operador integral
Z y=W
z=H
1
Z
yz
hT i = T dydz (3.3.2)
HW
z=0 y=0

En donde hT iyz sólo es función de x. Aplicando este operador de promediado a la ecuación

diferencial se obtiene
z=H ! y=W
d 2 hT iyz
Z  
1 ∂T ∂T 1 ∂T ∂T
Z
0= + − dz+ − dy (3.3.3)
dx2 HW ∂ y y=W ∂ y y=0 HW ∂ z z=H ∂ z z=0
z=0 y=0

Note que en el primero término fue posible intercambiar diferenciación e integración ya que los
límites de la integral no dependen de x. La ecuación anterior puede simplificarse tomando en cuenta
las condiciones de frontera dadas en las ecuaciones (3.3.1d) y (3.3.1e) a la siguiente expresión
d 2 hT iyz
0= (3.3.4a)
dx2
Para encontrar las condiciones de frontera correspondientes es necesario aplicar el operador de
promediado definido en la ecuación (3.3.2) a las ecuaciones (3.3.1b) y (3.3.1c) para obtener

En x = 0, hT iyz = T0 (3.3.4b)
yz
En x = L, hT i = TL (3.3.4c)
 Heat and Mass Transfer Fundamentals
Fall, 2009

50 Capítulo 3. Conducción de calor en estado estacionario

Figura 3.1: Evaluación del perfil de temperatura en una placa con temperaturas fijas en sus extremos.
Figure 2.4. Temperature profiles in a slab with specified temperatures.

Aquí se supuso que las temperaturas T0 y TL son constantes en las direcciones y y z.

De esta Eqs.
forma (2.3-6) andconseguido
se ha (2.3-7) can besimplificar
used to evaluate the constants of integration,
considerablemente A and B,para
el problema and obtener una
the temperature field is given by
ecuación diferencial ordinaria, lineal y homogénea sujeta a dos condiciones de frontera tipo
Dirichlet. La solución de este problema de valor a la frontera x es
T = T0 − (T0 − T1 ) (2.3-9)
L
x
hT iyz = T0 − (T0 − T1 ) (3.3.5)
Most often the engineer
L is interested in this result because it can be used to compute the
heat flux or the rate at which energy is transferred across the flat plate. This is given by
La evaluación de la solución se muestra en la Figura 3.1 en donde se puede apreciar que la pendiente
del perfil de temperatura aumenta aqmedida ∂T aumenta
que ⎛ k ⎞ la diferencia de temperaturas entre los dos
x = −k = ⎜ ⎟ (T0 − T1 ) (2.3-10)
extremos. Para colapsar estos resultados en una ∂x sola⎝línea
L ⎠ se proponen las siguientes definiciones
de variablesand
adimensionales
we see that the heat flux is constant everywhere in the slab. Shorter is better. and it is
directly
hT iyz −proportional
T1 to the
x thermal conductivity and the temperature difference at the walls
x∗ = to the slab thickness. It is helpful to think of Eq. (2.3-10) in the
T ∗ = and inversely;proportional (3.3.6)
T0 − T1 L
form
De esta forma, la ecuación (3.3.5) puede escribirse como
driving force
heat flux = (2.3-11)
resistance
T ∗ = 1 − x∗ (3.3.7)

Lo anterior sirve para ilustrar la utilidad de resolver el problema en términos de variables adimen-
sionales. El uso de estas variables se explora en17ejemplos posteriores y también en el ejercicio
3.4.
La temperatura promedio en toda la placa se define como
x=L
1
Z
xyz
hT i ≡ hT i = hT iyz dx (3.3.8)
L
x=0

Sustituyendo la ecuación (3.3.5) en la definición anterior, se obtiene el siguiente resultado

T1 + T0
hT i = (3.3.9)
2
Es decir, la temperatura promedio en toda la placa corresponde a la media aritmética.
Por último, el flux de calor por conducción en la dirección x es:
∂T
qx = −k (3.3.10)
∂x
 3.3 Conducción de calor en una placa 51

O bien, al promediar en las direcciones y y z:

dhT iyz T0 − T1 fuerza impulsora

hqx iyz = −k = = (3.3.11)
dx L/k resistencia

Donde es evidente que, a mayor longitud (o menor conductividad térmica), habrá mayores resisten-
cias a la transferencia de calor. Como comentario final de este problema, vale la pena mencionar
que la solución dada en la ecuación (3.3.5) satisface el principio del máximo para la ecuación de
Laplace:

Principio del máximo: Cuando la temperatura satisface la ecuación de Laplace, el valor máximo
(o mínimo) que puede adquirir no se encuentra en el interior del dominio, sino en las fronteras.

Este principio aplica también a la formulación original del problema dada en las ecuaciones (3.3.1).
Este principio matemático no aplica en el siguiente caso de estudio.

Ejercicio 3.4 — Problema adimensional. Utilice las definiciones dadas en las ecuaciones
(3.3.5) para escribir las ecuaciones (3.3.4) en forma adimensional. Resuelva el problema y
recupere la solución dada en la ecuación (3.3.7). Calcule además la forma adimensional de la
temperatura promedio y del flux de calor.


Ejercicio 3.5 — Condiciones tipo Cauchy. Reconsidere este problema para el caso en que
en x = L en lugar de especificar la temperatura, se especifica el valor del flux de calor, q1 . 

3.3.2 Placa con generación de calor uniforme

Reconsidere el problema de la placa cuyas caras se encuentran a temperaturas constantes con la
diferencia de que ahora la placa contiene una resistencia eléctrica interna. Por lo tanto, el problema
de valor a la frontera es:

∂ 2T ∂ 2T ∂ 2T Φ
0= + + + (3.3.12a)
∂ x2 ∂ y2 ∂ z2 k
En x = 0, T = T0 (3.3.12b)
En x = L, T = TL (3.3.12c)
∂T
En y = 0, H =0 (3.3.12d)
∂y
∂T
En z = 0,W =0 (3.3.12e)
∂z

O bien, tras aplicar el operador de promediado definido en la ecuación (3.3.2),

d 2 hT iyz hΦiyz
0= + (3.3.13a)
dx2 k
En x = 0, hT iyz = T0 (3.3.13b)
yz
En x = L, hT i = T1 (3.3.13c)

Dado que hΦiyz es constante, la solución de este problema es:

x hΦiyz L2  x x
hT iyz = T0 − (T0 − T1 ) + 1− (3.3.14)
L k 2 L L
 Π = 1 − Ω (1 − 2X ) (2.3-32)

Evaluation of Eqs. (2.3-29) and (2.3-32) is presented in Fig. 2.5 for both positive and
negative values of Ω; we have also included the results from taking Ω = 0 as a reference.
52 Capítulo 3. Conducción de calor en estado estacionario

a) b)
Figure
Figura 2.5a)a)Perfiles
3.2: Temperature and b) heat
de temperatura y b) flux distributions
del flux de calor eninforma
a flat adimensional,
plate with a heat
parasource (Ω
diferentes
valores de Ω. > 0) or sink (Ω < 0).

Interestingly, all the values of the dimensionless heat flux, Π, converge to a value of 1 at
O bien, definiendo las variables adimensionales:
the middle of the plate. This observation is confirmed by taking X = 0.5 in Eq. (2.3-32).
hT iyz − T1 ΦL2
= when Ω >0,
NoteT that ; xthe ; Ω = profiles exhibit a maximum and when Ω (3.3.15)
= temperature < 0, a
∗ ∗ x
T0 − T1 in the temperature
minimum appears L 2k(T0 − TThe
profiles. 1 ) position of the maximum or minimum

setemperature can be found

obtiene el siguiente by taking the derivative of Θ with respect to X in Eq. (2.3-29) and
resultado
finding the value of X that makes the resulting expression to be equal to zero; the resulting
T ∗ = [1 + Ωx∗ ] (1 − x∗ ) = 1 + (Ω − 1) x∗ − Ωx∗2 (3.3.16)
critical value is
Esta expresión es un polinomio de segundo orden cuyo único parámetro es el término fuente Ω.
1
La evaluación de esta expresión se muestra en la figura 1
X crit = − 3.2a). Como puede notarse, para valores
(2.3-33)
de Ω 6= 0, los perfiles incluyen un punto máximo (o2 mínimo)2Ω dentro del dominio. Lo cual está en
concordancia con el principio del máximo. El lugar donde se encuentra este punto crítico es:
∗ Ω−1 1 1
xcrit = = − (3.3.17)
2Ω 2 2Ω

Sustituyendo este resultado en la ecuación (3.3.16) lleva al siguiente resultado para la tempera-
tura máxima o mínima 21
 2
∗ Ω+1
Tmáx,mı́n =Ω (3.3.18)
2Ω
La temperatura promedio en la placa se define como
x=L
1
Z
hT i = hT iyz dx (3.3.19)
L
x=0

Por lo que al sustituir la ecuación (3.3.14) y llevando a cabo los pasos de integración correspon-
dientes, se obtiene el siguiente resultado
1 hΦiyz L2
hT i = T0 − (T0 − T1 ) + (3.3.20)
2 k 12
 3.3 Conducción de calor en una placa 53

O bien en forma adimensional

xZ∗ =1
1 Ω
hT ∗ i = T ∗ dx∗ = + (3.3.21)
2 6
x∗ =0

El flux de calor promedio se definió anteriormente como

dhT iyz
hqx iyz = −k (3.3.22)
dx
Por lo que en este caso el resultado es
(T0 − T1 )
 
yz yz L 2x
hqx i = − hΦi 1− (3.3.23)
L/k 2 L
Para escribir esta expresión totalmente en función de variables adimensionales se define
qx L
Π= (3.3.24)
k(T0 − T1 )
Por lo tanto,
dT ∗
Π=− = 1 − Ω(1 − 2x∗ ) (3.3.25)
dx∗
La evaluación de los perfiles del flux de calor adimensionales se muestra en la figura 3.2b) para
diferentes valores de Ω. Como puede notarse, todos los valores del flux de calor convergen a la
unidad en la mitad de la placa.

Ejercicio 3.6 — Criterio de la segunda derivada. Utilice el criterio de la segunda derivada

para demostrar que los valores críticos dados en la ecuación (3.3.17) corresponden a un máximo
o mínimo dependiendo del signo de Ω. 

Ejercicio 3.7 — Fuente dependiente de la posición. Repita los desarrollos de este problema
para el caso en que la fuente de calor dependa linealmente de la posición de acuerdo a:
 x
hΦiyz = Φ0 1 − (3.3.26)
L
donde Φ0 es el valor máximo de la fuente y es una constante. 

3.3.3 Transferencia de calor considerando convección externa

Una persona se encuentra en la calle en medio de una tormenta, donde la temperatura ambiente
está a Ta = 10 o C, mientras que en sus pulmones la temperatura del aire está a T` = 25 o C. En este
problema se busca predecir el perfil de temperatura desde el ambiente, pasando por las capas de
ropa y piel hasta llegar a los pulmones así como la tasa de transferencia de calor del cuerpo hacia el
exterior.
Este problema abarca en realidad una amplia cantidad de fases: el aire exterior, las capas de
ropa, la piel, músculos, huesos, las paredes de los pulmones y el aire dentro de los pulmones. Para
simplificar el tratamiento a las capas de ropa se les visualiza como una sola fase (fase-σ ) con una
conductividad térmica constante (kσ ) y espesor `σ . Además, la piel, músculos, huesos y las paredes
de los pulmones se consideran como otra fase (fase-λ ) con una conductividad kλ y espesor `λ . De
esta forma, el problema se reduce a cuatro fases, de las cuales los perfiles de temperatura sólo se
calcularán en las fases sólidas.
54 Capítulo 3. Conducción de calor en estado estacionario

Fijando el sistema coordenado para que comience en la pared de los pulmones y siguiendo
los desarrollos de los casos de estudio anteriores, las ecuaciones diferenciales que gobiernan la
temperatura promediada en los ejes y y z en las fases sólidas son

d 2 Tλ
En x ∈ (0, `λ ), 0= (3.3.27a)
dx2
d 2 Tσ
En x ∈ (`λ , `+
σ ), 0= (3.3.27b)
dx2
yz yz
donde para simplificar la nomenclatura se usaron `+ σ = `σ + `λ , Tλ = hT iλ y Tσ = hT iσ . Las
ecuaciones anteriores están sujetas a las siguientes condiciones de frontera:

dTλ
En x = 0, n · q = kλ = hλ ( Tλ |x=0 − T` ) (3.3.27c)
dx x=0
En x = `λ , Tλ = Tσ (3.3.27d)
dT dTσ
En x = `λ , kλ λ = kσ (3.3.27e)
dx dx
dTσ  
En x = `+
σ, n · q = −kσ = hσ Tσ |x=`+σ − Ta (3.3.27f)
dx x=`+
σ

La solución de este problema está dada por el siguiente par de ecuaciones

h   i
kλ T` k`σσ + k`λ + h1σ + hTa − (T` − Ta ) x
Tλ = λ λ
 (3.3.28)
`σ `λ
kλ hσ + kσ + k + h1
1
λ λ
+
 
`
x (T` − Ta ) − kσ h1σ + kσσ (T` − Ta )
Tσ = Ta −   (3.3.29)
kσ h1σ + k`σσ + k`λ + h1
λ λ

El flux de calor promedio es constante en todo el sistema y está dado por

T` − Ta fuerza impulsora
qλ = qσ = = (3.3.30)
1 `σ `λ 1 resistencia
+ + +
hσ kσ kλ hλ

donde qλ = hqx iyz yz

λ y qσ = hqx iσ . Entonces, si se denota a la superficie externa de la ropa en la
espalda de la persona como A, la tasa de transferencia de calor desde la persona hacia el ambiente
es:

Q̇ = A qσ |x=`+σ = UA (T` − Ta ) (3.3.31)

donde U es el coeficiente global de transferencia de calor y es el inverso de la suma de las

resistencias, esto es:

1 −1
 
1 `σ `λ
U= + + + (3.3.32)
hσ kσ kλ hλ

Esta expresión se asemeja a la resistencia equivalente de circuitos eléctricos en serie e indica que, si
el objetivo está en estudiar las resistencias a la transferencia de calor pueden agregarse o eliminarse
términos en U conforme sea necesario.
 3.4 Conducción unidimensional en estado estacionario en coordenadas
cilíndricas 55
aire
r2 − r1

r1 − r0

r0

Aislante 2 Aislante 1
Tubo a T0
Figura 3.3: Esquema de la sección transversal de un cilindro de paredes compuestas

Ejercicio 3.8 — Cálculo de la tasa de transferencia de calor. Tomando T` = 25 o C, Ta =

10 o C, `σ = 30 mm, `λ = 2 cm, hλ = hσ = 13.1 W/m2 K, kσ = 0.0364 W/mK, kλ = 0.282
W/mK.; calcule Q̇. 

Ejercicio 3.9 — Temperatura interna con un abrigo. Suponiendo que el valor de Q̇ del
problema anterior es el mismo si la persona se pone un abrigo de 3 cm de espesor y conductividad
térmica de 0.035 W/mK, calcule la temperatura interna del aire en los pulmones de la persona. 

3.4 Conducción unidimensional en estado estacionario en coordenadas ci-

líndricas
Considere un tubo de dos paredes compuestas como el mostrado en la Figura 3.3, el cual
en su compartimiento interior, lleva un tubo a temperatura constante T0 , la cual es mayor que la
temperatura del aire, Ta . Deduzca una expresión para el perfil de temperaturas en las paredes del
tubo bajo condiciones de estado estacionario. Además determine la tasa de transferencia de calor
en la superficie del tubo externo.
Dado que no hay una fuente de generación de calor en las paredes del tubo, la ecuación de
energía térmica en este caso se expresa como:

1 ∂ 2T ∂ 2T
 
1 ∂ ∂T
0= r + 2 + (3.4.1)
r ∂r ∂r r ∂ θ 2 ∂ z2
La cual puede aplicarse tanto en el tubo interno (dominio 1) como en el externo (dominio 2).
Antes de proseguir es conveniente analizar la posibilidad de simplificar esta ecuación diferencial.
Suponiendo que las superficies en z = 0 y z = H están aisladas, es decir,
∂T
En z = 0, H; =0 (3.4.2)
∂z
y que el transporte de calor es axisimétrico:
∂T ∂T
= (3.4.3)
∂θ θ =0 ∂θ θ =2π
56 Capítulo 3. Conducción de calor en estado estacionario

Determine el perfil de temperatura a través de las capas de aislante. Al igual que en el caso de
la placa, para simplificar el problema es conveniente definir una temperatura promedio. Como
en este caso se tienen condiciones de frontera homogéneas en las direcciones axial y angular, es
conveniente definir el siguiente operador de promediado:
θZ=2π z=H
1
Z
θz
hT i (r) = T (r, θ , z)dθ dz (3.4.4)
2πH
θ =0 z=0
Aplicando el operador de promediado a la ecuación diferencial en cada dominio y tomando en
cuenta las condiciones de frontera resulta que:
dhT1 iθ z
 
1 d
0= r ; r0 < r < r1 (3.4.5a)
r dr dr
dhT2 iθ z
 
1 d
0= r ; r1 < r < r2 (3.4.5b)
r dr dr
Las condiciones de frontera son:
En r = r0 , hT1 iθ z = T0 (3.4.5c)
θz θz
En r = r1 , hT1 i = hT2 i (3.4.5d)
dhT1 iθ z dhT2 iθ z
En r = r1 , − k1 = −k2 (3.4.5e)
dr dr
dhT2 iθ z
En r = r2 , − k2 = h1 (hT2 iθ z − Ta ) (3.4.5f)
dr
La solución del problema es
T0 − Ta r
hT1 iθ z = T0 −   ln (3.4.6a)
k1 r2 r1 r2 r2 1 r0
r2 k1 ln r0 + k2 ln r1 + h1
 
1 r2 r2
h1 + k2 ln r (T0 − Ta )
hT2 iθ z =   + Ta (3.4.6b)
r2 r1 r2 r2 1
k1 ln r0 + k2 ln r1 + h1

Con estos resultados se obtienen las siguientes expresiones para el flux de calor en cada material
aislante
T0 − Ta
hq1r iθ z = hq2r iθ z =   (3.4.7)
r2 r1 r2 r2 1 r
k1 ln r0 + k2 ln r1 + h1 r2

A diferencia del problema anterior, aunque el flux de calor es el mismo en ambos materiales, en
este caso no es constante, sino que disminuye con el radio. De esta forma, la tasa de transferencia
de calor en la superficie del segundo aislante es
2πr2 L (T0 − Ta )
Q̇ r2
=  (3.4.8)
r2
k1 ln rr10 + kr22 ln rr21 + h11
Este resultado puede expresarse en forma más compacta como sigue
Q̇ r2
= U2 A2 (T0 − Ta ) (3.4.9)
En donde A2 = 2πr2 L y U2 es el coeficiente global de transferencia de calor para la superficie
externa, el cual se define como
1 −1
 
r2 r1 r2 r2
U2 = ln + ln + (3.4.10)
k1 r0 k2 r1 h1
 3.5 Conducción de calor en coordenadas esféricas 57

Ejercicio 3.10 Determine la tasa de transferencia de calor basada en la superficie interna del
tubo, es decir, Q̇ r1
así como el coeficiente global de transferencia de calor para la superficie
interna. 

Ejercicio 3.11 — Resistencias internas a la transferencia de calor. Vuelva a resolver el

problema para el caso en que la superficie en r = r0 ofrezca resistencias a la transferencia de
calor, por lo que la condición de frontera dada en la ecuación (3.4.5c) se reemplaza por

dhT1 iθ z  
En r = r0 , k1 = h0 hT1 iθ z − T0 (3.4.11)
dz
Muestre que en este caso el coeficiente global de transferencia de calor es
 −1
r2 r2 r1 r2 r2 1
U2 = + ln + ln + (3.4.12)
r0 h0 k1 r0 k2 r1 h1


Ejercicio 3.12 — Nieve. Resuelva de nuevo el problema para el caso en el que el tubo aislante
2 esté cubierto de nieve como se muestra en el esquema de la figura.

aire
r2 − r1

r1 − r0

Nieve a Tn r0

Aislante 2 Aislante 1 Tubo a T0

3.5 Conducción de calor en coordenadas esféricas

Considere la transferencia de calor a través de una coraza esférica, la cual en el radio r0 se

mantiene a una temperatura T1 y en el radio r1 está expuesta al ambiente. Se desea determinar una
expresión para el perfil de temperatura así como del flux de calor considerando que la coraza ofrece
resistencia a la transferencia de calor en su superficie externa.
En este caso, la ecuación diferencial es:

∂ 2T
   
1 ∂ 2∂T 1 ∂ ∂T 1
0= r + senθ + (3.5.1a)
r2 ∂ r ∂r r2 senθ ∂ θ ∂θ r2 sen2 θ ∂ φ 2
58 Capítulo 3. Conducción de calor en estado estacionario

La cual está sujeta a las siguientes condiciones de frontera:

En r = r0 , T = T0 (3.5.1b)
∂T
En r = r1 , −k = h(T − Ta ) (3.5.1c)
∂r
Suponiendo que la transferencia de calor se de forma uniforme a través de la coraza, se pueden
imponer las siguientes condiciones de periodicidad:
∂T ∂T
= (3.5.1d)
∂ θ θ =0 ∂ θ θ =π
∂T ∂T
= (3.5.1e)
∂ φ φ =0 ∂ φ φ =2π

Para eliminar la dependencia de la temperatura con las coordenadas angular y azimutal, se define el
siguiente operador de promediado:
φZ=2π θZ=π
θφ 1
hT i = T senθ dθ dφ (3.5.2)
4π
φ =0 θ =0

Aplicando este operador a la ecuación diferencial considerando las condiciones de periodicidad se

obtiene:

2 dhT i
 θφ

1 d
0= 2 r (3.5.3a)
r dr dr
La cual está sujeta a las condiciones de frontera:

En r = r0 , hT iθ φ = T0 (3.5.3b)
dhT iθ φ  
En r = r1 , −k = h hT iθ φ − Ta (3.5.3c)
dr
Antes de resolver el problema, es conveniente escribirlo en términos de las siguientes variables y
parámetros adimensionales:
hT iθ φ − Ta r r0 hr1
T∗ = ; r∗ = ; r0∗ = ; Bi = (3.5.4)
T0 − Ta r1 r1 k
De esta forma, el problema adimensional es:
∗
 
1 d ∗2 dT
0 = ∗2 ∗ r (3.5.5a)
r dr dr∗

En r∗ = r0∗ , T∗ = 1 (3.5.5b)
dT ∗
En r∗ = 1, − ∗ = BiT ∗ (3.5.5c)
dr
Claramente, el problema está dirigido por la única fuente que se localiza en la frontera ξ = ξ0 . La
solución de este problema es
 
1
1 + Bi ∗ − 1
∗ r
T =   (3.5.6)
1
1 + Bi ∗ − 1
r0
 3.6 Una introducción al escalamiento: La aleta de enfriamiento 59

Figura 3.4: Ejemplos de aletas de enfriamiento rectangulares y cilíndricas.

Con este resultado puede calcularse el flux de calor como sigue

dhT iθ φ k (T0 − T1 ) dT ∗ Bik (T0 − T1 )

hqr iθ φ = −k =− ∗
= h  i (3.5.7)
dr r1 dr r2∗ r1 1 + Bi r1∗ − 1
0

Y de esta forma resulta que la tasa de transferencia de calor en la superficie externa de la esfera es
4πr0 Bik (T0 − T1 )
Q̇ = 4πr12 hqr iθ φ = (3.5.8)
1 + Bi 1 − r0∗
r=r1


r=r1

Es interesante tomar los valores extremos del número de Biot para este problema. Para el caso
en el que la coraza ofrezca demasiadas resistencias a la transferencia de calor (Bi → 0), la ecuación
anterior se reduce a:

Bi → 0, T∗ = 1 (3.5.9)

Es decir, la temperatura de la coraza se mantendrá a T0 y no disipará el calor. En el otro extremo, si

la coraza no ofrece resistencias a la transferencia de calor Bi  1, se obtiene que:
1
∗
−1
Bi  1, ∗
T = r (3.5.10)
1
−1
r0∗
El cual es el perfil de temperatura para el caso en el que en r = r1 , T = Ta (o bien, en r∗ = 1,
T ∗ = 0).
Por último, el flux de calor es en este caso
h (T0 − Ta ) 1
hqr iθ φ =   ∗2 (3.5.11)
1
1 + Bi r∗ − 1 r
0

La cual es una función parabólica decreciente con la posición.

3.6 Una introducción al escalamiento: La aleta de enfriamiento

Cuando una superficie posee una limitada tasa de disipación de calor al ambiente, se suelen
agregar aletas de enfriamiento para mejorar la tasa global de transferencia de calor. Algunos
ejemplos de aletas de enfriamiento rectangulares y cilíndricas se muestran en la Figura 3.4
Considere la transferencia de calor en estado estacionario a través de una aleta de enfriamiento
rectangular de largo, L (dirección x), anchura b (dirección z) y profundidad W (dirección y) (ver
figura 3.5). La ecuación de energía térmica bajo estas condiciones es
∂ 2T ∂ 2T ∂ 2T
+ 2 + 2 =0 (3.6.1)
∂ x2 ∂y ∂z
60 Capítulo 3. Conducción de calor en estado estacionario

Las condiciones de frontera correspondientes son

En x = 0, T = Tp (3.6.2a)
∂T
En x = L, −k = h(T − Ta ) (3.6.2b)
∂x
∂T
En y = 0, k = h(T − Ta ) (3.6.2c)
∂y
∂T
En y = W , −k = h(T − Ta ) (3.6.2d)
∂y
∂T
En z = −b/2, k = h(T − Ta ) (3.6.2e)
∂z
∂T
En z = b/2, −k = h(T − Ta ) (3.6.2f)
∂z

donde Tp y Ta son los valores de la temperatura en contacto con la superficie caliente y lejos de la
superficie de la aleta, respectivamente; mientras que h es el coeficiente interfacial de transferencia
de calor. Para reducir la dimensión del problema, se define el siguiente operador de promediado

z=b/2 y=W
1
Z Z
yz
hT i = T dydz (3.6.3)
bW
z=−b/2 y=0

Aplicando este operador a la ecuación diferencial resulta que

 z=b/2

y=W y=W z=b/2
d 2 hT iyz 1  ∂T ∂T
Z Z
+ dz + dy = 0 (3.6.4)
dx2 bW ∂y y=0 ∂z z=−b/2
z=−b/2 y=0

Figura 3.5: Esquema de una pared con aletas de enfriamiento con sus dimensiones y ubicación del
sistema coordenado para estudiar una sola aleta.
 3.6 Una introducción al escalamiento: La aleta de enfriamiento 61

O bien, tras aplicar las condiciones de frontera,

        
d 2 hT iyz h T |y=W − Ta + h T |y=0 − Ta h T |z=b/2 − Ta + h T |z=−b/2 − Ta
k − + =0
dx2 W b
(3.6.5)

Para simplificar esta ecuación, se supone que la temperatura a lo largo del perímetro de la sección
transversal de la aleta es la misma; esto es

T |y=W = T |y=0 = T |z=b/2 = T |z=−b/2 = Ts (3.6.6)

Bajo esta suposición, la ecuación (3.6.5) se reduce a:

d 2 hT iyz hP
− (Ts − Ta ) = 0 (3.6.7)
dx2 kAc

donde P = 2(b + W ) y Ac = bW son el perímetro y el área de la sección transversal de la aleta,

respectivamente.

3.6.1 Análisis basado en suponer que la temperatura en la superficie es igual a la

temperatura promedio
Hasta este punto, los desarrollos están basados en suposiciones razonables y fáciles de cumplir
en la práctica. Sin embargo, en la literatura se suele suponer que la temperatura en la superficie es
igual a la temperatura promedio, esto es

Ts ≈ hT iyz (3.6.8)

Bajo esta suposición, la ecuación (3.6.7) se puede escribir como

d 2 hT iyz hP
− (hT iyz − Ta ) = 0 (3.6.9a)
dx2 kAc
La cual está sujeta a las siguientes dos condiciones de frontera

En x = 0, hT iyz = Tp (3.6.9b)
dhT iyz
En x = L, −k = h (hT iyz − Ta ) (3.6.9c)
dx
Antes de resolver el problema, es conveniente reformularlo en términos de las siguientes variables
y números adimensionales

hT iyz − Ta x hL BiPL 2(b +W )hL2

T∗ = , x∗ = ; Bi = ; m2 = = (3.6.10)
Tp − Ta L k Ac bW k

De esta forma, el problema adimensional es

d2T ∗
− m2 T ∗ = 0 (3.6.11a)
dx∗2
En x∗ = 0, T ∗ = 1 (3.6.11b)
dT ∗
En x∗ = 1, − ∗ = BiT ∗ (3.6.11c)
dx
62 Capítulo 3. Conducción de calor en estado estacionario

La solución de este problema es

Bisenh [m (1 − x∗ )] + mcosh [m (1 − x∗ )]
T∗ = (3.6.12)
Bisenh (m) + m cosh (m)

Para obtener este resultado se usaron las siguientes identidades:

senh [m (1 − x∗ )] = senh (m) cosh (mx∗ ) − cosh (m) senh (mx∗ )

(3.6.13)
cosh [m (1 − x∗ )] = cosh (m) cosh (mx∗ ) − senh (m) senh (mx∗ )

Con este resultado se puede calcular la tasa de transferencia de calor que recibe la aleta como:

∂T
Q̇ f = W bqx = −kW b (3.6.14)
∂x x=0

O bien, promediando

Bi
dhT iyz Wb dT ∗ mkW b + tanh (m)
hQ̇ f iyz = −kW b = −k (Tp − Ta ) ∗ = (Tp − Ta ) m
dx L dx x∗ =0 L Bi
x=0 tanh (m) + 1
m
(3.6.15)
√
Note que mkWL
b
= PhkAc . Note que este resultado se reduce a la siguiente expresión cuando
Bi → 0, esto es, cuando la aleta está aislada:
p
Bi → 0, hQ̇ f iyz = PhkAc (Tp − Ta ) tanh (m) (3.6.16)

La eficiencia de una aleta se define como la razón entre la tasa de transferencia de calor y la máxima
la tasa de transferencia de calor que puede disipar la aleta, la cual corresponde al caso en el que
toda la aleta se encuentre a la temperatura Tp , de esta forma

hQ̇ f iyz
máx = PLh(Tp − Ta ) (3.6.17)

De esta forma, la eficiencia de la aleta es:

Bi
m2
+ tanh(m)
m
ηf = Bi (3.6.18)
m tanh (m) +1

Figura 3.6: Esquema de una superficie con varias aletas de enfriamiento.

 3.6 Una introducción al escalamiento: La aleta de enfriamiento 63

Más aún, para una placa que tenga n aletas de enfriamiento de espesor b, separadas a una
distancia d − b (ver figura 3.6), se tiene lo siguiente:

A f = nbW ; Au = n(d − b)W ; A = ndW (3.6.19)

Por lo que la eficiencia global de la aleta está dada por

Bi tanh (m)
Au + A f η f A − A f (1 − η f ) Af b b m2 + m
η= = = 1− (1 − η f ) = 1 − + (3.6.20)
A A A d d Bi
tanh (m) + 1
m
Aquí se tomó en cuenta que Au = A − A f . Para el caso en el que Bi → 0, el resultado anterior se
reduce a
b b tanh(m)
η = 1− + (3.6.21)
d d m

Ejercicio 3.13 Resuelva de nuevo el problema para el caso en el que la condición de frontera
en x = L sea

En x = L, T = TL (3.6.22)

Muestre que en este caso la tasa de transferencia de calor es

p cosh(m) − TL∗
hQ̇ f iyz = PhkAc (Tp − Ta ) (3.6.23)
senh(m)


Ejercicio 3.14 Algunos autores han sugerido que una aproximación razonable de la ecuación
(3.6.20) es considerar el caso para Bi → 0 utilizando la longitud corregida Lc∗ = 1 + b/(2L), de
tal manera que la eficiencia global de la aleta se aproxima como

b b tanh(mLc∗ )
η = 1− + (3.6.24)
d d mLc∗

Compare las predicciones de esta ecuación con las de la ecuación (3.6.20) y discuta la utilidad
de esta aproximación. 

3.6.2 Análisis basado en en la suposición de desviaciones periódicas

Regresando la atención a la ecuación diferencial (3.6.5), en lugar de suponer que las temperatu-
ras en las superficies son todas iguales a una temperatura superficial y que ésta a su vez es igual a la
temperatura promedio, se puede proponer la siguiente descomposición

T = hT iyz + T̃ (3.6.25)

donde T̃ se refiere a las desviaciones espaciales de la temperatura. Usando esta descomposición, la

ecuación (3.6.5) se escribe como
    
d 2 hT iyz 2h (hT i yz − Ta ) + h T̃ y=W
+ T̃ y=0
2h (hT i yz − Ta ) + h T̃ z=b/2
+ T̃ z=−b/2
k − + =0
dx2 W b
64 Capítulo 3. Conducción de calor en estado estacionario

(3.6.26)

O bien, tras hacer algunos pasos algebraicos

!
d 2 hT iyz hP h T̃ y=W
+ T̃ y=0
T̃ z=b/2
+ T̃ z=−b/2
− (hT i yz − Ta ) − + =0 (3.6.27)
dx2 kAc k W b

Note que esta ecuación corresponde a la ecuación (3.6.9a) más un nuevo término referente a las
desviaciones de la temperatura. Para calcular estas desviaciones se resta la ecuación anterior a la
ecuación que gobierna a T , es decir, a la ecuación (3.6.1), el resultado es
!
∂ 2 T̃ ∂ 2 T̃ ∂ 2 T̃ hP h T̃ y=W
+ T̃ y=0
T̃ z=b/2
+ T̃ z=−b/2
0= 2
+ 2 + 2 + (hT i yz − Ta ) + +
∂x ∂y ∂z kAc k W b
(3.6.28)

Suponiendo que

T̃ y=W
= T̃ y=0
= T̃ z=b/2
= T̃ z=−b/2
= T̃s (3.6.29)

y tomando como dominio de solución a la sección transversal de la aleta (es decir al plano yz), la
ecuación (3.6.28) se reduce a

∂ 2 T̃ ∂ 2 T̃ hP yz


0= + + hT i − Ta + T̃s (3.6.30a)
∂ y2 ∂ z2 kAc
La cual está sujeta a las siguientes condiciones de frontera

∂ T̃
= h hT iyz − Ta + T̃s


En y = 0, k (3.6.30b)
∂y
∂ T̃
= h hT iyz − Ta + T̃s


En y = W , −k (3.6.30c)
∂y
∂ T̃
= h hT iyz − Ta + T̃s


En z = −b/2, k (3.6.30d)
∂z
∂ T̃
= h hT iyz − Ta + T̃s


En z = b/2, −k (3.6.30e)
∂z
La solución formal de este problema es

T̃ = b hT i yz − Ta + T̃s


(3.6.31)

Donde la variable de cerradura b resuelve el siguiente problema de valor a la frontera

∂ 2 b ∂ 2 b hP
0= + + (3.6.32a)
∂ y2 ∂ z2 kAc
∂b
En y = 0, k =h (3.6.32b)
∂y
∂b
En y = W , −k =h (3.6.32c)
∂y
∂b
En z = −b/2, k =h (3.6.32d)
∂z
∂b
En z = b/2, −k =h (3.6.32e)
∂z
 3.7 Aletas de enfriamiento cilíndricas 65

Además, la variable de cerradura debe satisfacer la siguiente restricción promedio

hbiyz = 0 (3.6.32f)

Utilizando la ecuación (3.6.31), se puede llegar a

bs
T̃s = (hT i yz − Ta ) (3.6.33)
1 − bs
Con esta información, se deduce que el modelo promedio cerrado es

d 2 hT iyz hP
2
− (hT iyz − Ta ) (3.6.34)
dx kAc (1 − bs )

Claramente, cuando bs = 0 este resultado se reduce al dado en la ecuación (3.6.9a). Note que esta
modificación al modelo promedio no altera los resultados mostrados arriba; la única modificación
es que
BiPL
m2 = (3.6.35)
Ac (1 − bs )

3.7 Aletas de enfriamiento cilíndricas

Considere ahora que en lugar de usar una aleta rectangular se utiliza una aleta cilíndrica de
radio R y largo L. La ecuación que gobierna la transferencia de calor en cada punto de la aleta es

1 ∂ 2T ∂ 2T
 
1 ∂ ∂T
0= r + 2 + (3.7.1a)
r ∂r ∂r r ∂ θ 2 ∂ z2

La cual está sujeta a las siguientes condiciones de frontera

∂T
En r = 0, =0 (3.7.1b)
∂r
∂T
En r = R, −k = h(T − Ta ) (3.7.1c)
∂r
En z = 0, T = T0 (3.7.1d)
∂T
En z = L, −k = h(T − Ta ) (3.7.1e)
∂z
En este caso se define el siguiente operador de promediado
θZ=2π r=R
1
Z
θr
hT i = Trdrdθ (3.7.2)
πR2
θ =0 r=0

Considerando condiciones de periodicidad en la dirección angular, el resultado de aplicar este

operador a la ecuación de energía térmica es

d 2 hT iθ r 2h
0= − ( T |r=R − Ta ) (3.7.3)
dz2 Rk
Para simplificar esta ecuación se supondrá que

T |r=R ≈ hT iθ r (3.7.4)
 66 Capítulo 3. Conducción de calor en estado estacionario

Figura 3.7: Esquema de una aleta anular de espesor constante.

De esta forma, la ecuación (3.7.3) se escribe como

d 2 hT iθ r 2h  θ r 
0= − hT i − Ta (3.7.5a)
dz2 Rk
La cual está sujeta a las condiciones de frontera

En z = 0, hT iθ r = T0 (3.7.5b)
dhT iθ r  
En z = L, −k = h hT iθ r − Ta (3.7.5c)
dz
Para resolver el problema es conveniente expresarlo en términos de las siguientes variables y
parámetros adimensionales

hT iθ r − Ta z hL 2BiL
T∗ = ; z∗ = ; Bi = ; m2 = (3.7.6)
T0 − Ta L k R
De esta forma, el problema adimensional es

d2T ∗
0= − m2 T ∗ (3.7.7a)
dz∗2
En z∗ = 0, T ∗ = 1 (3.7.7b)
dT ∗
En z∗ = 1, − ∗ = BiT ∗ (3.7.7c)
dz
Este problema es idéntico al dado en las ecuaciones (3.6.11) por lo que pueden seguirse todos los
desarrollos mostrados arriba para obtener el perfil de temperatura promedio en la aleta.

3.8 Aletas de enfriamiento anulares

Considere la aleta de enfriamiento esquematizada en la figura 3.7. Bajo condiciones de estado
estacionario, la ecuación de energía térmica se reduce a

1 ∂ 2T ∂ 2T
 
1 ∂ ∂T
0= r + 2 + (3.8.1a)
r ∂r ∂r r ∂ θ 2 ∂ z2
3.8 Aletas de enfriamiento anulares 67

La cual está sujeta a las siguientes condiciones de frontera

En r = r0 , T = T0 (3.8.1b)
∂T
En r = r1 , −k = h (T − Ta ) (3.8.1c)
∂r
∂T
En z = −b/2, k = h (T − Ta ) (3.8.1d)
∂z
∂T
En z = b/2, −k = h (T − Ta ) (3.8.1e)
∂z
Para simplificar esta ecuación se define el siguiente operador de promediado
θZ=2π z=b/2
1
Z
θz
hT i = T dzdθ (3.8.2)
2πb
θ =0 z=−b/2

Tomando en cuenta la periodicidad en la dirección angular, el resultado de aplicar el operador de

promediado a la ecuación diferencial es

dhT iθ z
 
1 d 2h  θ z 
0= r − hT i − Ta (3.8.3)
r dr dr bk
Para llegar a este resultado se supuso que

T |z=±b/2 ≈ hT iθ z (3.8.4)

La ecuación diferencial promedio está sujeta a las siguientes condiciones de frontera

En r = r0 , hT iθ z = T0 (3.8.5)
dhT iθ z  
En r = r1 , −k = h hT iθ z − Ta (3.8.6)
dr
Para resolver este problema es conveniente reformularlo en términos de las siguientes variables y
parámetros adimensionales

hT iθ z − Ta r r0 hr1 Bi2r1
T∗ = ; r∗ = ; r0∗ = ; Bi = ; m2 = (3.8.7)
T0 − Ta r1 r1 k b
De esta forma, el problema adimensional es
∗
 
1 d ∗ dT
0= ∗ ∗ r − m2 T ∗ (3.8.8a)
r dr dr∗
En r∗ = r0∗ , T ∗ = 1 (3.8.8b)
dT ∗
En r∗ = 1, − ∗ = BiT ∗ (3.8.8c)
dr
Al desarrollar la derivada en la dirección radial en la ecuación (3.8.8a) y multiplicando la
ecuación resultante por r∗2 se obtiene que

d2T ∗ ∗ dT
∗
0 = r∗2 + r − m2 r∗2 T ∗ (3.8.9)
dr∗2 dr∗
Esta expresión tiene una forma compatible con la ecuación diferencial de Bessel modificada

d2y dy
x2 2
+ x − (αx2 − n2 )y = 0, n≥0 (3.8.10)
dx dx
68 Capítulo 3. Conducción de calor en estado estacionario

a) b)

Figura 3.8: Evaluación de las funciones de Bessel modificadas a) de primera y b) segunda clase.

Las soluciones generales de esta ecuación son

√ √
y = AIn ( αx) + BI−n ( αx), ∀n 6= 0, 1, 2, . . . (3.8.11a)
√ √
y = AIn ( αx) + BKn ( αx), ∀n = 0, 1, 2, . . . (3.8.11b)

donde In (x) y Kn (x) son las funciones de Bessel modificadas de primera y segunda clase, respecti-
vamente. La evaluación de estas funciones se muestra en la figura 3.8. Como puede notarse, La
función Bessel modificada de primera clase es una función creciente con la posición, mientras
que lo opuesto ocurre con la función modificada de segunda clase. Estas funciones satisfacen las
siguientes identidades
√
dI0 ( αx) √ √
= αI1 ( αx) (3.8.12a)
dx
√
dK0 ( αx) √ √
= − αK1 ( αx) (3.8.12b)
dx
Comparando las ecuaciones (3.8.9) y (3.8.10), se deduce que n = 0 y α = m2 , por lo tanto, la
solución general de la ecuación diferencial (3.8.9) es

T ∗ = c1 I0 (mr∗ ) + c2 K0 (mr∗ ) (3.8.13)

Usando la condición de frontera en r∗ = r0∗ se obtiene que

1 − c1 I0 (mr0∗ )
c2 = (3.8.14)
K0 mr0∗

Sustituyendo este resultado en la ecuación (3.8.13) se obtiene que

" #
∗ ∗ I0 (mr∗ ) K0 (mr∗ ) K0 (mr∗ )
T = c1 I0 (mr0 ) − + (3.8.15)
I0 mr0∗ K0 mr0∗ K0 mr0∗

Usando ahora la condición de frontera en r∗ = 1 se determina a la constante de integración faltante

como
mK (m) − BiK0 (m)
c1 =
1   (3.8.16)
Bi I0 (m) K0 mr0 − I0 mr0 K0 (m) + m I1 (m) K0 mr0∗ + I0 mr0∗ K1 (m)
∗ ∗




3.8 Aletas de enfriamiento anulares 69

De esta forma, la solución del problema se puede expresar como

Bi [I0 (m) K0 (mr∗ ) − I0 (mr∗ ) K0 (m)] + m [I1 (m) K0 (mr∗ ) + I0 (mr∗ ) K1 (m)]
T∗ =  (3.8.17)
Bi I0 (m) K0 mr0∗ − I0 mr0∗ K0 (m) + m I1 (m) K0 mr0∗ + I0 mr0∗ K1 (m)


 

Con este resultado puede calcularse la tasa de transferencia de calor que recibe la aleta como

∂T
Q̇ r=r0
= 2πr0 b qr |r=r0 = −2πr0 bk (3.8.18)
∂r r=r0

Promediando ambos lados de la ecuación anterior y recordando las definiciones adimensionales de

la posición radial y la temperatura se obtiene

dhT iθ z 2πr0 bk dT ∗
hQ̇iθ z r=r0
= −2πr0 bk =− (T0 − Ta ) ∗ (3.8.19)
dr r=r0 r1 dr r∗ =r0∗

Tomando la derivada de T ∗ respecto a r∗ en la ecuación (3.8.17), evaluando el resultado en r∗ = r0∗

y sustituyendo la expresión obtenida en la expresión anterior da lugar a

θz 2πr0 bk
hQ̇i = (T0 − Ta ) m
r=r0 r1
Bi [I0 (m) K1 (mr0∗ ) + I1 (mr0∗ ) K0 (m)] + m [I1 (m) K1 (mr0∗ ) − I1 (mr0∗ ) K1 (m)]
(3.8.20)
Bi I0 (m) K0 mr0∗ − I0 mr0∗ K0 (m) + m I1 (m) K0 mr0∗ + I0 mr0∗ K1 (m)


 



La tasa máxima de transferencia de calor se calcula como

θz
hQ̇imáx = A f h (T0 − Ta ) (3.8.21)

Donde el área de la aleta es

"  #
b 2 b 2
 
A f = 2π r12 − r0 2


+ 2πr1 b + 2πr0 b = 2π r1 + − r0 − (3.8.22)
2 2

De esta forma, puede por último calcularse la eficiencia de la aleta como sigue

θz
hQ̇i
r=r0 2r0∗
ηf = z
=  2  2 
hQ̇iθmáx b ∗ b
m 1 + 2r1 − r0 − 2r1

Bi [I0 (m) K1 (mr0∗ ) + I1 (mr0∗ ) K0 (m)] + m [I1 (m) K1 (mr0∗ ) − I1 (mr0∗ ) K1 (m)]
(3.8.23)
Bi I0 (m) K0 mr0∗ − I0 mr0∗ K0 (m) + m I1 (m) K0 mr0∗ + I0 mr0∗ K1 (m)


 


 4. Transferencia de calor en más de una dirección

4.1 El problema de Sturm-Liouville

Para los desarrollos que siguen es de especial importancia el siguiente tipo de ecuación diferen-
cial ordinaria
 
d dφ
p(x) + q(x)φ + λ 2 σ (x)φ = 0 (4.1.1)
dx dx

donde el dominio es x ∈ (a, b), λ es una constante conocida como valor propio, mientras que
p(x), q(x) y σ (x) son funciones conocidas de la posición. Esta ecuación diferencial está sujeta a
condiciones de frontera homogéneas y posee las siguientes propiedades (Haberman, 2012):
1. Todos los valores propios son reales.
2. Existe un número infinito de valores propios que puede ordenarse a partir del valor más
pequeño:

λ1 < λ2 < λ3 < . . . (4.1.2)

Donde λn → ∞ cuando n → ∞.
3. A cada valor propio λn le corresponde una función propia φn (x), la cual tiene n − 1 ceros en
el intervalo a < x < b.
4. Las funciones propias forman un conjunto completo que permite expresar a cualquier función
suave por tramos, f (x), como una combinación lineal de ellas:
∞
f (x) = ∑ An φn (x) (4.1.3)
n=1

5. Dos funciones propias (que correspondan a diferentes valores propios) son ortogonales
relativas a la función peso σ (x), esto es
x=b
Z
φn (x)φm (x)σ (x)dx = 0, ∀n 6= m (4.1.4)
x=a
 72 Capítulo 4. Transferencia de calor en más de una dirección

Como se verá más adelante, estas propiedades son cruciales en la solución de ecuaciones dife-
renciales parciales. No debe olvidarse que un problema de Sturm-Liouville se conforma por una
ecuación diferencial homogénea sujeta a condiciones de frontera que también son homogéneas.
Para problemas multidimensionales, donde la variable dependiente f es, por ejemplo, una función
de x, y y z lo que se hace es proponer una solución de este tipo

f = ∑∑∑ Aφ (x)ϕ(y)ψ(z) (4.1.5)

Si el problema tiene condiciones de frontera homogéneas, entonces las funciones propias se

obtienen en esas direcciones. De no ser este el caso, se pueden establecer de manera arbitraria en
qué direcciones se calcularán las funciones propias. De hecho, en principio es posible calcular las
funciones propias en todas las direcciones. Una vez determinadas las funciones propias, puede
usarse la fórmula de Green (ver ecuación 1.4.7):
Z Z
Ψ∇2 f − f ∇2 Ψ dV = n · (Ψ∇ f − f ∇Ψ) dA


3D : (4.1.6a)
V A
Z Z
2 2


2D : Ψ∇ f − f ∇ Ψ dA = n · (Ψ∇ f − f ∇Ψ) dσ (4.1.6b)
A σ
Z
Ψ∇2 f − f ∇2 Ψ dσ = n · (Ψ∇ f − f ∇Ψ)|σσsup


1D : in f
(4.1.6c)
σ

para determinar el valor de f . En las ecuaciones anteriores Ψ representa a la función (1D) o

funciones (2D y 3D) propias. Por ejemplo, en el caso tridimensional, si las funciones propias aplican
a las tres direcciones, la ecuación (4.1.6a) conduce a una ecuación algebraica para determinar al
coeficiente A en la ecuación (4.1.5). Sin embargo, si las funciones propias sólo se determinaron
en dos direcciones, entonces la ecuación (4.1.6a) conduce a una ecuación diferencial ordinaria. El
uso de esta técnica basada en funciones propias y la fórmula de Green se explora en las siguientes
secciones.

4.2 Coordenadas rectangulares

4.2.1 Transporte bidimensional
Considere una viga de sección transversal rectangular de base b, alto L y profundidad W , como
las vías que se encuentran en las líneas del metro o del tranvía. Suponga que todas las caras de la
viga se encuentran a una temperatura uniforme T0 y que únicamente la cara superior está a una
temperatura T1 > T0 . Plantee y resuelva el modelo de transferencia de calor en estado estacionario
considerando que no hay una fuente interna de transferencia de calor. Suponiendo que b, L  W , la
ecuación de energía térmica se reduce a:
∂ 2T ∂ 2T
0= + 2 (4.2.1a)
∂ x2 ∂y
La cual está sujeta a las siguientes cuatro condiciones de frontera,

En x = 0, b, T = T0 (4.2.1b)
En y = 0, T = T0 (4.2.1c)
En y = L, T = T1 (4.2.1d)

Para resolver este problema es conveniente adoptar los siguientes cambios de variables
T − T0 x y b
T∗ = ; x∗ = ; y∗ = ; b∗ = (4.2.2)
T1 − T0 L L L
4.2 Coordenadas rectangulares 73

Figura 4.1: Esquema de la sección transversal de una viga con las temperaturas especificadas en
cada cara.

De esta forma, el problema puede escribirse en su forma adimensional como:

∂ 2T ∗ ∂ 2T ∗
0= + ∗2 (4.2.3a)
∂ x∗2 ∂y
En x = 0, b , T ∗ = 0
∗ ∗
(4.2.3b)
∗ ∗
En y = 0, T =0 (4.2.3c)
∗ ∗
En y = 1, T =1 (4.2.3d)

Como puede notarse existen condiciones de frontera homogéneas en la dirección x, por lo que se
propone la siguiente solución
∞
T ∗ (x∗ , y∗ ) = ∑ An (y∗ )φn (x∗ ) (4.2.4)
n=1

Las funciones propias deben satisfacer un problema de Sturm-Liouville compatible con el problema
original, es decir, deben tener condiciones tipo Dirichlet. Las funciones propias deben resolver el
siguiente problema

d 2 φn
+ λn2 φn = 0 (4.2.5a)
dx∗2
En x∗ = 0, b∗ , φn = 0 (4.2.5b)

La solución general de esta ecuación es

φn = c1n sen (λn x∗ ) + c2n cos (λn x∗ ) (4.2.6)

Usando la condición de frontera en x∗ = 0 se deduce que c2n = 0 y usando la segunda condición de

frontera se obtiene que

0 = c1n sen (λn b∗ ) (4.2.7)

74 Capítulo 4. Transferencia de calor en más de una dirección

Si se concluyera que c1n = 0 entonces toda la solución del problema sería cero. Como este no es el
resultado correcto, se debe cumplir que

0 = sen (λn b∗ ) (4.2.8)

Esta ecuación se llama condición de valores propios pues permite encontrar los valores de λn , que
en este caso son
nπ
λn = , n = 1, 2, 3, . . . , ∞ (4.2.9)
b∗

Del análisis de estos párrafos se deduce que φn = sen (nπx∗ /b∗ ); note que no es necesaria la
constante ya que al multiplicar a An en la ecuación (4.2.4) sólo se obtiene otra función de y. Es
conveniente recordar que las funciones propias satisfacen la siguiente condición de ortogonalidad
∗ =b∗
xZ
φn φm dx∗ = 0, ∀n 6= m (4.2.10)
x∗ =0

Ahora que ya se conocen las funciones y valores propios puede usarse la fórmula de Green en
1D para obtener lo siguiente
∗ =b∗
xZ x∗ =b∗
∂ 2T ∗ d 2 φm ∂T∗
 
dφm
φm ∗2 − T ∗ ∗2 dx = φm ∗ − T ∗ ∗
∗
(4.2.11)
∂x dx ∂x dx x∗ =0
x∗ =0

En donde se usa el subíndice m para denotar a un valor particular de las funciones propias.
Sustituyendo las ecuaciones diferenciales y condiciones de frontera para T ∗ y φn en la ecuación
anterior resulta que
∗ =b∗
xZ
∂ 2T ∗
 
φm − T ∗ λm2 dx∗ = 0 (4.2.12)
∂ y∗2
x∗ =0

Sustituyendo la ecuación (4.2.4) en el resultado anterior se obtiene

 xZ∗ =b∗
d 2 An
∞ 
∑ − λm2 An φn φm dx∗ = 0 (4.2.13)
n=1 dy∗2
x∗ =0

Con base en la fórmula de ortogonalidad dada en la ecuación (4.2.10), resulta que todos los términos
de la sumatoria son cero excepto cuando n = m, esto es,

 xZ∗ =b∗
d 2 An

− λn2 An φn2 dx∗ = 0 (4.2.14)
dy∗2
x∗ =0

x∗R=b∗
Ya que φn2 dx∗ 6= 0, se deduce que la ecuación anterior se satisface porque
x∗ =0

d 2 An
− λn2 An = 0 (4.2.15)
dy∗2
4.2 Coordenadas rectangulares 75

Esta ecuación requiere dos condiciones de frontera, las cuales se deducen al sustituir la solución
dada en la ecuación (4.2.4) en las ecuaciones (4.2.3c) y (4.2.3d):
∞
En y∗ = 0, ∑ An (0)φn (x∗ ) = 0 (4.2.16a)
n=1
∞
En y∗ = 1, ∑ An (1)φn (x∗ ) = 1 (4.2.16b)
n=1

Multiplicando ambos lados de las ecuaciones anteriores por φm e integrando desde x∗ = 0 hasta
x∗ = b∗ se obtienen las siguientes expresiones
∗ =b∗
xZ
∞
∗
En y = 0, ∑ An (0) φn φm dx∗ = 0 (4.2.17a)
n=1
x∗ =0
∗ =b∗
xZ ∗ =b∗
xZ
∞
En y∗ = 1, ∑ An (1) φn φm dx∗ = φm dx∗ (4.2.17b)
n=1
x∗ =0 x∗ =0

Tomando en cuenta la propiedad de ortogonalidad de las funciones propias (ecuación 4.2.10), se

pueden destruir las sumatorias de las ecuaciones anteriores para así obtener

En y∗ = 0, An (0) = 0 (4.2.18a)
x∗ =b∗ x∗ =b∗
nπx∗
φn dx∗ dx∗
R R

sen ∗
∗ x∗ =0 x∗ =0
b∗ b
− nπ [(−1)n − 1]
En y = 1, An (1) = x∗R=b∗
= x∗R=b∗
= b∗
∗
dx∗ sen2 nπx dx∗

2
φn2 b∗
x∗ =0 x∗ =0
2 h i
= 1 + (−1)n+1 (4.2.18b)
nπ
La solución general de la ecuación (4.2.14) es

An = d1n senh (λn y∗ ) + d2n cosh (λn y∗ ) (4.2.19)

O bien, después de usar las condiciones de frontera,

2 h i senh (λ y∗ )
n
An (y∗ ) = 1 + (−1)n+1 (4.2.20)
nπ senh (λn )
Sustituyendo este resultado en la ecuación (4.2.4) se obtiene la forma final de la solución
h i
n+1
∞ 1 + (−1)
2 senh (λn y∗ )
T ∗ (x∗ , y∗ ) = ∑ sen (λn x∗ ) (4.2.21)
π n=1 n senh (λn )

Como puede notarse, la mitad de la serie es nula; esto se debe a que siempre que n sea par, el
resultado de la resta entre corchetes es cero. Para atender este problema se propone el siguiente
cambio de variables:

n + 1 = 2k, n, k = 1, 2, . . . (4.2.22)

De esta forma, la solución es ahora

4 ∞ 1 senh (λk y∗ )
T ∗ (x∗ , y∗ ) = ∑ (2k − 1) senh (λk ) sen (λk x∗ ) (4.2.23)
π k=1
76 Capítulo 4. Transferencia de calor en más de una dirección

en donde
(2k − 1) π
λk = (4.2.24)
b∗
Existe aún un problema con esta solución y es que la función de seno hiperbólico se desborda
cuando su argumento es suficientemente grande. Para resolver esta situación, se propone utilizar la
identidad de Euler:
ex − e−x
senh(x) = (4.2.25)
2
Haciendo las operaciones algebraicas correspondientes, la solución toma la siguiente forma

4 ∞ exp [−λk (1 − y∗ )] 1 − exp (−2λk y∗ )

T ∗ (x∗ , y∗ ) = ∑ sen (λk x∗ ) (4.2.26)
π k=1 2k − 1 1 − exp (−2λk )

A continuación se lista el programa de Matlab para la evaluación de la solución en valores

específicos de y∗ :

Programa para la evaluación de la solución dada en la ecuación (4.2.26).

clear; close all; clc;

format long
b= 1;
kmax= 1000;

nx = round(100*b);
x=linspace(0,b, nx);
T=zeros(nx,1);
y= 0.7;

for i=1:nx
suma=0;
for k=1:kmax
lam= (2*k-1)*pi/b;
suma=suma+ 2/(2*k-1)*exp(-lam*(1-y))*(1-exp(-2*lam*y))/...
((1-exp(-2*lam)))*sin(lam*x(i));
end
T(i)= 2*suma/pi;
end
plot(x,T)

La evaluación de la solución se muestra en la figura 4.2 para el caso en que b∗ = 1. Como puede
notarse, los perfiles son simétricos alrededor de x∗ = 0.5 y se van incrementando conforme aumenta
y∗ . Estos resultados son los esperados de acuerdo a la física del problema, ya que la temperatura
más alta se encuentra en y∗ = 1.

Ejercicio 4.1 — Ortogonalidad. Sustituya φn = sen(nπx∗ /b∗ ) y φm = sen(mπx∗ /b∗ ) en la

ecuación (4.2.10) y lleve a cabo el paso de integración para verificar que el resultado es cero
sólo cuando n 6= m. 
4.2 Coordenadas rectangulares 77

1 1

y ∗ = 0.9
0.8 0.8

0.6 0.6
T∗ 0.7
0.4 0.4
0.5
0.2 0.2
0.3
0.1
0 0
0 0.2 0.4 0.6 0.8 1
x∗
Figura 4.2: Perfiles de temperatura adimensionales en la sección transversal de una viga. Se
utilizaron 100 términos en la sumatoria y se tomó b∗ = 1.

Ejercicio 4.2 — Influencia de b∗ . Explore la influencia que tiene el parámetro b∗ en los

perfiles de temperatura. 

Ejercicio 4.3 — Método de separación de variables. Este método sólo es aplicable a pro-
blemas que tengan una ecuación diferencial homogénea con condiciones de frontera homogéneas
en todas las direcciones excepto en una como es el caso del problema estudiado en esta sección.
En este método se propone que la solución sea:

T ∗ (x∗ , y∗ ) = f (x∗ )g(y∗ ) (4.2.27)

Sustituya esta propuesta en la ecuación (4.2.3a) y deduzca la siguiente expresión

1 d2 f 1 d2g
= − = −λ 2 (4.2.28)
f dx∗2 g dy∗2

La igualdad a −λ 2 se debe a que un lado de la igualdad sólo es función de x∗ , mientras que

el otro lado sólo es función de y∗ y la única forma de satisfacer esta condición es que ambos
lados sean iguales a una constante. El signo negativo se debe a que las condiciones de frontera
homogéneas están en la dirección x∗ por lo que es en esta dirección donde se busca obtener
un problema de tipo Sturm-Liouville. Deduzca y resuelva el problema para f y demuestre que
el resultado corresponde a φn . Resuelva el problema para g y, tras utilizar las condiciones de
frontera y la condición de ortogonalidad de las funciones propias, demuestre que se recupera la
solución dada en la ecuación. (4.2.21) 

Ejercicio 4.4 — Solución alternativa 1. En los párrafos anteriores se llevó a cabo la solución
proponiendo funciones propias en la dirección x para acabar obteniendo una función de y, sin
 78 Capítulo 4. Transferencia de calor en más de una dirección

embargo, la metodología también admite proponer una solución del tipo

∞
T ∗ (x∗ , y∗ ) = ∑ Bm (x∗ )ϕm (y∗ ) (4.2.29)
m=1

Resuelva el problema para este caso. 

Ejercicio 4.5 — Solución alternativa 2. Una alternativa más que puede considerarse es obte-
ner funciones propias en ambas direcciones con el fin de generar no una ecuación diferencial,
sino una ecuación algebráica para la incógnita. Para explorar esta alternativa, considere la
siguiente propuesta de solución
∞ ∞
T ∗ (x∗ , y∗ ) = ∑ ∑ Anm φn (x∗ )ϕm (y∗ ) (4.2.30)
n=1 m=1

Resuelva el problema para este caso. 

Ejercicio 4.6 — Temperatura promedio y flux de calor. Partiendo de la solución dada en la

ecuación (4.2.26), calcule la temperatura promediada en la dirección x∗ y grafique el resultado
como función de y∗ . Además, determine el flux de calor en la superficie y∗ = 1.


Ejercicio 4.7 — Fuente oscilatoria. Resuelva de nuevo el problema para el caso en el que la
temperatura en la superficie y∗ = 1 obedezca a la siguiente función oscilatoria
1
En y∗ = 1, T∗ = (1 + sen(8πx∗ )) (4.2.31)
2


4.2.2 Transporte por conducción en tres dimensiones

Considere la transferencia de calor en una habitación que tiene la forma de un prisma de sección
transversal rectangular. El techo de esta habitación cuenta con dos lámparas de las cuales sólo
una se encuentra encendida y aporta un flux de calor conocido. Además, del lado de la habitación
donde se encuentra encendida la lámpara está una ventana abierta y del lado opuesto se encuentra
la puerta de entrada. Por último, en el costado derecho está la entrada a la cocina y al calentador
de agua junto con un mueble con una televisión y del lado opuesto sólo hay estantes con libros.
Determine el perfil de temperatura en la habitación en estado estacionario. Desprecie los muebles y
los efectos convectivos.
4.2 Coordenadas rectangulares 79

El planteamiento del modelo matemático que describe la situación física descrita arriba es
∂ 2T ∂ 2T ∂ 2T
+ 2 + 2 =0 (4.2.32a)
∂ x2 ∂y ∂z
∂T
x = 0, k = h(T − Ta ) (4.2.32b)
∂x
∂T
x = L, =0 (4.2.32c)
∂x
∂T
y = 0, k = q0 (4.2.32d)
∂y
∂T
y=W =0 (4.2.32e)
∂y
∂T
z = 0, =0 (4.2.32f)
∂z

∂T q1 , x ∈ [0, L/2]
z = H −k = (4.2.32g)
∂z 0, x ∈ (L/2, L]
Donde q0 es el flux combinado de la televisión, la cocina y el calentador de agua, mientras que q1
es el flux de calor de la lámpara encendida. Antes de comenzar con la solución, es conveniente
plantear el problema en términos de las siguientes variables y parámetros adimensionales
T − Ta r hL qi L W H
T∗ = ; r∗ = ; Bi = ; q∗i = , i = 0, 1; W∗ = , H∗ = (4.2.33)
Ta L k kTa L L
De esta forma, el problema adimensional es
∂ 2T ∗ ∂ 2T ∗ ∂ 2T ∗
+ ∗2 + ∗2 = 0 (4.2.34a)
∂ x∗2 ∂y ∂z
∂ T ∗
x∗ = 0, = BiT ∗ (4.2.34b)
∂ x∗
∂T∗
x∗ = 1, =0 (4.2.34c)
∂ x∗
∂T∗
y∗ = 0, = q∗0 (4.2.34d)
∂ y∗
∂T∗
y∗ = W ∗ =0 (4.2.34e)
∂ y∗
∂T∗
z∗ = 0, =0 (4.2.34f)
∂ z∗
∂T∗
 ∗
q1 , x∗ ∈ [0, 1/2]
z∗ = H ∗ − ∗ = (4.2.34g)
∂z 0, x∗ ∈ (1/2, L]
Para resolver el problema se propone la siguiente solución
∞ ∞
T∗ = ∑ ∑ Anm (z∗ )φn (x∗ )ϕm (y∗ ) (4.2.35)
n=1 m=1

Donde las funciones propias satisfacen los siguientes problemas de Sturm-Liouville

d 2 φn
+ λn2 φn = 0 (4.2.36a)
dx∗2
dφn
x∗ = 0, = Biφn (4.2.36b)
dx∗
dφn
x∗ = 1, =0 (4.2.36c)
dx∗
80 Capítulo 4. Transferencia de calor en más de una dirección

d 2 ϕn
+ βm2 ϕm = 0 (4.2.37a)
dy∗2
dϕm
y∗ = 0, =0 (4.2.37b)
dy∗
dϕm
y∗ = W ∗ , =0 (4.2.37c)
dy∗

La solución general del primer problema es

φn = Asen(λn x∗ ) + Bcos(λn x∗ ) (4.2.38)

Aplicando la condición de frontera en x∗ = 0 se obtiene

Aλn = BiB (4.2.39)

Por lo que la solución es ahora

 
∗ Bi ∗ ∗
φn (x ) = B sen(λn x ) + cos(λn x ) (4.2.40)
λn

Usando ahora la condición de frontera en x∗ = 1 se obtiene la siguiente condición de valores propios

ya que B 6= 0

Bi
cos(λn ) = sen(λn ) (4.2.41)
λn

Dirigiendo la atención al problema para ϕm , su solución es

ϕm (y∗ ) = Ccos(βm y∗ ) (4.2.42)

en donde βm = (2m − 1)π/2W ∗ . De esta forma, la solución del problema puede expresarse como
∞ ∞
T∗ = ∑∑ fnm (z∗ )φn (x∗ )ϕm (y∗ ) (4.2.43)
n=1 m=1

donde fnm = Anm BC. Para determinar a la función fnm se utiliza la fórmula de Green en 2D:

xZ∗ =1 y∗Z=W ∗
∂ 2T ∗ ∂ 2T ∗ d2φp d 2 ϕq
   
∗
φ p ϕq + ∗2 − T ϕq ∗2 + φ p ∗2 dy∗ dx∗
∂ x∗2 ∂y dx dy
x∗ =0 y∗ =0
y∗Z=W ∗ x∗ =1 xZ∗ =1  y=W ∗
∂T∗ ∂T∗

∗ dφ p ∗ ∗ dϕq
= ϕq φ p ∗ − T dy + φ p ϕq ∗ − T dx∗ (4.2.44)
∂x dx∗ x∗ =0 ∂y dy∗ y=0
y∗ =0 x∗ =0

Sustituyendo las ecuaciones diferenciales y condiciones de frontera para T ∗ y las funciones propias
en la ecuación anterior se obtiene la siguiente expresión

xZ∗ =1 y∗Z=W ∗ xZ∗ =1

∂ 2T ∗


∗ ∗ ∗
2 2
− λ p + βq T φ p ϕq dy dx = ϕq (0)q0 φ p dx∗ (4.2.45)
∂ z∗2
x∗ =0 y∗ =0 x∗ =0
4.2 Coordenadas rectangulares 81

Sustituyendo la ecuación (4.2.43) en las integrales del lado izquierdo de la ecuación da como
resultado
 xZ∗ =1 y∗Z=W ∗ xZ∗ =1
d 2 fnm
∞ ∞ 
− λ p2 + βq2 fnm φn φ p dx∗ ϕm ϕq dy∗ = ϕq (0)q0 φ p dx∗


∑∑ dz ∗2
n=1 m=1
x∗ =0 y∗ =0 x∗ =0
(4.2.46)
Esta ecuación puede simplificarse considerablemente al considerar las fórmulas de ortogonalidad
de las funciones propias
xZ∗ =1
φn φ p dx∗ = 0, ∀n 6= p (4.2.47a)
x∗ =0
y∗Z=W ∗

ϕm ϕq dy∗ = 0, ∀m 6= q (4.2.47b)
y∗ =0

De esta forma las sumatorias en la ecuación (4.2.46) se reducen a un sólo término cada una y esta
ecuación se puede escribir como
x∗R=1
ϕm (0)q0 φn dx∗
d 2 fnm 2 2

x∗ =0
− λn + βm fnm = ∗ y∗ =W ∗ ≡ γnm (4.2.48)
dz∗2 x R=1
φn2 dx∗ ϕm2 dy∗
R
x∗ =0 y∗ =0

Note que la constante γnm contiene la influencia de la fuente q0 . Para deducir las condiciones de
frontera para esta ecuación diferencial se sustituye la ecuación (4.2.26) en las ecuaciones (4.2.34f)
y (4.2.34g)
∞ ∞
d fnm
z∗ = 0, ∑ ∑ φn (x∗ )ϕm (y∗ ) dz∗ =0 (4.2.49a)
n=1 m=1
q∗1 , x∗ ∈ [0, 1/2]
∞ ∞ 
∗ ∗ ∗ d fnm ∗
z =H − ∑ ∑ φn (x )ϕm (y ) ∗ = (4.2.49b)
n=1 m=1 dz 0, x∗ ∈ (1/2, L]

Tomando en cuenta las condiciones de ortogonalidad, se obtienen las siguientes condiciones de

frontera
d fnm
z∗ = 0, =0 (4.2.50a)
dz∗
y∗ =W ∗ x∗ =1/2
dy∗ φm dx∗
R R
ϕm
d fnm y∗ =0 x∗ =0
z∗ = H ∗ − ∗ = q∗1 y∗ =W ∗ ≡ Q∗1 (4.2.50b)
dz x∗R=1
ϕm2 dy∗ φn2 dx∗
R
y∗ =0 x∗ =0

La solución general de la ecuación (4.2.48) es

fnm = Asenh(ωnm z∗ ) + Bcosh(ωnm z∗ ) (4.2.51)

donde ωnm2 = λ 2 + β 2 . Usando la condición de frontera en z∗ = 0 se obtiene que A = 0, mientras

n m
que al usar la condición de frontera en z =∗ = H ∗ , el resultado es
Q∗1
B=− (4.2.52)
ωnm senh(ωnm H ∗ )
 82 Capítulo 4. Transferencia de calor en más de una dirección

De esta forma, la solución del problema es

∞ ∞
Q∗1 cosh(ωnm z∗ )
T∗ = − ∑ ∑ ωnm senh(ωnm H ∗ ) φn (x∗ )ϕm (y∗ ) (4.2.53)
n=1 m=1

4.3 Conducción de calor transitoria en una barra cilíndrica

En las horas pico en días calurosos las temperaturas dentro de los vagones del metro pueden
oscila entre los 30o C y 36o C. Las barras de acero que utilizan los pasajeros para sostenerse tienen
una difusividad térmica de 1.2 × 10−5 m2 /s, un radio de 5 cm y una longitud de 2 m. Si la barra
estaba inicialmente a una temperatura de 25o C, se desea examinar la dinámica de la transferencia
de calor tanto radial como axial en la barra. Considere que la conductividad térmica del acero es
de 16 W/mK y su coeficiente interfacial de transferencia de calor es de 50 W/m2 K. Suponga que
en los extremos de la barra la temperatura es T1 = 30o C, mientras que en su superficie lateral, la
temperatura es Ta = 36o C.
Para las condiciones que describe el problema, la ecuación de energía térmica y las condiciones
de frontera e iniciales correspondientes son
∂ 2T
 
∂T α ∂ ∂T
= r +α 2 (4.3.1a)
∂t r ∂r ∂r ∂z
∂T
r = 0, =0 (4.3.1b)
∂r
∂T
r = R, −k = h(T − Ta ) (4.3.1c)
∂r
∂T
z = 0, k = h(T − T1 ) (4.3.1d)
∂z
∂T
z = L, −k = h(T − T1 ) (4.3.1e)
∂z
t = 0, T = T0 (4.3.1f)
Note que ya se ha considerado transporte axisimétrico por lo que no se escribió la dependencia
angular de la temperatura. Antes de resolver el problema, es conveniente (aunque no obligatorio),
reformularlo en términos de las siguientes variables y parámetros adimensionales
T − T0 r z tα
T∗ = ; r ∗ = ; z∗ = ; t ∗ = 2 ;
Ta − T0 R R R
hR T1 − T0 L
Bi = = 0.16; T1∗ = = 0.45; L∗ = = 40 (4.3.2)
k Ta − T0 R
De esta forma, el problema adimensional es
∂T∗ ∗ ∂ 2T ∗
 
1 ∂ ∗∂T
= r + (4.3.3a)
∂t ∗ r∗ ∂ r∗ ∂ r∗ ∂ z∗2
∂T∗
r∗ = 0, =0 (4.3.3b)
∂ r∗
∂T∗
r∗ = 1, − ∗ = Bi(T ∗ − 1) (4.3.3c)
∂r
∂ T ∗
z∗ = 0, = Bi(T ∗ − T1∗ ) (4.3.3d)
∂ z∗
∂T∗
z∗ = L∗ , − ∗ = Bi(T ∗ − T1∗ ) (4.3.3e)
∂z
∗ ∗
t = 0, T = 0 (4.3.3f)
4.3 Conducción de calor transitoria en una barra cilíndrica 83

a) b)

Figura 4.3: Esquema de las funciones Bessel de a) primera y b) segunda clase.

A primera vista, pareciera que esta elección de variables adimensionales para la temperatura no
es adecuada pues no dio lugar a ninguna condición de frontera homogénea. Aunque esto es un
impedimento para usar el método de separación de variables, no lo es para el método usado en
párrafos anteriores. La solución del problema se propone como una doble expansión en términos de
las funciones propias en cada dirección, esto es,
∞ ∞
T ∗ (r∗ , z∗ ,t ∗ ) = ∑ ∑ Anm (t ∗ )φn (r∗ )ϕm (z∗ ) (4.3.4)
n=1 m=1

Donde Anm (t ∗ ) es una función del tiempo que se debe determinar junto con las funciones propias en
cada dirección. Para ello, se proponen los problemas de Sturm-Liouville para las funciones propias
en las direcciones radial y axial. A continuación se presenta el análisis en la dirección radial:
Problema en la dirección radial
 
1 d ∗ dφn
r + λn2 φn = 0 (4.3.5a)
r∗ dr∗ dr∗
dφn
r∗ = 0, =0 (4.3.5b)
dr∗
dφn
r∗ = 1, − ∗ = Biφn (4.3.5c)
dr
La ecuación diferencial corresponde a la ecuación de Bessel

d2y dy
x2 + x + ax2 − n2 y = 0,


2
n≥0 (4.3.6)
dx dx
La solución general de esta ecuación está dada por
√ √
y = AJn ( ax) + BJ−n ( ax), ∀n 6= 0, 1, 2, . . . (4.3.7a)
√ √
y = AJn ( ax) + BYn ( ax), ∀n = 0, 1, 2, . . . (4.3.7b)
√ √
En donde Jn ( ax) y Yn ( ax) son las funciones Bessel de primera y segunda clase, respectivamente
(ver figura 4.3). Estas funciones satisfacen las siguientes fórmulas de derivación
d  √
n √  √
n √ √
ax Jn ( ax) = ax Jn−1 ( ax) a (4.3.8a)
dx
d h √
−n √ i √
−n √ √
ax Jn ( ax) = − ax Jn+1 ( ax) a (4.3.8b)
dx
Para el caso de la ecuación (4.3.5a), la solución es

φn (r∗ ) = AJ0 (λn r∗ ) + BY0 (λn r∗ ) (4.3.9)

84 Capítulo 4. Transferencia de calor en más de una dirección

Como puede apreciarse de la figura 4.3b), Y0 (λn r∗ ) → −∞ cuando r∗ → 0, por ello, se decide que
B = 0 y la ecuación anterior se reduce a

φn (r∗ ) = AJ0 (λn r∗ ) (4.3.10)

Aplicando la condición de frontera en r∗ = 1 y tomando en cuenta que A 6= 0 resulta que

λn J1 (λn ) = BiJ0 (λn ) (4.3.11)

Esta ecuación es la condición de los valores propios. Con esta información se concluye el análisis
de la función propia en la dirección radial y se dirige la atención al análisis en la dirección axial.
Problema en la dirección axial
d 2 ϕm
+ βm2 ϕm = 0 (4.3.12a)
dz∗2
dϕm
z∗ = 0, = Biϕm (4.3.12b)
dz∗
dϕm
z∗ = L∗ , − ∗ = Biϕm (4.3.12c)
dz
La solución general de esta ecuación diferencial es

ϕm = Csen(βm z∗ ) + Dcos(βm z∗ ) (4.3.13)

Usando la primera condición de frontera se obtiene que

Cβm = BiD (4.3.14)

Por lo que la solución es ahora

 
Bi ∗ ∗
ϕm = D sen(βm z ) + cos(βm z ) (4.3.15)
βm

Usando ahora a la segunda condición de frontera, tomando en cuenta que D 6= 0, resulta lo siguiente

   
Bi Bi
Bi sen(βm ) + cos(βm ) + βm cos(βm ) − sen(βm ) = 0 (4.3.16)
βm βm

La cual es otra condición de valores propios. De esta forma, la solución del problema es
∞ ∞
T ∗ (r∗ , z∗ ,t ∗ ) = ∑∑ fnm (t ∗ )φn (r∗ )ϕm (z∗ ) (4.3.17)
n=1 m=1

Donde fnm (t ∗ ) = Anm (t ∗ )AD. Para determinar la ecuación diferencial que gobierna a esta función,
se utiliza la fórmula de Green en 2D:
z∗Z=L∗ rZ∗ =1
∗ ∂ 2T ∗ d 2 ϕq
       
∗ 1 ∂ ∗ ∗∂T ∗ ϕq d ∗ dφ p
φ p (r )ϕq (z ) ∗ ∗ r + ∗2 − T r + φ p ∗2 r∗ dr∗ dz∗
r ∂r ∂ r∗ ∂z r∗ dr∗ dr∗ dz
z∗ =0 r∗ =0
z∗Z=L∗ r∗ =1 rZ∗ =1 z∗ =L∗
∗ ∂T∗
  
∗ ∗ ∂ T ∗ dφ p ∗ ∗ ∗ dϕq
= ϕq (z ) r φ p ∗ − T dz + φ p (r ) ϕq ∗ − T dr∗
∂r dr∗ ∗
r =0 ∂ z dz∗
∗
z =0
z∗ =0 r∗ =0
(4.3.18)
4.3 Conducción de calor transitoria en una barra cilíndrica 85

Sustituyendo las ecuaciones diferenciales y condiciones de frontera tanto para T ∗ como para φ p y
φq se obtiene lo siguiente
z∗Z=L∗ rZ∗ =1 z∗Z=L∗
∂T∗

+ λ p2 + βq2 T ∗ φ p ϕq r∗ dr∗ dz∗ = Biφ p (1) ϕq (z∗ )dz∗


∂t ∗
z∗ =0 r∗ =0 z∗ =0
r∗ =1
Z
+ BiT1∗ [ϕq (L∗ ) + ϕq (0)] φ p (r∗ )dr∗ (4.3.19)
r∗ =0

Sustituyendo la ecuación (4.3.17) en el resultado anterior da lugar a

∞ ∞   z∗Z=L∗ rZ∗ =1 z∗Z=L∗

d fnm
+ λ p2 + βq2 fnm ϕm ϕq dz∗ φn φ p r∗ dr∗ = Biφ p (1) ϕq (z∗ )dz∗


∑∑ dt ∗
n=1 m=1
z∗ =0 r∗ =0 z∗ =0
r∗ =1
Z
+ BiT1∗ [ϕq (L∗ ) + ϕq (0)] φ p (r∗ )dr∗ (4.3.20)
r∗ =0

O bien, tomando en cuenta las fórmulas de ortogonalidad de las funciones propias:

rZ∗ =1
φn φ p r∗ dr∗ = 0, ∀n 6= p (4.3.21a)
r∗ =0
z∗Z=L∗
ϕm ϕq dz∗ = 0, ∀m 6= q (4.3.21b)
z∗ =0

se obtiene la siguiente ecuación diferencial ordinaria para fnm

z∗R
=L∗ r∗R=1
Biφn (1) ϕm (z∗ )dz∗ + BiT1∗ [ϕm (L∗ ) + ϕm (0)] φn (r∗ )dr∗
d fnm z∗ =0 r∗ =0
+ λn2 + βm2 fnm =


dt ∗ z∗R
=L∗ r∗R=1
ϕm2 dz∗ φn2 r∗ dr∗
z∗ =0 r∗ =0
(4.3.22)
Note que el lado derecho de la ecuación anterior está expresado en términos de cantidades conocidas.
Por ello, con el fin de simplificar la notación, se define a la siguiente constante
z∗R
=L∗ r∗R=1
Biφn (1) ϕm (z∗ )dz∗ + BiT1∗ [ϕm (L∗ ) + ϕm (0)] φn (r∗ )dr∗
z∗ =0 r∗ =0
γnm = z∗ R
=L∗ r∗R=1
(4.3.23)
ϕm2 dz∗ φn2 r∗ dr∗
z∗ =0 r∗ =0

De esta forma, la ecuación (4.3.22) se escribe como

d fnm
+ λn2 + βm2 fnm = γnm


∗
(4.3.24)
dt
Para deducir la condición inicial de esta ecuación diferencial se sustituye la ecuación (4.3.17) en la
ecuación (4.3.3f) para obtener
∞ ∞
t ∗ = 0, ∑∑ fnm (t ∗ )φn (r∗ )ϕm (z∗ ) = 0 (4.3.25)
n=1 m=1
86 Capítulo 4. Transferencia de calor en más de una dirección

De donde, al usar las condiciones de ortogonalidad, se obtiene que

t ∗ = 0, fnm = 0 (4.3.26)

La ecuación (4.3.24) puede reformularse, por medio del método del factor integrante, como sigue

d
exp (λn2 + βm2 )t ∗ fnm = γnm exp (λn2 + βm2 )t ∗
 
 
∗
(4.3.27)
dt
Integrando semidefinidamente la ecuación anterior, tomando en cuenta la condición inicial, resulta
que
γnm
fnm (t ∗ ) = 1 − exp −(λn2 + βm2 )t ∗
 

(4.3.28)
λn + βm2
2

De esta forma, la solución del problema se puede expresar como

∞ ∞
γnm
T ∗ (r∗ , z∗ ,t ∗ ) = φ (r∗ )ϕm (z∗ ) 1 − exp −(λn2 + βm2 )t ∗
 

∑∑ 2 2 n
(4.3.29)
n=1 m=1 λn + βm

Como puede notarse, la solución consiste de una porción que depende del tiempo y otra que no
depende del tiempo (parte estacionaria). Esta última se obtiene al tomar el límite de la ecuación
anterior cuando t ∗ → ∞:
∞ ∞
γnm
T ∗ (r∗ , z∗ ) = ∑∑ 2
φ (r∗ )ϕm (z∗ )
2 n
(4.3.30)
λ
n=1 m=1 n + βm

Es importante notar que el numerador de la constante γnm contiene la influencia de las fuentes del
problema adimensional.

Ejercicio 4.8 — Solución alternativa. Un método más corto al mostrado aquí consiste en
proponer que la solución sea
∞ ∞
T∗ = ∑∑ fnm (t ∗ )Ψnm (r∗ , z∗ ) (4.3.31)
n=1 m=1

Donde Ψnm es una función propia en 2D que satisface el siguiente problema de Sturm-Liouville:

∇2 Ψnm + ωnm2
Ψnm = 0 (4.3.32a)
∂ Ψnm
r∗ = 0, =0 (4.3.32b)
∂ r∗
∂ Ψnm
r∗ = 1, − = BiΨnm (4.3.32c)
∂ r∗
∂ Ψnm
z∗ = 0, = BiΨnm (4.3.32d)
∂ z∗
∂ Ψnm
z∗ = L∗ , − = BiΨnm (4.3.32e)
∂ z∗
Resuelva este último problema con el método de separación de variables y demuestre que
Ψnm (r∗ , z∗ ) = φn (r∗ )ϕm (z∗ ) y ωnm
2 = λ 2 + β 2 . Posteriormente, utilice la solución dada en la
n m
ecuación (4.3.31) junto con la fórmula de Green en 2D para recuperar la solución dada en la
ecuación (4.3.29). 
 4.4 Transporte en medios semi-infinitos 87

4.4 Transporte en medios semi-infinitos

Considere la dinámica de la transferencia de calor en una placa que inicialmente se encontraba a
una temperatura uniforme T0 y que repentinamente una de sus caras se calienta de manera uniforme
y adquiere una temperatura T1 > T0 . Si las otras caras están aisladas, determine los perfiles de
temperatura a tiempos cortos, esto es, a escalas de tiempo cercanas a la condición inicial.
De acuerdo al planteamiento del problema, la ecuación de energía térmica en la placa es:

∂T
ρc p = k∇2 T (4.4.1)
∂t
Para reducir la dimensión del problema definimos a la temperatura promedio como:
y=H
Z z=H
1
Z
yz
hT i = T dydz (4.4.2)
HW
y=0 z=0

Dado que todas las caras están aisladas excepto una, tras aplicar el operador de promediado, la
ecuación de energía térmica se reduce a:

∂ hT iyz ∂ 2 hT iyz
=α (4.4.3a)
∂t ∂ x2
La cual está sujeta a las siguientes condiciones de frontera e inicial:

En x = 0, hT iyz = T1 (4.4.3b)
yz
En x → ∞, hT i = T0 (4.4.3c)
yz
Cuando t = 0, hT i = T0 (4.4.3d)

Dado que el análisis está restringido a tiempos cortos, la segunda condición de frontera se utiliza
para representar que el resto de la placa aún no ha sido afectada por el cambio repentino en la
temperatura que ocurrió en la superficie x = 0. Para resolver este problema es conveniente definir a
la temperatura adimensional como:

hT iyz − T0
T∗ = (4.4.4)
T1 − T0

De esta forma, la formulación del problema en términos de la temperatura adimensional es:

∂T∗ ∂ 2T ∗
=α (4.4.5a)
∂t ∂ x2
En x = 0, T ∗ = 1 (4.4.5b)
∗
En x → ∞, T = 0 (4.4.5c)
∗
Cuando t = 0, T = 0 (4.4.5d)

Aunque este problema puede, en principio, resolverse de manera numérica o analítica con el método
de separación de variables; en los siguientes párrafos se examina un tipo de solución que permitirá
transformar el problema de valor inicial y a la frontera a un problema de valor a la frontera con una
ecuación diferencial ordinaria. Se obtendrá la solución exacta del problema y posteriormente se
explorará la solución aproximada con el método de von Kármán-Pohlhausen de manera aproximada.
Con estas herramientas se contará con las herramientas necesarias para estudiar el problema de
capa límite entre un fluido y un obstáculo.
 88 Capítulo 4. Transferencia de calor en más de una dirección

4.4.1 Solución exacta

La transformada de Laplace es un operador integral que se aplica por lo regular a problemas
transitorios y se define como
Z∞
L { f } ≡ f (s) = e−st f (t) dt (4.4.6)
t=0

donde s es un número complejo que es constante en la posición y sus unidades son 1/s. Algunas
fórmulas útiles de transformadas de Laplace se presentan a continuación:

f (t) f (s)
∂f
∂t s f − f (0)
u=t
R
f (u)g(t − u) du f (s)g(s)
u=0
a
a (cte) s √
exp(−a2 /(4t))
√ exp(−a
√ s)
πt  s√
a exp(−a s)
erfc 2√ s
√ t  √
a exp(−a
√ √ s)
exp[b(bt + a)]erfc b t + 2√ t s( s+b)

En la tabla anterior la segunda fila se conoce como teorema de la convolución; mientras que en las
últimas filas, la función error complementaria se define como
ξZ=η
2
erfc(η) = 1 − erf(η) = 1 − √ exp(−ξ 2 )dξ (4.4.7)
π
ξ =0

Dado que
ξZ=∞
√
exp(−ξ 2 )dξ = π (4.4.8)
ξ =−∞

se tiene que
ξZ=∞
2
erfc(η) = √ exp(−ξ 2 )dξ (4.4.9)
π
ξ =η

Con este resultado es fácil deducir que erfc(−∞)=2; erfc(0)=1 y erfc(∞)=0; es decir, es una función
decreciente del argumento η.
Note que, al usar la regla de Leibniz,
x=b(x) x=b(x)
d df db da
Z Z
f (x) dx = dx + f (b) − f (a) (4.4.10)
dx dx dx dx
x=a(x) x=a(x)

se obtiene que

derfc(η) 2
= − √ exp(−η 2 ) (4.4.11)
dη π
4.4 Transporte en medios semi-infinitos 89

Aplicando el operador de transformada de Laplace a la ecuación (4.4.5a) resulta que

d2T ∗ s ∗
− T =0 (4.4.12a)
dx2 α
La cual está sujeta a las siguientes condiciones de frontera
En x = 0, T ∗ = 1/s (4.4.12b)
En x → ∞, T∗ =0 (4.4.12c)
La solución general de la ecuación diferencial es
r   r 
∗
s s
T = A exp x + B exp − x (4.4.13)
α α
Para satisfacer la condición de frontera en x → ∞, debe cumplirse que A = 0. De esta forma,
considerando la condición de frontera en x = 0 resulta que B = 1/s y la solución es
exp − αs x
p

∗
T = (4.4.14)
s
Cuya transformada inversa es
 
∗ x
T = erfc √ (4.4.15)
2 αt
En este momento es conveniente definir al espesor de la capa límite, δT ,como el lugar donde ocurren
el 99 % de los cambios. En este caso, dicho espesor corresponde al valor de x para el cual T ∗ = 0.01.
Usando tablas de la función error1 puede determinarse que el valor de 2√xαt que corresponde a
T ∗ = 0.01 es:
x δ
√ = √ = 1.83 (4.4.16)
2 αt T ∗ =0.01 2 αt
Es decir,
√
δT = 3.66 αt (4.4.17)
Con la solución dada en la ecuación (4.4.15) se puede obtener el flux de calor en la superficie
x = 0 como
∂ hT iyz ∂T∗ k(T1 − T0 )
hqx iyz |x=0 = −k = −k(T1 − T0 ) = √ (4.4.18)
∂ x x=0 ∂ x x=0 παt
Al comparar este resultado con la ley de enfriamiento de Newton:
qx |x=0 = h(T1 − T0 ) (4.4.19)
permite concluir que:
k
h= √ (4.4.20)
παt
Este resultado es importante pues permite calcular los valores del coeficiente interfacial de transfe-
rencia de calor en una interfase como función del tiempo. Note que si este resultado se escribe en
forma adimensional, lo que se obtiene es la definición del número de Nusselt:
h`re f `re f
Nu = =√ (4.4.21)
k παt
donde `re f es una longitud de referencia. Antes de examinar un método aproximado de obtener la
solución del problema, se examina otro caso de estudio.
1 Ver por ejemplo: https://www.geophysik.uni-muenchen.de/∼malservisi/GlobaleGeophysik2/erf_tables.pdf
 90 Capítulo 4. Transferencia de calor en más de una dirección

4.4.2 Flux de calor en la frontera

Reconsidere la situación física descrita arriba para el caso en el que se especifique el flux de
calor en x = 0; esto es,
∂ hT iyz
En x = 0, k = −q0 (4.4.22)
∂x
Como ahora no existe la temperatura T1 , la definición de T ∗ es
hT iyz − T0
T∗ = (4.4.23)
T0
Bajo estas condiciones, el problema es ahora
∂T∗ ∂ 2T ∗
= α (4.4.24a)
∂t ∂ x2
∂T∗ q0
En x = 0, k =− (4.4.24b)
∂x T0
∗
En x → ∞, T = 0 (4.4.24c)
∗
Cuando t = 0, T =0 (4.4.24d)

Aplicando la transformada de Laplace a la ecuación diferencial se recupera la ecuación (4.4.12a),

cuya solución general está dada en la ecuación (4.4.13). La cual se reduce a la siguiente expresión
al tomar en cuenta que la temperatura es finita cuando x → ∞:
x √
 
∗
T = B exp − √ s (4.4.25)
α
La condición de frontera dada en la ecuación (4.4.24b) se escribe en el dominio de Laplace como
sigue
∂T∗ q0
En x = 0, =− (4.4.26)
∂x kT0 s
Aplicando esta condición de frontera resulta que
√ exp − √x √s
 
q0 α α
T∗ = √ (4.4.27)
kT0 s s
Para regresar al dominio del tiempo el resultado anterior se utiliza el teorema de la convolución
√ Zu=t
  √ 
  exp − √x s 
q0 α 1 
α
T∗ = L −1 L −1 √ du (4.4.28)
kT0 s t−u  s 
u=0 u

O bien, tomando en cuenta los resultados de la tabla de transformada de Laplace proporcionada

arriba,
 
p Zu=t exp − x2
q0 α/π 4αu
T∗ = √ du (4.4.29)
kT0 u
u=0

Para resolver esta integral se recurre al método de integración por partes y el resultado es
 
 Zu=t
q0 α/π  √
p
2 2 2
  
x x x
T∗ = 2 t exp − − 3/2
exp − du (4.4.30)
kT0 4αt 2αu 4αu
u=0
 4.4 Transporte en medios semi-infinitos 91

En este punto vale la pena hacer el cambio de variable

x x
ξ= √ ; dξ = − √ 3/2 du (4.4.31)
2 αu 4 αu
De manera que
ξ = 2√xαt
 
Zu=t √ ξZ=∞
x2 x2
 
2x x π 2
Z
exp − du = − √ exp(−ξ 2 )dξ = √  √ exp(−ξ 2 )dξ 

2αu3/2 4αu α α π
u=0 ξ =∞ ξ = 2√xαt
r  
π x
=x erfc √ (4.4.32)
α 2 αt
En la última igualdad se utilizó la definición de la función error complementaria dada en la ecuación
(4.4.9). Sustituyendo este resultado en la ecuación (4.4.30) se obtiene que
p
x2
   
∗ 2q0 αt/π q0 x x
T = exp − − erfc √ (4.4.33)
kT0 4αt kT0 2 αt

4.4.3 Ley de enfriamiento en la frontera

Por último, considere el caso en el que la superficie en x = 0 es aplicable una ley de enfriamiento
de Newton, esto es,
∂ hT iyz
En x = 0, k = h(T − T∞ ) (4.4.34)
∂x
Si se considera que T∞ > T0 , entonces la temperatura adimensional está dada por
hT iyz − T0
T∗ = (4.4.35)
T∞ − T0
Por lo que el problema en términos de la temperatura adimensional es

∂T∗ ∂ 2T ∗
= α (4.4.36a)
∂t ∂ x2
∂T∗
En x = 0, k = h(T ∗ − 1) (4.4.36b)
∂x
En x → ∞, T ∗ = 0 (4.4.36c)
∗
Cuando t = 0, T =0 (4.4.36d)

La solución escrita en la ecuación (4.4.25) sigue siendo válida en este caso y está sujeta a la
siguiente condición de frontera

∂T∗
 
∗
1
En x = 0, k =h T − (4.4.37)
∂x s
De esta forma, se deduce que la solución es
 √
exp − √xα s
T∗ = b √ (4.4.38)
s( s + b)
En donde, para simplificar la notación se definió
√
h α
b= (4.4.39)
k
 92 Capítulo 4. Transferencia de calor en más de una dirección

Usando el método de fracciones parciales (ver ejercicio 4.9) puede demostrarse que
b 1 1
√ = −√ √ (4.4.40)
s( s + b) s s( s + b)
Por lo tanto, la ecuación (4.4.38) se puede expresar como
 √  √
exp − √xα s exp − √xα s
T∗ = − √ √ (4.4.41)
s s( s + b)
Usando las fórmulas dadas en la tabla de transformada de Laplace se deduce que el resultado es
   2   √ 
∗ x h αt hx h αt x
T = erfc √ − exp + erfc + √ (4.4.42)
2 αt k2 k k 2 αt
Con esta solución se puede calcular el flux en la superficie x = 0 definido como

∂T∗
hqx iyz |x=0 = −k = h (1 − T ∗ |x=0 ) (4.4.43)
∂x x=0

Ya que erfc(0) = 0, el resultado es

 2   √ 
yz h αt h αt
hqx i |x=0 = h exp erfc (4.4.44)
k2 k

√
Ejercicio 4.9 — Método de fracciones parciales. Adopte el cambio de variable s = w;
de tal forma que

b b A B C
√ = = + + (4.4.45)
s( s + b) w2 (w + b) w w2 w + b

Demuestre que A = −1/b; B = 1 y C = 1/b. De esta forma, la ecuación anterior se escribe

como
 
b 1 1 1 1
√ = − √ −√ (4.4.46)
s( s + b) s b s s+b

Lleve a cabo las operaciones correspondientes para expresar el resultado anterior en la forma
dada en la ecuación (4.4.40).


4.4.4 Solución aproximada par el método de von Kármán-Pohlhausen

Una forma alternativa y aproximada de resolver el problema consiste en observar que la solución
del problema es parabólica, por lo que se propone una solución de tipo polinómica:

x x2 x3
T ∗ = a0 + a1 + a2 2 + a3 3 (4.4.47)
δT δT δT
Este problema se resolverá únicamente en el dominio definido por la capa límite para el primer
caso discutido en los párrafos anteriores. Por tanto, la solución debe satisfacer:

En x = 0, T∗ = 1 (4.4.48a)
∗
En x = δ , T ≈0 (4.4.48b)
4.4 Transporte en medios semi-infinitos 93

Con sólo dos condiciones de frontera no es posible calcular todas las constantes del polinomio.
Para atender este problema, se impone que, después de la capa límite, no hay transferencia de calor,
por tanto:
∂T∗
En x = δ , =0 (4.4.48c)
∂x
Por último, dado que en x = 0 la temperatura es constante, la ecuación diferencial permite deducir
que

∂ 2T ∗
En x = 0, =0 (4.4.48d)
∂ x2
Considerando estas condiciones de frontera resulta que a0 = 1, a1 = −3/2, a2 = 0, a3 = 1/2 y
entonces:
3 x 1 x3
T∗ = 1− + (4.4.49)
2 δT 2 δT3
Para que esta solución esté completa es necesario conocer el espesor de la capa límite. Para
determinar esta incógnita se sustituirá la solución dada en la ecuación anterior, en la versión
integrada en la posición de la ecuación diferencial:
x=δ
Z T
∂T∗ ∂T∗ 3α
dx = −α = (4.4.50)
∂t ∂x x=0 2δT
x=0

La derivada temporal puede expresarse, usando la regla de la cadena, como sigue

∂T∗ 3 x x3 dδ
 
= − (4.4.51)
∂t 2 δT2 δT4 dt
Integrando en el espacio el resultado anterior se obtiene que
x=δ
Z T
∂Θ 3 dδ
dx = (4.4.52)
∂t 8 dt
x=0

Entonces, la ecuación diferencial para el espesor de la capa límite es:

dδT
δT = 4α (4.4.53)
dt
La cual está sujeta a la condición inicial:

Cuando t = 0, δT = 0 (4.4.54)

ya que en un inicio no había transferencia de calor.

La solución de este problema es:
√ √
δT = 2 2αt = 2.8 αt (4.4.55)

Aunque este resultado es claramente distinto del resultado de la solución exacta, vale la pena
continuar y calcular el flux de calor en x = 0:
∂Θ 3k
qx |x=0 = −k(T1 − T0 ) = √ (T1 − T0 ) (4.4.56)
∂x x=0 5.6 αt
94 Capítulo 4. Transferencia de calor en más de una dirección

lo cual lleva a concluir que:

0.536k
h= √ , modelo aproximado (4.4.57)
αt
Este resultado es de hecho muy similar al obtenido anteriormente:

0.564k
h= √ , modelo exacto (4.4.58)
αt
De lo anterior se deduce que este método es especialmente práctico para predecir de manera
aproximada el coeficiente interfacial de transferencia de calor.

Ejercicio 4.10 — Método aproximado con flux de calor en la frontera. Utilice el método
aproximado explicado aquí para predecir la dinámica del espesor de la capa límite para el caso
en el que se especifica el flux de calor en la superficie x = 0. 

Ejercicio 4.11 — Método aproximado con ley de enfriamiento. Aplique el método apro-
ximado explicado aquí para el caso en el que se puede aplicar una ley de enfriamiento de Newton
en la superficie x = 0. 
 5. Soluciones numéricas

La evaluación de la solución numérica del problema se muestra en la Figura 4.2. Como puede
notarse el perfil de temperatura es simétrico alrededor del eje x. Además en la Figura ?? se muestra
la evaluación de la temperatura promedio hΘi y del promedio de la componente normal a las
superficies de entrada y de salida. Como es de esperarse, a medida que aumenta la superficie de la
viga, aumentan tanto la temperatura promedio como el flux de calor.
Suponga ahora que en la superficie de la derecha de la viga se coloca una fuente oscilante de
calor de manera que la temperatura en esta superficie obedece la siguiente expresión:

En X = X0 , T (Y ) = 1 + 0.5 sin(8πY ) (5.0.1)

Los resultados de tomar esta condición de frontera se muestran en la Figura 5.1. Como puede
notarse la influencia de las oscilaciones en la frontera llega a afectar los perfiles de temperatura y
del flux de calor cerca de la segunda mitad de la viga. Después de ahí, los perfiles son similares a
los observados en el caso anterior.

5.1 Coordenadas cilíndricas

5.1.1 Transporte axisimétrico
Considere ahora un cilindro de radio R en donde se lleva a cabo la transferencia de calor por
conducción de manera axisimétrica. En la base del tubo (z = 0) la temperatura es constante y
conocida (T0 ) y en el extremo opuesto (z = L) la temperatura es constante pero no es conocida.
Además, la superficie externa está expuesta al ambiente y ofrece resistencias a la transferencia
de calor. Se desea evaluar la dependencia de la temperatura y el flux de calor con el coeficiente
interfacial de transferencia de calor h.
En este caso, la ecuación de energía térmica es:

∂ 2T
 
1 ∂ ∂T
0= r + 2 (5.1.1a)
r ∂r ∂r ∂z
La cual está sujeta a las siguientes condiciones de frontera:
96 Capítulo 5. Soluciones numéricas

Figura 5.1: Perfiles adimensionales de temperatura y de la componente y del flux de calor a

diferentes posiciones. En todos los casos X0 =0.5
 5.1 Coordenadas cilíndricas 97

∂T
En r = 0, =0 (5.1.1b)
∂r
∂T
En r = R, − k = h(T − Ta ) (5.1.1c)
∂r
En z = 0, T = T0 (5.1.1d)
∂T
En z = L, =0 (5.1.1e)
∂z
Para resolver este problema es conveniente expresarlo en términos de las siguientes variables y
parámetros adimensionales:

T − Ta r z L hR
Θ= ; ξ= ; ζ= ; λ= ; Bi = (5.1.2)
T0 − Ta R R R k
De esta forma, el problema de valor a la frontera es ahora:

∂ 2Θ
 
1 ∂ ∂Θ
0= ξ + 2 (5.1.3a)
ξ ∂ξ ∂ξ ∂ζ
∂Θ
En ξ = 0, =0 (5.1.3b)
∂ξ
∂Θ
En ξ = 1, − = BiΘ (5.1.3c)
∂ξ
En ζ = 0, Θ = 1 (5.1.3d)
∂Θ
En ζ = λ , =0 (5.1.3e)
∂ζ

En la Figura 5.2 se muestra la dependencia de los perfiles de temperatura con el número de Biot,
el cual relaciona el transporte por convección en la frontera respecto al transporte por conducción
dentro del cilindro. En todos los casos se tomó λ = 2. Como puede notarse, cuando este parámetro
es mucho menor que la unidad el cilindro se encuentra prácticamente a la misma temperatura que
en la base pues el calor no se disipa en la frontera. De esta forma, a medida que se incrementan los
valores del número de Biot se establecen más fácilmente los perfiles de temperatura a lo largo del
sistema.

5.1.2 Transporte no axisimétrico

Reconsidere el problema visto anteriormente con la diferencia que ahora sólo un hemisferio del
tubo está expuesto al ambiente y el otro hemisferio está sometido a una fuente variable de calor. De

Figura 5.2: Perfiles de temperatura para Bi = 1 × 10−2 , 0.1, 1, 100.

98 Capítulo 5. Soluciones numéricas

Figura 5.3: Dependencia de la temperatura promedio tomada en las superficies lateral y superior y
del flux de calor en la dirección z con el número de Biot.

esta forma, tenemos que:

∂T
En r = R, 0 ≤ θ ≤ π, −k = h(T − Ta ) (5.1.4a)
∂r
∂T Q̇
En r = R, π ≤ θ ≤ 2π, −k = (5.1.4b)
∂r A

Y el resto del problema permanece inalterado. En forma adimensional la última condición de

frontera puede expresarse como:

∂Θ Q̇R 1 Ω
En ξ = 1, π ≤ θ ≤ 2π, − = = (5.1.5)
∂ξ k(T0 − Ta ) A A

La dependencia de los perfiles de temperatura con la fuente de calor se muestra en la Figura 5.4
para un valor de Bi = 10. Como es de esperarse, en este caso la temperatura aumenta de manera
lineal con la fuente de calor, lo cual puede verificarse al tomar el promedio de la temperatura en
cualquier superficie del cuerpo y evaluarlo para diferentes valores de Ω.
 5.2 Problema de aplicación 99

5.2 Problema de aplicación

Los tubos de cobre que llevan agua a 180o F están embebidos en el espacio que se encuentra
entre el piso y el techo de varios edificios. El diámetro externo de los tubos es de 1 in y se supone
despreciable la resistencia externa a la transferencia de calor. Si la superficie superior se encuentra
a 60o F y la superficie inferior está a 83o F. La conductividad térmica del material puede tomarse
como 0.42 [Btu/(h ft o F)]. Determine los perfiles de temperatura y los valores del flux de calor en
las superficies superior e inferior. A continuación se muestra un esquema del sistema y también
los perfiles de temperatura. Los valores de la integral del flux de calor en las superficies superior e
inferior son 97.6 y 67.85 W/m, respectivamente.

5.3 Aleta de enfriamiento revisada

En nuestro estudio previo de la aleta de enfriamiento hicimos varias suposiciones para obtener
soluciones aproximadas. En este problema se busca tener una visión más realista de las aletas de
enfriamiento. Para ello consideremos una placa vertical de largo L, altura H y ancho W , en donde
 100 Capítulo 5. Soluciones numéricas

H = 5L y W = 2L. En este problema se desea comparar los perfiles de temperatura con y sin aletas
para evaluar el beneficio que significa usar aletas de enfriamiento. Para ello considere 5 aletas de
enfriamiento, cada una de altura 0.5L (o bien menor), ancho 2L y largo 2L. Considerando que las
superficies inferior y lateral izquierda de la placa se encuentran a temperatura constante T0 y que
el resto de las superficies están expuestas al ambiente. Evalúe el desempeño con y sin aletas de
enfriamiento para diversos valores del número de Biot. Suponga que las aletas de enfriamiento
tienen una conductividad térmica que es el doble de la conductividad térmica de la placa y que el
coeficiente interfacial de transferencia de calor es el mismo en la placa y en las aletas.
Las ecuaciones diferenciales que rigen la transferencia de calor en la placa y en las aletas son:

0 = ki ∇2 Ti , i = p, a (5.3.1a)

En las superficies inferior y lateral izquierda (Ω1 ), la condición de frontera es:

En Ω1 , Tp = T0 (5.3.1b)

En el resto de las superficies (Ω2 ) aplica la ley de enfriamiento:

En Ω2 , −n · ki ∇Ti = h(Ti − Ta ), i = a, p (5.3.1c)

Para llevar a cabo el análisis es conveniente definir las siguientes variables y parámetros
adimensionales:
r Ti − Ta ka hL
r∗ = , Θi = , κ p = , Bi = (5.3.2)
L T0 − Ta kp kp

De esta forma, el problema adimensional es:

0 = κi ∇∗2 Θi , i = p, a (5.3.3a)
En Ω∗1 , Θp = 1 (5.3.3b)
En Ω∗2 , −n · κi ∇∗ Θi = BiΘi , i = a, p (5.3.3c)

Note que κ p = 1. En la Figura 5.5 se muestran ejemplos de los perfiles de temperatura y de la

temperatura promedio en la superficie de la placa donde están las aletas para diferentes valores del
número de Biot considerando dos alturas de las aletas de enfriamiento. Como es de esperarse, a
mayor superficie de contacto, mejor disipación de calor.
 Ejemplo 5.1 Reconsidere el caso en el que el alto de la aleta sea 0.1L pero ahora en lugar de 5
aletas, utilice 10 aletas. Contraste sus resultados con los obtenidos en este ejercicio. 

 Ejemplo 5.2 Estudie el mismo problema con aletas de enfriamiento en coordenadas cilíndricas. 

5.4 Conducción de calor transitoria

Considere la transferencia de calor en una placa que inicialmente se encuentraba a una tem-
peratura uniforme T0 ; posteriormente la superficie en x = 0 se expone al ambiente presentando
resistencias interfaciales a la transferencia de calor, mientras que el resto de las superficies se
encuentran aisladas. La placa tiene un largo L, una altura H y un ancho W . Siendo H = 5L y
W = 8L.
La ecuación diferencial que rige este proceso es:

∂T
ρCp = k∇2 T (5.4.1a)
∂t
5.4 Conducción de calor transitoria 101

Figura 5.4: Dependencia de los perfiles de temperatura con la fuente de calentamiento Ω = 25, 50,
75 y 100 m2 . En todos los casos se tomó Bi = 10.

a) b)

c) d)

Figura 5.5: a) Ejemplos de los campos adimensionales de temperatura con y sin aleta para Bi = 1
b) Dependencia de la temperatura promedio en la superficie de la placa calculada con y sin aleta
de enfriamiento para una altura de aleta de 0.5L como función del número de Biot. En c) y d) se
presentan los resultados correspondientes a una altera de alto 0.1L.
 102 Capítulo 5. Soluciones numéricas

La cual está sujeta a las siguientes condiciones de frontera:

∂T
En x = 0, k = h(T − Ta ) (5.4.1b)
∂x
El resto de las superficies están aisaldas:

En ∂ Ω, n · ∇T = 0 (5.4.1c)

Por último, la condición inicial es:

Cuando t = 0, T = T0 (5.4.1d)

Para escribir este problema en forma adimensional se proponen las siguientes definiciones:
tk r T − Ta hL
t∗ = ; r∗ = ; Θ= ; Bi = (5.4.2)
ρCp L2 L T0 − Ta k
De esta forma, el problem adimensional es:
∂Θ
= ∇∗2 Θ (5.4.3a)
∂t ∗
∂Θ
En x∗ = 0, = BiΘ (5.4.3b)
∂ x∗
En ∂ Ω∗ , n · ∇∗ Θ = 0 (5.4.3c)
∗
Cuando t = 0, Θ = 1 (5.4.3d)

Usando el proceso de promediado ya visto, y suponiendo que la temperatura en la superficies es

igual a la temperatura promedio, el modelo aproximado es:
dhΘi
+ BihΘi = 0 (5.4.4a)
dt ∗
Cuando t ∗ = 0, hΘi = 1 (5.4.4b)

Cuya solución es:

∗
hΘi = e−Bit (5.4.5)

En la Figura 5.6 comparamos esta solución aproximada con los resultados de la solución numérica
del problema tridimensional en estado transitorio para distintos números de Biot. Como puede
notarse esta aproximación sólo es razonable cuando este parámetro es bajo.

5.5 Conducción y convección en estado transitorio

Considere una esfera (fase-σ ) de radio R que se encuentra en contacto con un fluido (fase-
β ) en un tubo de radio R2 = 5R y altura H = 15R. Inicialmente, la esfera se encontraba a una
temperatura T 0, la cual es el mayor que la temperatura del fluido T1 . En un tiempo > 0, el fluido
cambia su temperatura de entrada a la media aritmética entre T0 y T1 . Determine el tiempo al cual
la temperatura promedio del sólido adquiere la nueva temperatura del fluido. Considere que el
Reynolds vale 1000 y que el número de Péclet puede valer 1, 10, 100, 1000 y 10000. Puede suponer
que el flujo es Newtoniano, incompresible y en estado estacionario y que la difusividad térmica del
sólido es el doble de la del fluido. Además, puede suponer que las paredes del tubo están aisladas,
que a la salida del tubo la temperatura del fluido es constante y que la esfera no presenta resistencias
a la transferencia de calor con el fluido.
5.5 Conducción y convección en estado transitorio 103

1
Bi = 0.1
0.8 10 1

0.6
hΘi
0.4

0.2 102

10−4 10−3 10−2 10−1 100 101

t∗

Figura 5.6: Comparación de la dinámica de la temperatura promedio de la placa obtenida de manera

aproximada (líneas punteadas) y exacta (líneas continuas).

El problema está descrito por las ecuaciones gobernantes de la transferencia de cantidad de

movimiento y calor. Para resolver el problema de mecánica de fluidos es necesario considerar las
ecuaciones de continuidad y Navier-Stokes:

∇ · vβ = 0, en la fase − β (5.5.1a)
2
ρβ vβ · ∇vβ = −∇pβ + µβ ∇ vβ , en la fase − β (5.5.1b)

Las condiciones de frontera son:

En z = 0, pβ = pin (5.5.1c)
En z = H, pβ = pout (5.5.1d)
En r = R, R2 , vz = 0 (5.5.1e)

Mientras que, para la transferencia de calor, la ecuación de energía térmica en cada fase es:

∂ Tσ
ρσ c pσ = kσ ∇2 Tσ , en la fase-σ (5.5.2a)
 ∂t 
∂ Tβ
ρβ c pβ + vβ · ∇Tβ = kβ ∇2 Tβ , en la fase-β (5.5.2b)
∂t

Las condiciones de frontera son:

1
En z = 0, Tβ = (T0 + T1 ) (5.5.2c)
2
∂ Tβ
En z = H, =0 (5.5.2d)
∂z
En r = R, Tβ = Tσ (5.5.2e)
∂ Tβ ∂ Tσ
En r = R, −kβ = −kσ (5.5.2f)
∂r ∂r
∂ Tβ
En r = R2 , −kβ =0 (5.5.2g)
∂r
104 Capítulo 5. Soluciones numéricas

Por último, las condiciones iniciales son:

Cuando t = 0, Tσ = T0 (5.5.2h)
Cuando t = 0, Tβ = T1 (5.5.2i)
Para escribir las ecuaciones en forma adimensional, definimos las siguientes variables y parámetros:
vβ ∗ pβ − pout Ti − T1 ∗ r ∗ tαβ ασ kσ
v∗β = , pβ = , Θi = , r = ,t = 2 , ασ β = , kσ β =
vre f pin − pout T0 − T1 R R αβ kβ
ρβ vre f R vre f R
Re = , Pe = (5.5.3)
µβ αβ
De esta forma, el problema en forma adimensional es:
∇∗ · v∗β = 0, en la fase − β (5.5.4a)
Rev∗β · ∇∗ v∗β = −∇∗ p∗β + ∇∗2 v∗β , en la fase − β (5.5.4b)
Las condiciones de frontera son:
En z∗ = 0, p∗β = 1 (5.5.4c)
∗
En z = 15, p∗β = 0 (5.5.4d)
En r∗ = 1, 5, v∗z = 0 (5.5.4e)
Note que para escribir estas ecuaciones se tomó:
(pin − pout ) R
vre f = (5.5.5)
µβ
La forma adimensional del problema de transferencia de calor es:
∂ Θσ
= ασ β ∇2 Θσ , en la fase-σ (5.5.6a)
∂t ∗
∂ Θβ
+ Pev∗β · ∇+ Θβ = ∇∗2 Θβ , en la fase-β (5.5.6b)
∂t ∗
Las condiciones de frontera son:
1
En z∗ = 0, Θβ = (5.5.6c)
2
∂ Θ
En z∗ = H,
β
=0 (5.5.6d)
∂ z∗
En r∗ = 1, Θβ = Θσ (5.5.6e)
∂ Θβ ∂ Θσ
En r∗ = 1, − ∗ = −kσ β + (5.5.6f)
∂r ∂r
∂ Θ
En r∗ = 5,
β
=0 (5.5.6g)
∂ r∗
Por último, las condiciones iniciales son:
Cuando t ∗ = 0, Θσ = 1 (5.5.6h)
∗
Cuando t = 0, Θβ = 0 (5.5.6i)
El resultado de la dinámica de la temperatura promedio de la esfera se muestra en la Figura 5.7.
 Ejemplo 5.3 Resuelva de nuevo el problema considerando resistencias a la transferencia de calor

y considere sólo Pe = 1000 y lleva a cabo un barrido paramétrico del número de Biot. 
5.5 Conducción y convección en estado transitorio 105

Figura 5.7: Dinámica de la temperatura promedio en la esfera para diferentes números de Péclet
tomando Re=1000.
 6. Capa límite

6.1 Ecuaciones de capa límite de energía y cantidad de movimiento para

flujo laminar
Este capítulo está dedicado al análisis de los perfiles de temperatura involucrando convección
y conducción en la capa límite. Para el análisis se utiliza el método de von-Kármán-Pohlhausen.
Primeramente se considera un problema en el que se especifica la temperatura en la superficie y
posteriormente se aborda el caso en el que se especifica el flux de calor en la superficie.

6.1.1 Planteamiento y simplificación

Considere el caso mostrado en la Figura 6.1 en donde un flujo Newtoniano, laminar, incom-
presible y estacionario que tiene una presión p∞ , velocidad v = u∞ ex y temperatura T∞ se hace
pasar alrededor de una placa estacionaria que se encuentra a una temperatura T0 . Suponiendo que la
transferencia de calor se da en estado estacionario, las ecuaciones de conservación de masa total,
cantidad de movimiento y energía térmica son:

∂ vx ∂ vy ∂ vz
+ + =0 (6.1.1a)
∂x ∂y ∂z
 2
∂ vx ∂ 2 vx ∂ 2 vx
  
∂ vx ∂ vx ∂ vx ∂p
ρ vx + vy + vz =− +µ + 2 + 2 (6.1.1b)
∂x ∂y ∂z ∂x ∂ x2 ∂y ∂z
 2
∂ vy ∂ vy ∂ 2 vy
2
  
∂ vy ∂ vy ∂ vy ∂p
ρ vx + vy + vz =− − ρg+ µ + 2 + 2 (6.1.1c)
∂x ∂y ∂z ∂y ∂ x2 ∂y ∂z
 2
∂ vz ∂ vz ∂ 2 vz
2
  
∂ vz ∂ vz ∂ vz ∂p
ρ vx + vy + vz =− +µ + 2 + 2 (6.1.1d)
∂x ∂y ∂z ∂z ∂ x2 ∂y ∂z
 2 2 2

∂T ∂T ∂T ∂ T ∂ T ∂ T
vx + vy + vz =α + 2 + 2 (6.1.1e)
∂x ∂y ∂z ∂ x2 ∂y ∂z

Antes de continuar con las condiciones de frontera, es conveniente explorar la posibilidad de reducir
el problema. Para ello, suponga que el largo, L aunque es mayor que el ancho, h de la placa, es
108 Capítulo 6. Capa límite

Figura 6.1: Esquema de los perfiles de velocidad y de temperatura en sus respectivas capas límites
para el caso del flujo alrededor de una placa.

mucho menor que la profundidad de la placa, W . Esto es:

h<L W (6.1.2)

Más aún, suponga por el momento que todas las componentes de la velocidad son del mismo orden
de magnitud e iguales al valor máximo de la velocidad u∞ . Bajo estas condiciones, se pueden
proponer plantear los siguientes estimados de orden de magnitud:

∂v u 
∞
=O (6.1.3a)
∂x L
∂v u 
∞
=O (6.1.3b)
∂z W
∂ 2v u 
∞
=O 2 (6.1.3c)
∂ x2 L
∂ 2v u 
∞
=O (6.1.3d)
∂ z2 W2
∂ 2T
 
T
=O (6.1.3e)
∂ x2 L2
∂ 2T
 
T
=O (6.1.3f)
∂ z2 W2

Por tanto, dado que L  W , se puede prescindir de la dependencia en z de todas las variables e
incluso de la componente vz del vector de velocidad. Entonces, el sistema de ecuaciones se reduce
6.1 Ecuaciones de capa límite de energía y cantidad de movimiento para flujo
laminar 109
a:
∂ vx ∂ vy
+ =0 (6.1.4a)
∂x ∂y
 2
∂ vx ∂ 2 vx
  
∂ vx ∂ vx ∂p
ρ vx + vy =− +µ + 2 (6.1.4b)
∂x ∂y ∂x ∂ x2 ∂y
 2
∂ vy ∂ 2 vy
  
∂ vy ∂ vy ∂p
ρ vx + vy =− − ρg+ µ + 2 (6.1.4c)
∂x ∂y ∂y ∂ x2 ∂y
 2 2

∂T ∂T ∂ T ∂ T
vx + vy =α + (6.1.4d)
∂x ∂y ∂ x2 ∂ y2

Más aún, tomando en cuenta que los espesores de las capas límite térmica δT e hidrodinámica δH ,
son mucho menores que la longitud característica en la dirección x (es decir L), se tiene que las
ecuaciones anteriores se reducen a:
∂ vx ∂ vy
+ =0 (6.1.5a)
∂x ∂y
∂ 2 vx
 
∂ vx ∂ vx ∂p
ρ vx + vy =− +µ 2 (6.1.5b)
∂x ∂y ∂x ∂y
∂ 2 vy
 
∂ vy ∂ vy ∂p
ρ vx + vy =− − ρg + µ 2 (6.1.5c)
∂x ∂y ∂y ∂y
∂T ∂T 2
∂ T
vx + vy =α 2 (6.1.5d)
∂x ∂y ∂y
En este momento es oportuno corregir el estimado de orden de magnitud respecto al vector de
velocidad. Dado que, antes de chocar con la placa la velocidad era v = u∞ ex , es razonable suponer
que vx = O(e∞ ), dado que x = O(L) y y = O(δH ), a partir de la ecuación de continuidad se tiene
que:
 
δH
vy = O u∞ (6.1.6)
L

Por lo tanto, a partir de las ecuaciones de cantidad de movimiento se deducen los siguientes
estimados de orden de magnitud para las dos componentes del gradiente de presión:
 2 
∂p u∞ u∞
=O ρ +µ 2 (6.1.7a)
∂x L δH
  2   
∂p δH u∞ u∞ δH ∂p
=O ρ +µ 2 = O (6.1.7b)
∂y L L δH L ∂x

Los resultados anteriores pueden usarse en la siguiente expresión:

∂p ∂p
dp = dx + dy (6.1.8)
∂x ∂y

para deducir que, con base a la restricción de escala δH  L, se tiene que:

dp ∂ p
≈ (6.1.9)
dx ∂x
Lo cual indica que la presión sólo puede cambiar en la dirección x,. Sin embargo, dado que en
y = δH , p = p∞ y la presión no cambia en la dirección y, entonces se puede concluir que la presión
 110 Capítulo 6. Capa límite

no cambia de valor en toda la capa límite (esto es, d p/dx = 0, en la capa límite). Lo anterior
hace innecesario resolver la ecuación de cantidad de movimiento para vy . Entonces, el problema
simplificado es:
∂ vx ∂ vy
+ =0 (6.1.10a)
∂x ∂y
∂ vx ∂ vx ∂ 2 vx
vx + vy =ν 2 (6.1.10b)
∂x ∂y ∂y
∂T ∂T ∂ 2T
vx + vy =α 2 (6.1.10c)
∂x ∂y ∂y
La velocidad y la presión están sujetas a las siguientes condiciones de frontera:

En x = 0, vx = u∞ , T = T∞ (6.1.11a)
En y = 0, vx = vy = 0, T = T0 (6.1.11b)
En y = δH , vx = u∞ , T = T∞ (6.1.11c)

6.1.2 Solución aproximada de las ecuaciones de capa límite

Campo de velocidad
Siguiendo el método de von Kármán-Pohlhausen, se propone el siguiente perfil de velocidad de
forma polinomial;
y y2 y3
vx = a + b +c 2 +d 3 (6.1.12a)
δH δH δH
donde δH es una función de x que debe determinarse. La velocidad vx está sujeta a las siguientes
condiciones de frontera

En y = 0, vx = 0 (6.1.12b)
En y = δH , vx = u∞ (6.1.12c)

Para calcular todos los coeficientes del polinomio se necesitan dos condiciones adicionales. Dado
que fuera de la capa límite el transporte se da por esfuerzos inerciales y el arrastre viscoso es
prácticamente nulo, es razonable imponer que
∂ vx
En y = δH , =0 (6.1.12d)
∂y
Por último, dado que en y = 0, vx = vy = 0, la ecuación diferencial (6.1.10b) indica que:

∂ 2 vx
En y = 0, =0 (6.1.12e)
∂ y2
Usando estas condiciones de frontera se deduce que:
vx 3 y 1 y3
= − (6.1.13)
u∞ 2 δH 2 δH3
Para determinar la componente y de la velocidad se utiliza la ecuación de continuidad, la cual al ser
integrada semidefinidamente desde y = 0, da como resultado:
ζZ=y
∂ vx
vy − vy |y=0 = − dζ (6.1.14)
∂x
ζ =0
6.1 Ecuaciones de capa límite de energía y cantidad de movimiento para flujo
laminar 111
O bien, al llevar a cabo las operaciones de diferenciación e integración correspondientes y tomando
en cuenta que vy |y=0 = 0:

3u∞ y2 1 y4 dδH
 
vy = − (6.1.15)
4 δH2 2 δH4 dx
Para determinar el espesor de la capa límite se integra la ecuación de cantidad de movimiento desde
y = 0 hasta y = δH para obtener, al tomar en cuenta la condición de frontera de que el esfuerzo
cortante es nulo en la capa límite:
y=δ
Z H 
∂ vx ∂ vx ∂ vx 3νu∞
vx + vy dy = −ν =− (6.1.16)
∂x ∂y ∂y y=0 2δH
y=0

Sustituyendo vx y vy en esta ecuación y llevando a cabo las operaciones correspondientes resulta

que:
39 dδH 3ν
= (6.1.17)
280 dx 2δH u∞
La cual, puede integrarse, sujeta a la condición de frontera
En x = 0, δH = 0 (6.1.18)
para obtener
p
δH = 4.64 νx/u∞ (6.1.19)
Con este resultado se concluye el análisis de la capa límite de transferencia de cantidad de mo-
vimiento. Los perfiles de velocidad dados en las ecuaciones (6.1.13) y (6.1.15) se utilizan en los
siguientes párrafos para determinar el espesor de la capa límite térmica.
Campo de temperatura
Para la temperatura se procede de manera similar y se comienza con la propuesta del siguiente
polinomio:
y y2 y3
T = a+b +c 2 +d 3 (6.1.20a)
δT δT δT
Para estimar las cuatro constantes, se usan las siguientes condiciones de frontera:
∂ 2T
En y = 0, T = T0 ; =0 (6.1.20b)
∂ y2
∂T
En y = δT , T = T∞ ; =0 (6.1.20c)
∂y
de tal forma que la solución queda como:
1 y3
 
3 y
T = T0 + (T∞ − T0 ) − (6.1.21)
2 δT 2 δT3
Esta solución, junto con las otras obtenidas para las componentes x y y del vector velocidad se
sustituyen en la forma integrada de la ecuación de energía térmica:
y=δ
Z T 
∂T ∂T ∂T 3
vx + vy dy = −α = −α(T∞ − T0 ) (6.1.22)
∂x ∂y ∂y y=0 2δT
y=0
112 Capítulo 6. Capa límite

para obtener:
y=δ
Z T
3y y3 y3 dδT 3 y2 1 y4 y2 dδH
     
y 1 2α
− − 2+ 4 + − − dy = −
δH δH3 δT δT dx 2 δH2 2 δH4 δT δT3 dx δT u∞
y=0
(6.1.23)
O bien, definiendo a la razón entre ambas capas límite como

R = δT /δH (6.1.24)

la ecuación (6.1.23) se puede escribir como

y=δ
Z T
y3 y3 3 y2 1 y4 3 y2
    
3y y 1 dδT 2α
R − 3 R3 − 2+ 4 + 2
R− 4
R − 3 dy = −
δT δT δT δT 2 δT 2 δT δT δT dx δT u∞
y=0
(6.1.25)
Suponiendo que R es constante y llevando a cabo los pasos de integración correspondientes, se
llega al siguiente resultado:
R3
 
R dδT α
− δT = (6.1.26)
10 140 dx u∞
Esta ecuación puede resolverse sujeta a la condición de frontera:

En x = 0, δT = 0 (6.1.27)

De esta forma, el resultado es:

R3
 
R αx
− δT2 = (6.1.28)
20 280 u∞
O bien, en términos de la capa límite hidrodinámica, tomando en cuenta el resultado dado en la
ecuación (6.1.19),
R5 20αx 20α 0.93
R3 − = = = (6.1.29)
14 u∞ δH2 2
(4.64) ν Pr

en donde el número de Prandtl se define como Pr = ν/α. En este momento, es conveniente notar
los siguientes límites para esta ecuación:
1. Si Pr → ∞, entonces R → 0, esto es: δT /δH → 0.
2. Si Pr = 1, entonces R ≈ 1 y δT = δH .
3. Para varios gases, Pr = 0.7 y la solución de la ecuación es: δT ≈ 1.13δH .
Para este último caso, una aproximación de la solución es:

R = Pr−1/3 (6.1.30)

O bien,
r
νx −1/3
δT = 4.64 Pr (6.1.31)
u∞
Con estos resultados se puede calcular el flux interfacial de calor como:
r
∂T 3k u∞ 1/3
qy |y=0 = −k = (T0 − T∞ ) = 0.323k Pr (T0 − T∞ ) (6.1.32)
∂ y y=0 2δT νx
6.1 Ecuaciones de capa límite de energía y cantidad de movimiento para flujo
laminar 113
La cual, al comprarla con la ley de enfriamiento de Newton:

qy |y=0 = hloc (T0 − T∞ ) (6.1.33)

da como resultado
r
u∞ 1/3
hloc = 0.323k Pr (6.1.34)
νx
O bien, al promediar sobre toda la placa:
x=L
1 k
Z
hhloc i = hloc dx = 0.646 Re1/2 Pr1/3 (6.1.35)
L L
x=0

donde Re = u∞ L/ν. Definiendo el número de Nusselt como:

hhloc iL
Nu = (6.1.36)
k
Entonces, el resultado final de este análisis es:

Nu = 0.646Re1/2 Pr1/3 (6.1.37)

El cual está en excelente concordancia con el resultado exacto que es: Nu = 0.664Re1/2 Pr1/3 .
Ejercicio 6.1 — Desarrollos sin suponer R constante. En la obtención de los resultados
anteriores se supuso que las razón R entre las capas límite térmica e hidrodinámica es constante.
En este problema se busca obtener un resultado que no esté limitado por esta aproximación.
Para ello, sustituya de nuevo la ecuación (6.1.24) en la ecuación (6.1.23) para obtener
y=δ
Z T
y3 y3 dδT
 
3y y
R − 3 R3 − 2+ 4
δT δT δT δT dx
y=0

3 y2 1 y4 3 y2
   
1 dδT δT dR
+ R− R − − dy
2 δT2 2 δT4 δT δT3 dx R dx
2α
=− (6.1.38)
δT u∞
Lleve a cabo los pasos de integración necesarios para recuperar el siguiente resultado

R3 R2
   
R dδT 3 1 dR α
− δT + − δT2 = (6.1.39)
10 140 dx 2 15 70 dx u∞

La cual es una generalización de la ecuación (6.1.26). Como no se conoce δT , es conveniente

re-escribir esta ecuación en términos de R y δH , ya que esta última si es conocida. Lleve a cabo
las operaciones algebraicas correspondientes para obtener la siguiente expresión

R5
 
3 2 4 4 dR 13
R − + 4R − R x = (6.1.40)
14 7 dx 14Pr

Para resolver esta ecuación suponga que R5 /14  R3 y que R4 /7  R2 , para que se reduzca a

dR 13
R3 + 4R2 x = (6.1.41)
dx 14Pr
 114 Capítulo 6. Capa límite

Muestre que la solución general de este problema es (consejo: use y = R3 para llevar a cabo la
solución)
13
R3 = Cx−3/4 + (6.1.42)
14Pr
Para calcular a la constante C suponga que la capa límite térmica se desarrolla a partir de un
cierto espesor ξ , por lo tanto la condición de frontera a usar es

En x = ξ , R=0 (6.1.43)

De esta forma, la solución para R es

"  3/4 #1/3
Pr−1/3 ξ
R= 1− (6.1.44)
1.026 x

La cual se reduce a la ecuación (6.1.30) para el caso en el que ξ = 0. Utilice este resultado para
deducir que la forma local del coeficiente interfacial de transferencia de calor es

0.323k uνx∞ Pr1/3

p
hloc =  1/3 (6.1.45)
1 − (ξ /x)3/4


6.2 Capa límite cuando se especifica el flux de calor en la superficie

Reconsidere el problema anterior solo que ahora se conoce el flux de calor en la superficie y es
constante, es decir,

∂T
En x = 0, −k = q0 (6.2.1)
∂x
El resto de las condiciones de frontera no se alteran, por lo que el perfil polinomial de temperatura
satisface la siguiente expresión (ver ejercicio 6.2)

y3
 
2 δT q0 q0
T − T∞ = − y− 2 (6.2.2)
3 k k 3δT

De donde se obtienen las siguientes derivadas parciales

2 2y3 dδT
 
∂T q0
= − (6.2.3a)
∂x k 3 3δT3 dx
y2
 
∂T q0
=− 1− 2 (6.2.3b)
∂y k δT

Tomando en cuenta que las expresiones para las componentes de la velocidad siguen siendo

1 y3 3
 
3 y
vx = u∞ R− R (6.2.4a)
2 δT 2 δT3
3u∞ y2 1 y4 3 dδT
 
vy = R− R (6.2.4b)
4 δT2 2 δT4 dx
6.2 Capa límite cuando se especifica el flux de calor en la superficie 115

Note que aquí se utilizó la relación δH = δT /R. Sustituyendo las ecuaciones (6.2.3) y (6.2.4) en la
forma integrada (en la dirección y) de la ecuación de energía térmica

y=δ
Z T 
∂T ∂T ∂T
vx + vy dy = −α (6.2.5)
∂x ∂y ∂y y=0
y=0

da como resultado la siguiente expresión

y=δ
Z T
y3 2 2y3 y2 y4 3 y2
    
dδT 3 y 3 α
R − 3 R3 − − R − R 1− 2 dy =
dx 2 δT 2δT 3 3δT3 4 δT2 2δT4 δT u∞
y=0
(6.2.6)

La cual se puede reducir al llevar a cabo los pasos de integración correspondientes para dar lugar a
la siguiente ecuación diferencial ordinaria

R R3
 
dδT α
− δT = (6.2.7)
5 70 dx u∞

La solución de esta ecuación, tomando en cuenta que en x = 0, δT = 0, es

R R3 δT2 αx
 
− = (6.2.8)
5 70 2 u∞

Sustituyendo en esta expresión δT = RδH da lugar a

R3 R5 1
− = (6.2.9)
10 140 (4.64)2 Pr

R5 R3
Para situaciones en las que R < 1, es razonable suponer que 140  10 y por lo tanto se obtiene que

R = 0.774Pr−1/3 (6.2.10)

Con este resultado y tomando en cuenta a la ecuación (6.1.19) se deduce que el espesor de la capa
límite térmica está dado por
−1/2
δT = 3.59xPr−1/3 Rex (6.2.11)

donde Rex = u∞ x/ν. Con este resultado se puede determinar la temperatura en y = 0 como

2.4q0 x −1/3 −1/2

T (y = 0) − T∞ = Pr Rex (6.2.12)
k
De donde puede despejarse el flux de calor q0 como sigue

0.417k 1/3 1/2

q0 = Pr Rex (T (y = 0) − T∞ ) (6.2.13)
x
Haciendo el símil con la ley de enfriamiento de Newton, se deduce que

hloc x 1/2
Nuloc = = 0.417Pr1/3 Rex (6.2.14)
k
116 Capítulo 6. Capa límite

Por último, definiendo al coeficiente promedio de transferencia de calor como

x=L
1
Z
hhi = hloc dx (6.2.15)
L
x=0

se obtiene el siguiente resultado

0.835k 1/3 1/2

hhi = Pr Re (6.2.16)
L
Con el cual se deduce que el número de Nusselt promedio es

Nu = 0.835Pr1/3 Re1/2 (6.2.17)

Note que este valor es aproximadamente 30 % superior al que se dedujo para el caso en el que la
temperatura es constante en la superficie.

Ejercicio 6.2 — Perfil de temperatura. Partiendo del polinomio

y y2 y3
T = a+b +c 2 +d 3 (6.2.18a)
δT δT δT

sujeto a las siguientes condiciones de frontera:

∂T ∂ 2T
En y = 0, −k = q0 ; =0 (6.2.18b)
∂x ∂ y2
∂T
En y = δT , T = T∞ ; =0 (6.2.18c)
∂y
Demuestre que se obtiene el perfil de temperatura dado en la ecuación (6.2.2). 
 7. Modelos macroscópicos

7.1 Introducción
Con este teme se cierra el curso y se fundamenta en la forma macroscópica de la ecuación de
energía térmica

d
Z Z Z
ρCp T dV + n · [ρCp T (v − w) + q] dA = Φ dV (7.1.1)
dt
VA AA VA

La clave está en reconocer las suposiciones y condiciones de frontera correspondientes que permitan
la aplicación de esta ecuación. Se aborda una variedad de problemas de interés en ingeniería para
ilustrar el uso de esta ecuación.

7.2 Transferencia de calor en un panel fotovoltaico

Considere un solo panel fotovoltaico como el que se muestra en la Figura 7.1, dentro del panel
no hay convección ni radiación por lo que la ecuación de energía térmica se reduce a:

∂T
ρc p = ∇ · (k∇T ) (7.2.1)
∂t
La cual está sujeta a las siguientes condiciones de frontera:

nnc · q = 0, en Anc (7.2.2)

nc · q = − nc · qr + h(T − Tamb ) , en Ac (7.2.3)
| {z } | {z } | {z }
conducción radiación pérdidas por convección

En donde Ac es la superficie expuesta del panel y Anc las superficies que no tienen contacto directo
con la radiación solar. Como puede notarse, en la superficie expuesta ingresa calor por radiación
y al mismo tiempo se disipa calor por convección al ambiente, mientras que en las superficies no
expuestas se supone que no hay transferencia de calor.
Breviario sobre modelado Francisco J. Valdés Parada
Septiembre, 2012

2.4.4. Respuesta transitoria de un colector solar

118 Capítulo 7. Modelos macroscópicos
El uso de equipos de almacenamiento y aprovechamiento de energía solar han ido en
Para el caso bajo estudio no hay convección ni generación de calor dentro del panel fotovoltaico.
aumento enecuación
Por ello, la las últimas décadas.
macroscópica Su usotérmica,
de energía va desde relojes
para una regiónhasta
fija, seprototipos
reduce a de casas
d
Z Z
enterasρC
que
p pueden + n · q dA
T dVfuncionar =0
mediante este tipo de energía. (7.2.4)
dt
Vf Af

Los colectores
donde de energía
ρC se supuso constante.solar (Figura
Se desea 2.4-3)
encontrar funcionan,
un modelo casilaen
que describa su totalidad,
temperatura bajo
promedio
p
en el panel como función del tiempo. Para ello, se define el siguiente operador de promediado
condiciones
volumétrico transitorias en particular en los periodos de arranque durante las primeras
1
Z
horas hT
dei =
la mañana.
T dV El objetivo es determinar la dependencia temporal(7.2.5)
de la
V
V
temperatura de la superficie de un colector solar a lo largo de un año en el Valle de
donde V denota el dominio ocupado por el panel fotovoltaico (el cual no cambia con el tiempo) y
V su volumen. Dividiendo ambos lados de la ecuación (7.2.4) entre V , se obtiene que
México.
dhT i 1
Z
ρCp dV + n · q dA = 0 (7.2.6)
El sistema consiste
dt en
V una celda solar, dado que no hay convección ni generación de
AA

calor al interior
El primer términodeyaella,
está la
en ecuación de energía
su forma final, térmica
por lo que [Ec.se(1.5.5)]
la atención seelreduce
centra en a término,
segundo
el cual puede descomponerse como sigue:
1 1 1ρ C ∂T = k ∇ 2T (2.4.40)
Z Z Z
n · q dA = nc · q dA + {
p nnc · q dA (7.2.7)
V V V14 2 4∂3t transporte
AA Ac A nc difusivo
acumulación

Considerando las condiciones de frontera, la ecuación anterior toma la forma siguiente:

Evidentemente la temperatura T es función
Z de la posición y el tiempo, por lo que para
1 1 1
Z Z
n · q dA = − nc · qr dA + h(T − Tamb ) dA (7.2.8)
V V V
simplificar
AA el modelo se A
buscará
c
desarrollar
Ac una expresión en el cual la temperatura sólo
Para expresar este resultado en una forma más conveniente se define la potencia promedio incidente
dependa del tiempo. Para ello se define la temperatura promedio en el volumen como,
al panel solar como:
1 1
Z

V V∫
S=
Ac
nc · qr dA T = T dV (2.4.41)
(7.2.9)
Ac

Figura
Figura 2.4-3: Colector
7.1: Esquema delsolar
panelinstalado en la
fotovoltaico azotea en
instalado dellaedificio B de
azotea del la UAM-Iztapalapa.
edificio B de la UAM-
Iztapalapa.

84
7.2 Transferencia de calor en un panel fotovoltaico 119

Además, en la superficie del panel se descompone la temperatura en su valor promedio y desviacio-

nes por lo que:
1 1 1
Z Z Z
h(T − Tamb ) dA = h(hT i − Tamb ) dA + hT̃ dA (7.2.10)
V V V
Ac Ac Ac

Para simplificar este resultado se supone que

1
Z
T̃ dA  hT i (7.2.11)
Ac
Ac

y que el coeficiente interfacial de transferencia de calor, h, es constante. De esta forma, la ecuación

(7.2.8) se escribe como
1
Z
n · q dA = −ac S + ac h(hT i − Tamb ) (7.2.12)
V
AA

donde ac ≡ Ac /V representa el área interfacial por unidad de volumen. Bajo estas condiciones, el
modelo macroscópico es:
dhT i
ρc p = ac S − hac (hT i − Tamb ) (7.2.13)
dt
Para describir a la función S se pueden utilizar datos experimentales o bien modelos como la
ecuación de Stefan-Boltzman. En el caso del valle de México, S oscila a lo largo del año pero una
aproximación razonable es de 20 × 106 J/m2 día.

Ejercicio 7.1 — Sobre las desviaciones de la temperatura. Demuestre que, si se acepta la

siguiente suposición,
1
Z
T dA ≈ hT i
Ac
Ac

entonces puede concluirse que

1
Z
T̃ dA  hT i
Ac
Ac

Ejercicio 7.2 En un día típico en la CDMX, es normal que la radiación térmica obedezca a una
relación polinomial

S = a + bt + ct 2 + dt 3 (7.2.14)

La cual obedece a las siguientes restricciones temporales

Cuando t = 0, S=0 (7.2.15)

Cuando t = tmáx /2, S = Smáx (7.2.16)
dS
Cuando t = tmáx /2, =0 (7.2.17)
dt
Cuando t = tmáx , S = 0 (7.2.18)
 120 Capítulo 7. Modelos macroscópicos

Figura 7.2: Esquema de un reactor tipo tanque agitado

Utilice estas restricciones para demostrar que S puede expresarse como

t2
 
t
S = 4Smáx − 2 (7.2.19)
tmáx tmáx

Bajo esta premisa, puede deducirse que

t2
 
t
Ta = Ta mı́n + 4Ta máx − 2 (7.2.20)
tmáx tmáx

donde Ta mı́n y Ta máx representan los valores mínimo y máximo de la temperatura ambiente.
Sustituya las ecuaciones (7.2.19) y (7.2.20) en la ecuación (7.2.13) y resuélvala sujeta a la
condición inicial

Cuando t = 0, hT i = T0 (7.2.21)

7.3 Tanque agitado homogéneo

7.3.1 Planteamiento del problema
Considere un tanque agitado como el mostrado en la Figura 7.2 en donde ingresa una corriente
de fluido a una temperatura Ten y velocidad ven constantes. Suponga que dentro del tanque tiene
lugar una reacción química homogénea, la cual constituye una fuente de generación de calor que
puede modelarse como φ = RA ∆H. Siendo RA y ∆H la tasa de reacción y la entalpía de reacción,
respectivamente. Desarrolle una ecuación que describa en forma macroscópica la dinámica de la
temperatura en el tanque.
El sistema es el líquido contenido en el tanque agitado. Suponiendo que el flujo de entrada es
igual al flujo de salida, es posible considerar una región fija para el análisis. Para las condiciones
que especifica el problema, la forma macroscópica de la ecuación de energía térmica se expresa
como:
dhT i 1
Z
ρc p + n · (ρc p vT − k∇T ) dA = hRA ∆Hi (7.3.1)
dt V
A

En donde la temperatura promedio se define como

1
Z
hT i = T dV (7.3.2)
V
V
7.3 Tanque agitado homogéneo 121

Las superficies (A ) con las que está en contacto el fluido son las siguientes: superficies de entrada
Ae , de salida As , paredes A p y agitador Aa . Suponiendo que las paredes del tanque están aisladas
y que el transporte de calor a las entradas, salidas y en el agitador se da predominantemente por
convección, se tienen las siguientes condiciones de frontera:

En Ae , T = Ten , v = −ven , ne , ne · k∇T = 0 (7.3.3)

En As , T = Tsal , v = vsal ns , ns · k∇T = 0 (7.3.4)
En Aa , na · (ρc p vT − k∇T ) = 0 (7.3.5)
En A p , n p · (ρc p vT − k∇T ) = 0 (7.3.6)

La integral de superficie puede descomponerse en cuatro integrales para cada superficie. Usando
las condiciones de frontera se obtiene lo siguiente

dhT i 1 1 1
Z Z
− Ten ven dA + Tsal vsal dA = hRA i∆H (7.3.7)
dt V V ρc p
Ae Ae

Aquí se consideró que la entalpía de reacción puede suponerse como constante. Para avanzar en los
desarrollos se adoptan las siguientes suposiciones:
1. La temperatura de entrada es constante.
2. La temperatura de salida es aproximadamente igual a la temperatura promedio en el tanque.
Esta suposición es razonable bajo condiciones de mezclado perfecto.
Bajo estas condiciones el modelo promedio puede escribirse como:

dhT i Ten hT i 1
Z Z
− ven dA + vsal dA = hRA i∆H (7.3.8)
dt V V ρc p
Ae Ae

Dado que el flujo volumétrico de entrada se supuso igual al de salida se tiene que
Z Z
ven dA = vsal dA ≡ Q̇ (7.3.9)
Ae As

donde Q̇ es el flujo volumétrico; además se tomó en cuenta la suposición de flujo incompresible.

Usando este resultado, la ecuación macroscópica de energía térmica toma la forma:

dhT i 1 hRA i∆H

= (Ten − hT i) + (7.3.10)
dt tR ρc p

donde tR = V /Q̇ es el tiempo de residencia en el reactor.

Ejercicio 7.3 — Solución del modelo. Resuelva la ecuación diferencial (7.3.10) para el caso
hRA i∆H
en el que ρc p sea constante considerando las siguiente función de entrada:

2
Ten (t) = Tin [1 + 0.5sen(wt)] (7.3.11)
3
en donde Tin y w representan la temperatura máxima de entrada y la frecuencia de la función de
entrada, respectivamente.

 122 Capítulo 7. Modelos macroscópicos

7.4 Tanque agitado no homogéneo

Reconsidere el problema anterior con la diferencia de que ahora en el tanque agitado se
encuentra una resistencia eléctrica que se encuentra a una temperatura Ts ; además dentro del tanque
hay partículas esféricas a las cuales el fluido les cede el calor que gana de la resistencia eléctrica.
Desarrolle una ecuación macroscópica de la temperatura del fluido en el tanque.
El planteamiento del problema no es muy distinto del anterior, la principal diferencia es que no
hay una fuente interna de generación de calor, por lo que la ecuación macroscópica para el fluido
en el tanque en la escala microscópica es:

dhT i 1
Z
ρc p + n · (ρc p vT − k∇T ) dA = 0 (7.4.1a)
dt V
A

En este caso las superficies con las que tiene contacto el fluido son: entradas (Ae ), salidas (As ),
paredes (A p ), agitador (Aa ), esferas (Aσ ) y resistencia eléctrica (Ar ). Las condiciones de frontera
en cada una de estas superficies son:

En Ae , T = Ten , v = −ven , ne , ne · k∇T = 0 (7.4.1b)

En As , T = Tsal , v = vsal ns , ns · k∇T = 0 (7.4.1c)
En Aa , n · (ρc p vT − k∇T ) = 0 (7.4.1d)
En A p , n p · k∇T = 0, v = 0 (7.4.1e)
En Aσ , nσ · (ρc p vT − k∇T ) = −nσ · kσ ∇Tσ (7.4.1f)
En Aσ , nσ · kσ ∇Tσ = hσ ( Tσ |r=R − hT i) (7.4.1g)
En Ar , v = 0, −nr · k∇T = hr (hT i − Tr ) (7.4.1h)

Utilizando las condiciones de frontera, la ecuación (7.4.1a) se transforma en:

 
dhT i ρc p 
Z Z
ρc p − Ten ven dA − Tsal vsal dA
dt V
Ae As
hσ Aσ hr Ar
= ( Tσ |r=R − hT i) + (Tr − hT i) (7.4.2)
V V
Suponiendo que la temperatura de entrada es constante y que la temperatura de salida es igual a la
temperatura promedio en el tanque se tiene que:

dhT i ρc p hσ Aσ hr Ar
ρc p = (Ten − hT i) + ( Tσ |r=R − hT i) + (Tr − hT i) (7.4.3)
dt tR V V
Como puede notarse, esta ecuación está acoplada a la temperatura en las partículas, las cuales se
han supuesto que son idénticas en tamaño y forma. La ecuación de energía térmica en cada partícula
es:
∂ Tσ
ρσ c pσ = ∇ · (kσ ∇Tσ ) (7.4.4)
∂t

Ejercicio 7.4 — Distribución de tamaños de partícula. En este ejercicio se hizo la supo-

sición de que todas las partículas eran idénticas en tamaño; muestre cómo se modificaría el
modelo si existiera una distribución de N tamaños de partícula. 
 7.5 Transferencia de calor en un lecho fijo 123

Ejercicio 7.5 — Calentamiento no uniforme pt. 1. Reconsidere el ejercicio suponiendo

ahora que la temperatura en la resistencia no es constante, sino que cambia linealmente con la
posición. ¿Cómo se modifica el modelo final? 

Ejercicio 7.6 — Calentamiento no uniforme pt. 2. Muestre cómo se modificaría el modelo

macroscópico si ahora se especifica el flux de calor en la superficie de la resistencia eléctrica
por lo que la temperatura en dicha superficie no es constante. 

7.5 Transferencia de calor en un lecho fijo

Considere ahora la transferencia de calor en un tubo de radio R y largo Lpor el cual circula un
fluido (fase-β ) alrededor de partículas esféricas (fase-σ ) idénticas de radio R, las cuales tienen una
fuente interna de generación de calor. Desarrolle un modelo macroscópico que describa la dinámica
de la transferencia de calor en el fluido a lo largo del tubo.
Para analizar el sistema, considere una región de promediado V de forma cilíndrica, cuyo radio
es R y su altura es r0 . Esta altura debe ser tal que contenga suficientes partículas de tal manera
que no importa donde se coloque su centro geométrico, el resultado sea el mismo, al menos en
cuanto a fracción volumétrica de las fases involucradas. Una forma de asegurar esto es imponiendo
la siguiente restricción de escala

máx(`β , `σ )  r0  L (7.5.1)

En donde `β y `s igma representan las longitudes características de las fases fluida y sólida en el
nivel de escala del continuo. De esta forma, V = Vβ + Vσ , donde Vi (i = β , σ ) son independientes
de la dirección axial.
La ecuación de energía térmica en el fluido en el tubo es
d
Z Z h i
(ρCp )β Tβ dV + nβ · (ρCp )β Tβ (vβ − wβ ) + qβ dA = 0 (7.5.2)
dt
Vβ Aβ

Dado que Vβ es constante, la ecuación anterior puede adaptarse a una región fija como sigue

∂ hTβ iβ 1
Z  
(ρCp )β + nβ · (ρCp )β Tβ vβ + qβ dA = 0 (7.5.3)
∂t Vβ
Aβ

En donde se usó la siguiente definición de la temperatura promedio

1
Z
hTβ iβ = Tβ dV (7.5.4)
Vβ
Vβ

La superficie Aβ puede descomponerse en las superficies de entrada y salida (Aβ es ), la superficie

de la pared del tubo (Aβ p ) y la interfase con las partículas sólidas (Aβ σ ). La velocidad del fluido
en las últimas dos interfases es nula, por lo que la ecuación (7.5.3) puede escribirse como sigue

∂ hTβ iβ 1
Z   1
Z
(ρCp )β + nβ es · (ρCp )β Tβ vβ + qβ dA + nβ σ · qβ dA = 0 (7.5.5)
∂t Vβ Vβ
Aβ es Aβ σ

Aquí, se supuso que las paredes del tubo están aisladas y por ello

En Aβ p , nβ p · kβ ∇Tβ = 0 (7.5.6)
124 Capítulo 7. Modelos macroscópicos

En este punto, es conveniente considerar la siguiente expresión que forma parte del Teorema del
promediado volumétrico (ver ecuación (8) en Slattery, 1967) para un tensor B
Z Z
∇· B dV = n · B dA (7.5.7)
Vβ Aβ es

Este resultado puede extenderse a un vector a como sigue

Z Z
∇· a dV = n · a dA (7.5.8)
Vβ Aβ es

Usando este teorema en la ecuación (7.5.5) se obtiene que

∂ hTβ iβ D Eβ 1
Z
(ρCp )β + ∇ · (ρCp )β Tβ vβ + qβ + nβ σ · qβ dA = 0 (7.5.9)
∂t Vβ
Aβ σ

Antes de continuar es conveniente notar la siguiente identidad, la cual está implícita en los desarro-
llos hechos hasta el momento,
Teorema 7.5.1 — Teorema del promediado espacial para un vector.

N
1 1
Z Z
nβ · a dA = h∇ · aiβ = ∇ · haiβ + ∑ nβ j · a dA (7.5.10)
Vβ Vβ j=1
Aβ Aβ j

En donde N representa el número de fases con las que está en contacto la fase-β en la región de
promediado, sin incluir las superficies de entrada y salida. La ecuación anterior se conoce como
el teorema del promediado espacial y puede expresarse para el gradiente de un escalar ψ como
sigue

Teorema 7.5.2 — Teorema del promediado espacial para un escalar.

N
1
Z
h∇ψiβ = ∇hψiβ + ∑ nβ j ψ dA (7.5.11)
Vβ j=1
Aβ j

Usando la ley de Fourier, la ecuación (7.5.9) se escribe como

∂ hTβ iβ β β 1
Z
(ρCp )β + (ρCp )β ∇ · Tβ vβ = kβ ∇ · ∇Tβ + nβ σ · kβ ∇Tβ dA (7.5.12)
∂t Vβ
Aβ σ

Note que se despreciaron los cambios espaciales de (ρCp )β y kβ . Usando la forma escalar del
teorema del promediado espacial al primer término en el lado derecho de la ecuación anterior se
obtiene que:

∂ hTβ iβ β
(ρCp )β + (ρCp )β ∇ · Tβ vβ
∂t  
1  1
Z Z
2 β
= kβ ∇ Tβ + kβ ∇ ·  nβ σ Tβ dA + nβ σ · kβ ∇Tβ dA (7.5.13)
Vβ Vβ
Aβ σ Aβ σ
7.5 Transferencia de calor en un lecho fijo 125

Aquí se consideró que no hay un gradiente de temperatura en la interfase fluido-pared. En este

punto es conveniente descomponer a la velocidad del fluido de acuerdo a (Gray, 1975)

vβ = hvβ iβ + ṽβ (7.5.14)

Sustituyendo esta descomposición en la ecuación (7.5.13) se obtiene que

∂ hTβ iβ   1
Z
= kβ∗ ∇2
β β β
(ρCp )β + (ρCp )β ∇ · Tβ vβ Tβ + nβ σ · kβ ∇Tβ dA (7.5.15)
∂t Vβ
Aβ σ

En donde kβ∗ es el coeficiente de dispersión total de calor en la fase fluida y se define como

∗ β β kβ Z β
k ∇ Tβ = kβ ∇ Tβ + nβ σ Tβ dA − (ρCp )β ∇ · Tβ ṽβ (7.5.16)
Vβ
Aβ σ

La ecuación (7.5.15) casi está completa, ya que todos los términos excepto el último en el lado
derecho están escritos en términos de la temperatura promedio. Para atender esto, se consideran las
siguientes condiciones de frontera interfaciales

En Aβ σ , −nβ σ kβ ∇Tβ = −nβ σ kσ ∇Tσ (7.5.17a)

En Aβ σ , nβ σ kσ ∇Tσ = h(Tσ − Tβ ) (7.5.17b)

De esta forma, la ecuación (7.5.15) puede expresarse como sigue

∂ hTβ iβ   1
Z
= kβ∗ ∇2
β β β
(ρCp )β + (ρCp )β ∇ · Tβ vβ Tβ + h(Tσ − Tβ ) dA (7.5.18)
∂t Vβ
Aβ σ

Para cerrar el modelo es necesario contar con información sobre las partículas, las cuales obedecen
a la siguiente forma de la ecuación de energía térmica
∂ hTσ iσ 1
Z
(ρCp )σ + nβ · qσ dA = hΦσ iσ (7.5.19)
∂t Vσ
Aσ

Donde Φσ es conocida, por lo que hΦσ iσ se puede considerar como un dato. Note que en la
ecuación anterior no se consideró la convección ya que no la hay dentro de las partículas del lecho.
Usando la identidad dada en la ecuación (7.5.10) para el término conductivo se obtiene que
∂ hTσ iσ 1
Z
(ρCp )σ + ∇ · hqσ iσ − nβ σ · qσ dA = hΦσ iσ (7.5.20)
∂t Vσ
Aβ σ

O bien, usando la ley de Fourier,

∂ hTσ iσ 1
Z
(ρCp )σ = ∇ · (kσ h∇Tσ iσ ) − nβ σ · kσ ∇Tσ dA + hΦσ iσ (7.5.21)
∂t Vσ
Aβ σ

Aplicando la forma del teorema del promediado dada en la ecuación (7.5.11) en el término conduc-
tivo se obtiene
∂ hTσ iσ 1
Z
(ρCp )σ = kσ e f ∇2 hTσ iσ − nβ σ · kσ ∇Tσ dA + hΦσ iσ (7.5.22)
∂t Vσ
Aβ σ
126 Capítulo 7. Modelos macroscópicos

En donde kσ e f es el coeficiente de conductividad térmica efectiva y se define como

kσ
Z
σ σ
kσ e f ∇hTσ i = kσ ∇hTσ i + nσ β Tσ dA (7.5.23)
Vσ
Aβ σ

Sustituyendo la condición de frontera interfacial dada en la ecuación (7.5.17b) en el término integral

de la ecuación (7.5.22) da como resultado
∂ hTσ iσ 1
Z
= kσ e f ∇2 hTσ iσ −


(ρCp )σ h Tσ − Tβ dA + hΦσ iσ (7.5.24)
∂t Vσ
Aβ σ

Para avanzar en los desarrollos es conveniente descomponer las temperaturas en ambas fases como

T j = hT j i j + T̃ j , j = β,σ (7.5.25)

De esta forma, se tiene que

Z
  Z

h Tσ − Tβ dA = Aβ σ h hTσ iσ − hTβ iβ + h T̃σ − T̃β dA (7.5.26)
Aβ σ Aβ σ

Entonces, para condiciones en las cuales sea razonable suponer que

Z
 
h T̃σ − T̃β dA  Aβ σ h hTσ iσ − hTβ iβ (7.5.27)
Aβ σ

las formas cerradas de las ecuaciones de transporte de calor en ambas fases son

Modelo de dos ecuaciones

∂ hTβ iβ   Aβ σ  
= kβ∗ ∇2 Tβ
β β β
(ρCp )β + (ρCp )β ∇ · Tβ vβ + h hTσ iσ − hTβ iβ
∂t Vβ
(7.5.28a)
∂ hTσ iσ Aβ σ  
(ρCp )σ = kσ e f ∇2 hTσ iσ − h hTσ iσ − hTβ iβ + hΦσ iσ (7.5.28b)
∂t Vσ

Este conjunto de ecuaciones representa el modelo cerrado de dos ecuaciones para modelar el sistema
y es adecuado cuando el interés está en la temperatura promedio del fluido o de las partículas en el
sistema. Debe quedar claro que en las ecuaciones anteriores el operador ∇ sólo tiene la componente
axial.
Una alternativa al modelo anterior es el modelo de equilibrio local térmico, el cual se basa en
suponer que es razonable proponer lo siguiente
1
Z
hTβ iβ ≈ hTσ iσ ≈ hT i = T dV (7.5.29)
V
V

De esta forma, multiplicando a la ecuación (7.5.28a) por εβ = Vβ /V y a la ecuación (7.5.28b) por

εσ = Vσ /V y sumando las ecuaciones resultantes se obtiene

Modelo de equilibrio local térmico

∂ hT i 
β

hρCp i + εβ (ρCp )β ∇ · hT i vβ = hki∇2 hT i + εσ hΦσ iσ (7.5.30)
∂t
7.5 Transferencia de calor en un lecho fijo 127

En donde se usaron las siguientes definiciones

hρCp i = εβ (ρCp )β + εσ (ρCp )σ (7.5.31a)

hki = εβ kβ∗ + εσ kσ e f (7.5.31b)

Ejercicio 7.7 — Modelo de no-equilibrio y de una ecuación. El modelo dado en la ecua-

ción (7.5.30) es atractivo desde el punto de vista que en una sola ecuación se abarcan todos
los procesos de transporte de calor que ocurren en la región de promediado. Sin embargo, está
limitado por la aproximación dada en la ecuación (7.5.29) que corresponde al equilibrio local
térmico. Para atender esta situación, en este problema se busca desarrollar un modelo de una
sola ecuación que no involucre esta suposición. Para ello, comience con la ecuación de energía
térmica aplicada a la región V , la cual al dividirla entre V da como resultado

1 ∂ 1
Z Z
(ρCp ) T dV + n · [(ρCp ) T v + q] dA = εσ hΦσ iσ (7.5.32)
V ∂t V
V A

Aquí se consideró que sólo hay una fuente de calor en la fase sólida. Muestre que el primer
término del lado izquierdo de la ecuación anterior puede descomponerse como

1 ∂
Z
∂ hTβ iβ ∂ hTσ iσ
(ρCp ) T dV = εβ (ρCp )β + εσ (ρCp )σ (7.5.33)
V ∂t ∂t ∂t
V

Más aún, el término convectivo sólo tiene lugar en la fase fluida, por lo que debe demostrar que
la siguiente ecuación es válida
1
Z
β
n · (ρCp ) T v dA = εβ (ρCp )β ∇ · Tβ vβ (7.5.34)
V
A

Mientras que el término conductivo se debe demostrar que se expresa como

1
Z
n · q dA = εβ ∇ · hqβ iβ + εσ ∇ · hqσ iσ (7.5.35)
V
A

Note que aquí ya se consideró la condición de frontera de continuidad del flux de calor interfacial
y el aislamiento de la pared. Utilice la ley de Fourier en cada fase y el teorema del promediado
para obtener que
 
1 1
Z Z
n · q dA = −∇ · kβ εβ ∇hTβ iβ + kσ εσ ∇hTσ iσ + nβ σ (kβ Tβ − kσ Tσ ) dA (7.5.36)
 
V V
A Aβ σ

Sustituya estos resultados, tomando en cuenta la descomposición de la velocidad dada en la

 128 Capítulo 7. Modelos macroscópicos

ecuación (7.5.14), en la ecuación (7.5.32) para obtener

∂ hTβ iβ ∂ hTσ iσ 
β β

εβ (ρCp )β + εσ (ρCp )σ + εβ (ρCp )β ∇ · Tβ vβ
 ∂t ∂t 
1
Z
β
= ∇· kβ εβ ∇hTβ iβ + kσ εσ ∇hTσ iσ + nβ σ (kβ Tβ − kσ Tσ ) dA − εβ (ρCp )β Tβ ṽβ

V

Aβ σ

+ εσ hΦσ iσ (7.5.37)

Para avanzar descomponga las temperaturas de cada fase de acuerdo a

T j = hT i + T̂ j , j = β,σ (7.5.38)

para obtener la siguiente expresión

∂ hT i 
β

hρCp i + εβ (ρCp )β ∇ · hT i vβ = ∇ · (kne ∇hT i) + εσ hΦσ iσ (7.5.39)
∂t
Aquí se supuso que

∂ hT̂β iβ ∂ hT̂σ iσ ∂ hT i
εβ (ρCp )β + εσ (ρCp )σ  hρCp i (7.5.40)
∂t ∂t ∂t
y además se definió

1
Z
kne ∇hT i = (kβ εβ +kσ εσ )∇hT i+kβ εβ ∇hT̂β iβ +kσ εσ ∇hT̂σ iσ + nβ σ (kβ Tβ −kσ Tσ ) dA
V
Aβ σ
 
β β β
− εβ (ρCp )β Tβ ṽβ + T̂β vβ (7.5.41)

La ecuación (7.5.39) es el modelo resultante y guarda la misma estructura matemática que el

modelo de equilibrio local térmico dado en la ecuación (7.5.30); sin embargo, la definición de
kne dista mucho de la correspondiente para hki. 

7.6 Transferencia de calor en intercambiadores de calor de doble tubo

Considere un sistema como el dibujado en la figura 7.3, donde se tienen dos corrientes fluidas
transportándose en las secciones interna y externa de un intercambiador de calor de doble tubo.
Suponiendo condiciones de estado estacionario, las ecuaciones que gobiernan el transporte de calor
microscópicamente en cada corriente son

∂ 2 TI
 
∂ TI kI ∂ ∂ TI
ρI CpI vzI = r + kI 2 (7.6.1a)
∂z r ∂r ∂r ∂z
∂ 2 TII
 
∂ TII kII ∂ ∂ TII
±ρII CpII vzII = r + kII (7.6.1b)
∂z r ∂r ∂r ∂ z2

Aquí se supuso que la única componente relevante del vector de velocidad es la componente
z, la cual es función de r, además se supuso simetría axial, por lo que no se consideró el término
de conducción angular. Además, en la ecuación (7.6.1b), el signo + corresponde a una operación
co-corriente y el signo − a la operación en contracorriente.
7.6 Transferencia de calor en intercambiadores de calor de doble tubo 129

De los cursos de mecánica de fluidos, se pueden deducir las siguientes expresiones para los
perfiles de velocidad:
"  2 #
I ξ
vzI = hvzI i 2 1 − (7.6.2a)
ξi
| {z }
fI (ξ )

" #
2 − 1) ln ξ /ξ + ξ 2 − ξ 2 ln ξ


(ξe e e e
vzII = hvzII iII (7.6.2b)
(1 + ξe2 ) ln ξe + 1 − ξe2
| {z }
fII (ξ )

donde se usaron las siguientes definiciones

r ri re
ξ= ; ξi = ; ξe = (7.6.3)
Ri Ri Ri
Así como los siguientes operadores de promediado
r=r
Z i
I 2
hvzI i = 2 vzI rdr (7.6.4a)
ri
r=0
r=R
Z i
2
II
hvzII i = 2 vzII rdr (7.6.4b)
Ri − re2
r=re

La condición de frontera interfacial es:

 
∂ TI ∂ TII hTI iI − hTII iII  
−kI ri = −kII re =   = reUe hTI iI − hTII iII (7.6.5)
∂r ∂r 1 1 re 1
re he + k ln ri + ri hi
r=ri r=re

Figura 7.3: Intercambiador de calor de doble tubo con acercamiento en la sección transversal.
130 Capítulo 7. Modelos macroscópicos

Dirigiendo la atención por el momento al fluido del tubo interno (región-I), se aplica el operador
de promediado definido en la ecuación (7.6.4a) para obtener

dhTI vzI iI d 2 hTI iI

  
1 ∂ ∂ TI
ρI CpI = kI r + kI (7.6.6)
dz r ∂r ∂r dz2

Aquí se supusieron constantes las propiedades del fluido y se tomó en cuenta el hecho de que los
límites de integración no dependen de z. El primer término del lado derecho de la ecuación anterior
puede desarrollarse como sigue:
   r=r
Z i   r=r
Z i  
1 ∂ ∂ TI 2kI 1 ∂ ∂ TI 2kI ∂ TI 2 ∂ TI
kI r = 2 r rdr = 2 ∂ r = 2 kI ri (7.6.7)
r ∂r ∂r ri r ∂r ∂r ri ∂r ri ∂r r=ri
r=0 r=0

Sustituyendo la condición de frontera entre las dos corrientes de fluidos en la ecuación (7.6.6)
resulta, tras multiplicar ambos lados por πri2

dhTI vzI iI d 2 hTI iI

πri2 ρI CpI = −2πreUe (hTI iI − hTII iII ) + πri2 kI (7.6.8)
dz dz2

Por otro lado, el producto hTI vzI iI puede descomponerse en

hTI vzI iI = hTI iI hvzI iI + hT̃I iI hvzI iI + hTI iI hṽzI iI + hT̃I ṽzI iI (7.6.9)

donde T̃I y ṽz,I representan las desviaciones espaciales de la temperatura y la velocidad, respecti-
vamente. Para obtener este resultado, se utilizó la descomposición espacial ψ = hψiI + ψ̃, donde
ψ = TI , vzI y se tomó en cuenta que las cantidades promedio son constantes en el dominio de
integración. Los valores de dichas cantidades satisfacen las siguientes restricciones

hT̃I iI = 0; hṽzI iI = 0 (7.6.10)

Por lo que,

hTI vzI iI = hTI iI hvzI iI + hT̃I ṽzI iI (7.6.11)

De esta forma, la ecuación (7.6.8) puede expresarse como sigue:

dhTI iI d 2 hTI iI
πri2 hvzI iI ρI CpI = −2πreUe (hTI iI − hTII iII ) + πri2 kI∗ (7.6.12)
dz dz2
donde se definió el término de dispersión de calor como

d 2 hTI iI d 2 hTI iI dhT̃I ṽzI iI

kI∗ = k I − ρ C
I pI (7.6.13)
dz2 dz2 dz
Para simplificar más la ecuación (7.6.12) se hace la suposición que la transferencia de calor por
convección es mucho mayor que la dispersión de calor

d 2 hTI iI dhTI iI
kI∗  hvzI iI
ρ C
I pI (7.6.14)
dz2 dz
Bajo estas condiciones, se puede reducir la ecuación (7.6.12) a

dhTI iI
ṁI CpI = −2πreUe (hTI iI − hTII iII ) (7.6.15)
dz
7.6 Transferencia de calor en intercambiadores de calor de doble tubo 131

donde se usó la definición del flujo másico

ṁI = πri2 hvzI iI ρI (7.6.16)

Siguiendo pasos similares, se obtiene, para la región-II:

dhTII iII
±ṁII CpII = 2πreUe (hTI iI − hTII iII ) (7.6.17)
dz
donde

ṁII = π(R2i − re2 )hvzII iII ρII (7.6.18)

En resumen, las ecuaciones promedio para cada región son:

dhTI iI 2πreUe
=− (hTI iI − hTII iII ) (7.6.19a)
dz ṁI CpI
dhTII iII 2πreUe
= (hTI iI − hTII iII ) (7.6.19b)
dz ±ṁII CpII

Restando la ecuación (7.6.19b) a la ecuación (7.6.19a) da lugar a la ecuación diferencial de la

fuerza impulsora de la transferencia de calor en el sistema, que es la diferencia de temperaturas
(hTI iI − hTII iII ):
 
d I II 1 1
(hTI i − hTII i ) = −2πreUe + (hTI iI − hTII iII ) (7.6.20)
dz ṁI CpI ±ṁII CpII

Para expresar la ecuación (7.6.20) de una forma independiente de la dirección del flujo, se recurre a
la definición de la tasa total de transferencia de calor,

Q̇ = ṁI CpI (T2I − T1I ) = ∓ṁII CpII (T2II − T1II ) (7.6.21)

donde los subíndices “2” y “1” se refieren a las condiciones de salida y entrada, respectivamente;
mientras que los superíndices II y I denotan a las regiones II y I, respectivamente. De esta forma,
por ejemplo, T2I representa la temperatura de salida en el tubo interno. A partir de la ecuación
anterior, se deduce lo siguiente

Q̇
ṁI CpI = (7.6.22a)
(T2I − T1I )
Q̇
∓ṁII CpII = (7.6.22b)
(T2II − T1II )

Sustituyendo estas expresiones en la ecuación (7.6.20) resulta

1 d 2πreUe  I
(hTII iII − hTI iI ) = − (T2 − T1I ) − (T2II − T1II )

(7.6.23)
(hTII iII I
− hTI i ) dz Q̇
Para escribir esta ecuación en una forma más sencilla, se utilizan las siguientes definiciones

∆T = hTII iII − hTI iI (7.6.24a)

∆T1 = T1II − T1I (7.6.24b)
∆T2 = T2II − T2I (7.6.24c)
132 Capítulo 7. Modelos macroscópicos

Entonces, la ecuación (7.6.23) se escribe como

1 d∆T 2πreUe
=− (∆T1 − ∆T2 ) (7.6.25)
∆T dz Q̇
Esta ecuación se integra desde z = 0, donde ∆T = ∆T1 , hasta z = L, donde ∆T = ∆T2 . El resultado
es
∆T2 Ue Ae
ln = (∆T2 − ∆T1 ) (7.6.26)
∆T1 Q̇
donde Ae = 2πre L. De esta ecuación, se puede despejar fácilmente a la tasa total de calor como
sigue

Q̇ = Ue Ae ∆Tln (7.6.27)

aquí se usó la definición de la media logarítmica de temperaturas:

∆T2 − ∆T1 ∆T1 − ∆T2
∆Tln = = (7.6.28)
ln ∆T2
∆T1 ln ∆T1
∆T2

Note de la ecuación (7.6.22) que, para el caso en que ṁI CPI = ṁII CPII y se opere bajo condiciones
de flujo en contra-corriente, se tiene que ∆T1 = ∆T2 y, por tanto, la media logarítmica de temperatura
estaría indeterminada. En este caso, es conveniente tomar el límite de la ecuación (7.6.28) cuando
∆T1 → ∆T2 y usar la regla de L’Hôpital para obtener

lı́m ∆Tln = ∆T1 = ∆T2 (7.6.29)

∆T1 →∆T2

y entonces

Q̇ = Ue Ae (T II − T I ) (7.6.30)

donde T II − T I = ∆T1 = ∆T2 .

 7.7 Transferencia de calor en un pozo geotérmico 133

7.7 Transferencia de calor en un pozo geotérmico

Un pozo geotérmico consiste en un tubo (o un arreglo de tubos) que se insertan en el subsuelo
con el fin de calentar agua para posteriormente generar electricidad. Desarrolle una ecuación
macroscópica para describir la transferencia de calor en el medio poroso.
Dentro del tubo hay dos fluidos que son el agua en fase líquida (fase-β ) y en fase gaseosa
(fase-γ), cuyas ecuaciones de conservación de energía térmica son:
∂ Tβ

ρβ c pβ + ∇ · ρβ c pβ vβ Tβ − kβ ∇Tβ = 0, en la fase-β (7.7.1a)
∂t
∂ Tγ

ργ c pγ + ∇ · ργ c pγ vγ Tγ − kγ ∇Tγ = 0, en la fase-γ (7.7.1b)
∂t
Las cuales están acopladas en la interfase, mediante las siguientes

nβ γ · (−kβ ∇Tβ + kγ ∇Tγ ) = nβ γ · ρβ (vβ − w)∆Hvap , en la interfase β − γ (7.7.1c)

Tβ = Tγ = Tsat , en la interfase β − γ (7.7.1d)

Además, están acopladas con el transporte de calor con el medio poroso de acuerdo a las siguientes
condiciones de frontera:

nβ ω · kβ ∇Tβ = hβ (Tβ ω − Tω ), en la interfase β − ω (7.7.1e)

−nβ ω · kβ ∇Tβ = −nβ ω · kω ∇Tω (7.7.1f)

−nγω · kγ ∇Tγ = hγ (Tγω − Tω ), en la interfase γ − ω (7.7.1g)

nγω · kγ ∇Tγ = −nγω · kω ∇Tω (7.7.1h)

En el medio poroso la ecuación de energía térmica es:

∂ Tω
ρω c pω = ∇ · (kω ∇Tω ) (7.7.1i)
∂t
En la superficie superior del medio poroso (Aωa ), podemos imponer la ley de enfriamiento de
Newton

−nωa · kω ∇Tω = ha (Tω − Ta ), en Aωa (7.7.1j)

En las demás fronteras lejos del pozo, podemos suponer que la temperatura es constante (T∞ ), por
tanto ya no se da más transferencia de calor y tenemos:

nω∞ · kω ∇Tω = 0, en Aω,∞ (7.7.1k)

Dado que nuestro objetivo es desarrollar una ecuación de transporte de calor macroscópica para la
temperatura en el medio poroso, definimos el operador de promediado:
1
Z
hTω i = Tω dV (7.7.2)
Vω
Vω

Aplicando este operador a la ecuación del medio poroso obtenemos:

dhTω i 1
Z
ρω c pω = nω · kω ∇Tω dA (7.7.3)
dt Vω
Aω
134 Capítulo 7. Modelos macroscópicos

Usando las condiciones de frontera interfaciales tenemos:

dhTω i 1 1 1
Z Z Z
ρω c pω = hβ (Tβ ω − Tω ) dA + hγ (Tγω − Tω ) dA + ha (Ta − Tω ) dA
dt Vω Vω Vω
Aβ ω Aγω Aωa
(7.7.4)

Suponiendo que los coeficientes interfaciales de transferencia de calor son constantes en cada
interfase y que en cada interfase es razonable aproximar Tω ≈ hTω i, entonces:

dhTω i hβ Aβ ω hγ Aγω ha Aaω

ρω c pω = (hT iβ ω − hTω i) + (hT iγω − hTω i) + (hT iaω − hTω i)
dt Vω Vω Vω
(7.7.5)

donde
1
hT iiω = Tiω dA, i = β γ, a (7.7.6)
Aiω
 Bibliografía

Aissa, A., Abdelouahab, M., Noureddine, A., Elganaoui, M., Pateyron, B. (2015). Ranz and
marshall correlations limits on heat flow between a sphere and its surrounding gas at high
temperature. Thermal Science 19, 1521–1528.

Arfken, G.B., Weber, H.J., Harris, F.E. (2012). Mathematical Methods for Physicists. Elsevier
LTD, Oxford.

Aris, R. (1990). Vectors, Tensors and the Basic Equations of Fluid Mechanics. Dover Publications
Inc.

Aris, R. (1994). Mathematical Modelling Techniques (Dover Books on Computer Science). Dover
Publications.

Auriault, J.L., Boutin, C., Geindreau, C. (2009). Homogenization of Coupled Phenomena in

Heterogenous Media. ISTE LTD.

Batchelor, G.K. (2000). An Introduction to Fluid Dynamics. Cambridge University Press.

Bear, J. (2018). Modeling Phenomena of Flow and Transport in Porous Media. Springer Internatio-
nal Publishing.

Bergman, T.L., Lavine, A.S., Incropera, F.P., DeWitt, D.P. (2017). Fundamentals of Heat and Mass
Transfer, 8th Edition. Wiley.

Besson, U. (2012). The history of the cooling law: When the search for simplicity can be an
obstacle. Science & Education 21, 1085–1110.

Bird, R.B., Stewart, W.E., Lightfoot, E.N. (2006). Transport Phenomena, Revised 2nd Edition.
John Wiley & Sons, Inc.

Bottaro, A. (2019). Flow over natural or engineered surfaces: an adjoint homogenization perspective.
Journal of Fluid Mechanics 877.
136 Capítulo 7. Modelos macroscópicos

Callen, H.B. (1985). Thermodynamics and an Introduction to Thermostatistics. John Wiley &
Sons.

Cattaneo, C. (1958). Sur une forme de i’équation de la chaleur éliminant le paradoxe d’une
propagation instantanée. Comptes Rendus de l’Acadèmie des Sciences 247, 431–433.

Fourier, J.B.J. (1878). The Analytical Theory of Heat. Cambridge University Press.

Gauss, C.F. (1813). Theoria attractionis corporum sphaeroidicorum ellipticorum homogeneorum

methodo nova tractata. Commentationes societatis regiae scientiarium Gottingensis recentiores
2, 355–378.

Gershenfeld, N. (1998). The Nature of Mathematical Modeling. Cambridge University Press.

Gray, W., Gray, G. (2017). Introduction to Environmental Modeling. Cambridge University Press.

Gray, W.G. (1975). A derivation of the equations for multi-phase transport. Chemical Engineering
Science 30, 229–233.

Gray, W.G., Miller, C.T. (2014). Introduction to the Thermodynamically Constrained Averaging
Theory for Porous Medium Systems. Springer International Publishing.

Groot, S.R.D., Mazur, P. (2011). Non-equilibrium thermodynamics. Dover Publications. dover ed

edición.

Haberman, R. (2012). Applied Partial Differential Equations with Fourier Series and Boundary
Value Problems (Featured Titles for Partial Differential Equations). Pearson. fifth edición.

Herrera, I., Pinder, G. (2012). Mathematical Modeling in Science and Engineering: An axiomatic
approach. Wiley-Blackwell.

Himmelblau, D.M., Bischoff, K.B. (1968). Process Analysis and Simulation. Deterministic Systems.
John Wiley & Sons Inc.

Katsikadelis, J.T. (2016). The Boundary Element Method for Engineers and Scientists: Theory and
Applications. ACADEMIC PR INC.

Lapidus, L., Pinder, G.F. (1999). Numerical Solution of Partial Differential Equations in Science
and Engineering. Wiley-Interscience.

Maxwell, J.C. (1867). On the dynamical theory of gases. Philosophical Transactions of the Royal
Society of London 157, 49–88.

Narasimhan, T.N. (1999). Fourier’s heat conduction equation: History, influence, and connections.
Reviews of Geophysics 37, 151–172.

Navier, M. (1822). Mémoire sur les lois du mouvement des fluides. volumen 6. l’Académie Royale
des Sciences.

Ranz, W., Marshall, W. (1952). Evaporation from drops. part i. Chemical Engineering Progress 48,
141–146.

Slattery, J.C. (1967). Flow of viscoelastic fluids through porous media. AIChE Journal 13,
1066–1071.

Slattery, J.C., Sagis, L., Oh, E.S. (2006). Interfacial Transport Phenomena. SPRINGER VERLAG
GMBH.
7.7 Transferencia de calor en un pozo geotérmico 137

Stokes, G.G. (1880). On the theories of the internal friction of fluids in motion, and of the
equilibrium and motion of elastic solids, en: Mathematical and Physical Papers vol.1. Cambridge
University Press, pp. 75–129.

Valdés-Parada, F., Alvarez-Ramírez, J., Ochoa-Tapia, J. (2007). Análisis de problemas de transporte

de masa y reacción mediante funciones de Green. Revista Mexicana de Ingeniería Química 6,
283–294.

Vernotte, P. (1958). Les paradoxes de la théorie continue de l’équation de la chaleur. Comptes

Rendus de l’Académie des Sciences 246, 3154–3155.

Whitaker, S. (1983). Fundamental Principles of Heat Transfer. Krieger Publishing Company.

Whitaker, S. (1988). Levels of simplification: The use of assumptions, restrictions and constraints
in engineering analysis. Chemical Engineering Education 22, 104–108.

Whitaker, S. (1992a). Introduction to Fluid Mechanics. Krieger Pub Co.

Whitaker, S. (1992b). The species mass jump condition at a singular surface. Chemical Engineering
Science 47, 1677–1685.

Whitaker, S. (1999). The Method of Volume Averaging. Springer Netherlands.

Whitaker, S. (2009). Newton’s laws, Euler’s laws and the speed of light. Chemical Engineering
Education 43, 96–103.

Wood, B.D. (2009). The role of scaling laws in upscaling. Advances in Water Resources 32,
723–736.

Wood, B.D., Valdés-Parada, F.J. (2013). Volume averaging: Local and nonlocal closures using a
Green’s function approach. Advances in Water Resources 51, 139–167.