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

net/publication/334451439

Transferencia de calor: Notas de clase

Book · July 2019

CITATIONS
0

1 author:

Francisco J. Valdes-Parada
Metropolitan Autonomous University
131 PUBLICATIONS 1,119 CITATIONS

SEE PROFILE

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

Porous Media View project

All content following this page was uploaded by Francisco J. Valdes-Parada on 14 July 2019.

The user has requested enhancement of the downloaded file.

Transferencia de calor:
2122092
Notas de clase

Francisco J. Valdés Parada y Stephen Whitaker

[email protected]; [email protected]
Última actualización: 17 de julio de 2017
2
 Índice

Página

1. Introducción 5
1.1. Hipótesis del continuo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2. Tres tipos de derivadas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3. Teoremas del transporte y de la divergencia . . . . . . . . . . . . . . . . . . . . . . . 9
1.4. Repaso de las ecuaciones gobernantes en mecánica de fluidos . . . . . . . . . . . . . 10
1.4.1. Conservación de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4.2. Conservación de cantidad de movimiento . . . . . . . . . . . . . . . . . . . . 12

2. Ecuación de energı́a 15
2.1. Ecuación de energı́a mecánica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2. Primera ley de la termodinámica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3. Ecuación de energı́a total . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.4. Ecuación de energı́a térmica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
2.5. Ley de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.6. Mecanismos de transferencia de calor . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.1. Conducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.2. Convección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.3. Radiación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.7. Condiciones de frontera . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.8. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

3. Transferencia de calor en estado estacionario 29

3.1. Conducción de calor en una placa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
3.1.1. Temperaturas superficiales especificadas . . . . . . . . . . . . . . . . . . . . . 29
3.1.2. Transferencia de calor considerando convección externa . . . . . . . . . . . . 31
3.1.3. Placa con generación de calor uniforme . . . . . . . . . . . . . . . . . . . . . 33
3.2. Conducción unidimensional en estado estacionario en coordenadas cilı́ndricas . . . . 35
3.3. Conducción de calor en coordenadas esféricas . . . . . . . . . . . . . . . . . . . . . . 38
3.4. Una introducción al escalamiento: La aleta de enfriamiento . . . . . . . . . . . . . . 40

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

4.1. Coordenadas rectangulares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.1.1. Transporte bidimensional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.2. Coordenadas cilı́ndricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46

3
 4.2.1. Transporte axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2.2. Transporte no axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
4.3. Problema de aplicación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.4. Aleta de enfriamiento revisada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
4.5. Conducción de calor transitoria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
4.6. Conducción y convección en estado transitorio . . . . . . . . . . . . . . . . . . . . . . 56

5. Capa lı́mite 59
5.1. Calentamiento súbito de una placa semi-infinita . . . . . . . . . . . . . . . . . . . . . 59
5.1.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
5.1.2. Solución exacta por el método de similaridad . . . . . . . . . . . . . . . . . . 60
5.1.3. Solución aproximada par el método de von Karman-Pohlhausen . . . . . . . . 63
5.2. Ecuaciones de capa lı́mite de energı́a y cantidad de movimiento para flujo laminar . 64
5.2.1. Planteamiento y simplificación . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.2.2. Solución aproximada de las ecuaciones de capa lı́mite . . . . . . . . . . . . . 67

6. Balances macroscópicos 71
6.1. Transferencia de calor en un panel fotovoltaico . . . . . . . . . . . . . . . . . . . . . 72
6.2. Tanque agitado homogéneo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6.2.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6.2.2. Promediado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
6.3. Tanque agitado no homogéneo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
6.4. Transferencia de calor en un pozo geotérmico . . . . . . . . . . . . . . . . . . . . . . 77
6.5. Transferencia de calor en un lecho fijo . . . . . . . . . . . . . . . . . . . . . . . . . . 78
6.6. Transferencia de calor en intercambiadores de calor de doble tubo . . . . . . . . . . . 80

Referencias 83

4
 1 | Introducción

Contenido
1.1. Hipótesis del continuo . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2. Tres tipos de derivadas . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3. Teoremas del transporte y de la divergencia . . . . . . . . . . . . . 9
1.4. Repaso de las ecuaciones gobernantes en mecánica de fluidos . . . 10
1.4.1. Conservación de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.4.2. Conservación de cantidad de movimiento . . . . . . . . . . . . . . . . 12

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 curso 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.
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. Normalmente,
en este curso 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 condiciones 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.
Para un mejor aprovechamiento del curso 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.

5
1.1. Hipótesis del continuo
Como sabemos, la materia está constituida por átomos y moléculas y éstos a su vez por partı́culas
subatómicas (protones, neutrones, electrones, etc.). El estudio del movimiento a nivel de escala
atómico está en el campo de la mecánica estadı́stica y normalmente queda fuera de los alcances de
los fenómenos de transporte, los cuales tienen su campo de aplicación a partir del nivel de escala que
le llamaremos del continuo. La diferencia entre estos dos niveles de escala consiste en que en el nivel
atómico los átomos y las moléculas están separados por espacios vacı́os, por lo que para describir,
por ejemplo, la temperatura de una masa de moléculas tendrı́amos que cuantificar la temperatura
de cada molécula en dicha masa. Por otra parte, en el nivel del continuo, las propiedades (como
temperatura, velocidad, concentración, etc.) son funciones continuas del espacio y el tiempo, esto
es, ya no identificamos los espacios vacı́os entre átomos y moléculas.
Entonces, en el nivel del continuo ya no hablamos de moléculas y vacı́o sino de fases, las cuales
a su vez están separadas entre sı́ mediante interfases. En este contexto, entendemos como fase
una porción de un material cuya estructura intermolecular no cambia, en una cierta extensión del
espacio que ocupa el material. Ası́, en un sistema homogéneo sólo hay una fase y en un sistema
heterogéneo hay más de una fase, separada mediante una o más interfases. Hablaremos con más
precisión acerca de las interfases más adelante.
Evidentemente, llevar a cabo estudios al nivel del continuo es menos realista que hacerlo al
nivel atómico o molecular. Esta pérdida de precisión muchas veces es aceptable en aplicaciones de
la ingenierı́a y se ve compensada además con la sencillez que implican los modelos matemáticos
a este nivel de escala. De manera general, puede establecerse que la hipótesis del continuo es
aplicable cuando la longitud caracterı́stica del instrumento de medición, `, que estemos utilizando
para estudiar (experimentalmente o teóricamente) un sistema sea mucho mayor1 que la trayectoria
libre media entre moléculas2 , `mol . Como lo indica Whitaker (1992), para que nuestro instrumento
de medición no pierda demasiada información, su longitud caracterı́stica debe ser también mucho
menor que la longitud caracterı́stica del sistema a estudiar, L, esto es:

`mol  `  L (1.1)

Una amplia variedad de sistemas de interés en ciencias básicas e ingenierı́a son de naturaleza
jerárquica; esto es, los procesos que tienen lugar en un cierto nivel de escala están determinados por
los procesos que se llevan a cabo en niveles de escala inferiores. En el estudio de sistemas jerárquicos,
es posible identificar un número discreto de niveles de escala que permiten definir diversos continuos.
Dichos niveles suelen ir desde el molecular, pasando por el nanoscópico (Slattery y col., 2004),
hasta llegar, por ejemplo, a la escala de campo o industrial (Whitaker, 2009). De esta forma, puede
razonarse que lo que ocurre en una interfase a nivel microscópico pueda ser determinante en el
desempeño de una planta quı́mica (ver Fig. 11 en Whitaker, 2009). En ocasiones, el pasar de un
nivel de escala a otro involucra promediar ecuaciones diferenciales parciales3 , que gobiernan el
transporte en un nivel de escala inferior (por ejemplo, del orden de micrómetros), para obtener
ecuaciones diferenciales ordinarias4 que describen el sistema en un nivel de escala superior (por
ejemplo, de metros). Con este mismo proceso es posible incluso pasar de una ecuación diferencial
1
Esto es, al menos uno o dos órdenes de magnitud mayor
2
La trayectoria libre media es la distancia promedio que puede recorrer una molécula respecto a las demás sin
incurrir en una colisión.
3
Aquellas en las que la variable dependiente es función de más de una variable independiente.
4
Aquellas en las que la variable dependiente es función de sólo una variable independiente.

6
parcial a una ecuación algebraica para describir un determinado fenómeno. De esta forma, no es
sorprendente que las ecuaciones que gobiernen el transporte de masa en una operación unitaria o en
un reactor sean considerablemente más sencillas que las ecuaciones que gobiernan el transporte en
un nivel de escala microscópico. La reducción del nivel de complejidad de los modelos matemáticos
con el incremento en la escala no es exclusiva de aplicaciones ingenieriles, por ejemplo, la ecuación
de Navier-Stokes para fluidos Newtonianos, que es válida a niveles de escala microscópicos, es
considerablemente más simple que la ecuación de Schrödinger a nivel molecular.
De acuerdo a diversos autores (Bear y Cheng, 2010; Pinder y Gray, 2008; Bird y col., 2007;
Whitaker, 1999), los modelos matemáticos en sistemas macroscópicos pueden deducirse a partir
de sus contrapartes microscópicas. Entre las herramientas necesarias para llevar a cabo estas
deducciones se encuentra el promediado y el escalamiento. Wood (2009) aclara la diferencia entre
la operación matemática de promediar (es decir, aplicar un operador de integración) y escalar. Esta
última se entiende como una reducción sistemática del número de grados de libertad involucrados en
el modelo de transporte de la escala más pequeña. Ası́, por ejemplo, si se estudia el flujo vehicular,
se puede eliminar información como el modelo o color de cada vehı́culo para describir de manera
macroscópica el aforo que tiene una cierta avenida.

Ejercicio 1.1.
En la desigualdad (1.1) se acota la longitud caracterı́stica del instrumento de medición entre dos
longitudes caracterı́sticas: `mf y L. Sin embargo, es natural pensar que existan distintas trayectorias
libres medias en las moléculas de una sustancia y que el sistema donde se encuentre tenga diversas
longitudes caracterı́sticas. Explique cómo eligirı́a de una distribución de valores de `mf y L, los que
utilizarı́a en la desigualdad (1.1). Justifique detalladamente su respuesta.

Ejercicio 1.2.
Proponga tres ejemplos donde aplique una separación de escalas caracterı́sticas que involucre tres
niveles de escala.

1.2. Tres tipos de derivadas

Un cuerpo puede definirse como un conjunto de partı́culas materiales o simplemente partı́culas,
las cuales contienen miles de moléculas (Slattery, 1999). Consideremos por el momento el
movimiento de una sola partı́cula p en el cuerpo. Denotemos el vector de posición que localiza
el centro de masa de dicha partı́cula p como rp , el cual, expresado en coordenadas cartesianas, se
escribe como:
rp (t) = x1 (t)e1 + x2 (t)e2 + x3 (t)e3 (1.2)

donde ei y xi (i = 1, 2, 3) representan los vectores unitarios en las direcciones x, y y z y los valores

correspondientes de las proyecciones sobre cada eje coordenado que localizan la posición de la
partı́cula p, respectivamente. Como es de esperarse las posiciones xi (t) (i = 1, 2, 3) son funciones
continuas del tiempo y por tanto son derivables. De hecho, la definición del vector de velocidad
de una partı́cula p, vp , se obtiene al derivar con respecto al tiempo ambos lados de la ecuación
anterior. El resultado es:
drp dx1 dx2 dx3
vp ≡ = e1 + e2 + e3 (1.3)
dt dt dt dt

7
De forma análoga, el vector de aceleración de una partı́cula p se define como:

dvp d2 x1 d2 x2 d2 x3
ap ≡ = e1 + e2 + e3 (1.4)
dt dt2 dt2 dt2
Consideremos ahora a un observador que tiene un movimiento independiente del movimiento
del flujo del fluido y que está interesado en medir una cierta cantidad A en diferentes posiciones y
tiempos. El vector de posición del observador es

ro = ro (ξ, t) (1.5)

donde ξ es el vector de posición de referencia del observador y ro = x1,o e1 + x2,o e2 + x3,o e3 . De

acuerdo a lo que aprendimos arriba, tenemos que la velocidad del observador (vo = w) es

dro
w= (1.6)
dr ξ

Las mediciones de la propiedad A que recolectó el observador son función entonces del vector de
posición del observador y del tiempo, es decir:

A = A[ro (ξ, t), t] = A[x1,o (ξ, t), x2,o (ξ, t), x3,o (ξ, t), t] (1.7)

Entonces, la razón de cambio en el tiempo de la propiedad A medida por el observador es:

dA dA dx1,o dA dx2,o dA dx3,o dA

= + + + (1.8)
dt ξ dx1,o x2,o ,x3,o ,t dt ξ dx2,o x1,o ,x3,o ,t dt ξ dx3,o x1,o ,x2,o ,t dt ξ dt ξ

donde usamos la regla de la cadena. Podemos reescribir la ecuación anterior tomando en cuenta la
definición de la velocidad del observador dada en la Ec. (1.6) como sigue:

dA dA dA dA dA
= w1 + w2 + w3 + (1.9)
dt ξ dx1,o x2,o ,x3,o ,t dx2,o x1,o ,x3,o ,t dx3,o x1,o ,x2,o ,t dt ξ

o bien
dA ∂A ∂A ∂A ∂A ∂A
= w1 + w2 + w3 + = + w · ∇A (1.10)
dt ξ ∂x1 ∂x2 ∂x3 ∂t ∂t

la cual se conoce como derivada total. De esta forma, cuando el observador se mueve con la misma
velocidad que el fluido (w = v), recuperamos la definición de la derivada material:

dA ∂A
= + v · ∇A (1.11)
dt ξ=R ∂t

Por último, cuando el observador está quieto (w = 0) tenemos la definición de la derivada parcial:

dA dA ∂A
= = (1.12)
dt ξ=O dt ξ=r ∂t

8
1.3. Teoremas del transporte y de la divergencia
En la deducción de los modelos matemáticos que utilizaremos requeriremos intercambiar
diferenciación temporal e integración espacial. Esta operación puede conseguirse mediante el
teorema general del transporte (ver sección 3.4 en Whitaker, 1992), el cual, aplicado a una propiedad
continua A puede escribirse como:
d ∂A
Z Z Z
AdV = dV + Aw · ndA (1.13)
dt ∂t
Va (t) Va (t) Aa (t)

donde Va define un dominio arbitrario que tiene forma y tamaño dependientes del tiempo y se
mueve a una velocidad w, de superficie Aa cuyo vector unitario normal a Aa es n.
Para el caso en que el dominio siempre tenga la misma cantidad de materia (es decir, es un
dominio material, Vm ) y se mueva a la velocidad del fluido v, el teorema general del transporte se
transforma en el Teorema del transporte de Reynolds:
d ∂A
Z Z Z
AdV = dV + Av · ndA (1.14)
dt ∂t
Vm (t) Vm (t) Am (t)

Por último, para el caso en que el dominio de integración está fijo en el espacio (y es por tanto
independiente del tiempo), el teorema general del transporte se reduce a:
d ∂A
Z Z
AdV = dV (1.15)
dt ∂t
V V

Con los teoremas anteriores se pueden intercambiar las operaciones de integración espacial y
diferenciación temporal. Para el caso en que el operador de integración volumétrica se aplica a una
derivada espacial, es conveniente usar el teorema de la divergencia para transformar esta integral
a una integral de superficie: Z Z
∇ · σdV = n · σdA (1.16)
Va Aa

este teorema puede aplicarse igualmente a dominios materiales o fijos e incluso a gradientes de
variables escalares.
Ejercicio 1.3.
Utilice el teorema de la divergencia aplicado a un vector a en un dominio material:
Z Z
∇ · a dV = a · n dA (1.17)
Vm (t) Am (t)

junto con la definición de la derivada material dada en la Ec. (1.11) para mostrar que el teorema
del transporte de Reynolds puede expresarse como:
Z  
d DA
Z
AdV = + A∇ · v dV (1.18)
dt Dt
Vm (t) Vm (t)

9
Ejercicio 1.4.
Demuestre que si A = 1, el teorema del transporte de Reynolds indica que

dVm
Z
= ∇ · vdV
dt
Vm (t)

R
donde Vm = dV es el volumen del dominio material.
Vm (t)

1.4. Repaso de las ecuaciones gobernantes en mecánica de fluidos

En ciencias básicas e ingenierı́a, las leyes fı́sicas que se utilizan para describir sistemas a diversos
niveles de descripción son leyes de tipo conservativas. Las cuales, en una buena parte, son leyes de
conservación de masa, energı́a y cantidad de movimiento. Las semejanzas entre estas leyes básicas
han sido remarcadas por varios autores (cf . Bird y col., 2007) y se pueden resumir, en forma tal
vez demasiado general, como sigue
( ) ( ) ( )
Tasa de Tasa de Tasa de
− +
entrada de ξ salida de ξ generación de ξ
( ) ( ) (1.19)
Tasa de Tasa de
− =
consumo de ξ acumulación de ξ
donde ξ puede representar masa, energı́a o cantidad de movimiento. La aplicación de estas leyes
de conservación a un sistema determinado lleva a ecuaciones que se suelen llamar ecuaciones de
balance5 . De esta manera, es posible desarrollar modelos matemáticos a partir de la aplicación
de un balance de materia, energı́a o cantidad de movimiento (o incluso combinaciones de ellos)
a un sistema en particular. Al estudio de los modelos resultantes de dichos balances se le conoce
como Fenómenos de transporte (Bird y col., 2007). En los siguientes párrafos se desarrollarán las
ecuaciones correspondientes a la conservación de masa, cantidad de movimiento y energı́a.

1.4.1. Conservación de masa

Si se trata con sistemas no relativistas, el principio de conservación de masa puede enunciarse
como sigue (Callen, 1985):
( )
La razón de cambio en el tiempo
=0 (1.20)
de la masa de un cuerpo

Comparando este enunciado con la ecuación de balance proporcionada anteriormente (Ec. 1.19), se
puede concluir que no hay entradas ni salidas en el sistema y que las tasas de generación y consumo
de la masa de un cuerpo se deben cancelar entre sı́. Lo anterior, es otra forma de enunciar que la
masa de un cuerpo no se crea ni se destruye, sólo se transforma.
5
La palabra balance se usa en el mismo sentido que en contadurı́a, esto es: los ingresos menos las pérdidas dan
como resultado el cambio neto del capital (Herrera y Pinder, 2012).

10
 Para expresar la Ec. (1.20) en lenguaje matemático, se empieza por definir la masa, m, de un
cuerpo como la integral de todos los diferenciales de masa que componen al cuerpo,
Z
m= dm (1.21)
Vm

donde Vm denota el dominio en el espacio ocupado por el cuerpo es un volumen material cuyas
propiedades discutimos anteriormente.
El diferencial de masa, dm está relacionado con la densidad mediante

dm = ρ dV (1.22)

Por lo que la Ec. (1.21) puede expresarse como sigue,

Z
m= ρ dV (1.23)
Vm

Además, la forma de representar matemáticamente a la razón de cambio en el tiempo de cualquier

cantidad es mediante la primera derivada con respecto del tiempo. Tomando en cuenta lo anterior,
es posible escribir en lenguaje matemático a la Ec. (1.20) como sigue,
dm d
Z
= ρ dV = 0 (1.24)
dt dt
Vm

Para poder avanzar es necesario intercambiar los operadores de diferenciación e integración; sin
embargo, como el volumen material puede ser función del tiempo, dicho intercambio no puede
hacerse de manera directa. Para ello, se recurre al teorema del transporte de Reynolds para obtener:
∂ρ
Z Z
dV + n · ρvdA = 0 (1.25)
∂t
Vm Am

Usando ahora el teorema de la divergencia en el último término y agrupando los términos resultantes
en uno solo, la ecuación anterior toma la forma,
Z  
∂ρ
+ ∇ · (ρv) dV = 0 (1.26)
∂t
Vm

En este momento, cabe preguntarse, ¿bajo qué condiciones es el resultado de una integral cero?,
algunas posibilidades son: 1) cuando el integrando es cero, 2) cuando los lı́mites de integración son
iguales, 3) cuando se integra en un punto. Dado que el volumen del cuerpo no es un punto y sus
lı́mites son iguales, se deduce que el integrando debe ser cero; por lo tanto,
∂ρ
+ ∇ · (ρv) = 0 (1.27)
∂t
La cual se conoce como ecuación de continuidad y es una forma matemática de expresar que la masa
de un cuerpo se conserva. El vector de velocidad v en la Ec. (1.27) puede determinarse a partir de
la ecuación de conservación de cantidad de movimiento como se verá más adelante. Notemos que,

11
al desarrollar la derivada espacial, la ecuación anterior puede escribirse en términos de la derivada
material como sigue:
Dρ
+ ρ∇ · v = 0 (1.28)
Dt
Vale la pena notar que, para casos en los que la densidad no cambie en la posición ni en el tiempo6 ,
la ecuación anterior se reduce a
∇·v =0 (1.29)
La cual indica que la velocidad de flujos incompresibles es solenoidal.
Ejercicio 1.5.
Utilice la ecuación de continuidad y la definición de la derivada material para deducir el teorema
modificado de Reynolds:
d DB
Z Z
ρBdV = ρ dV (1.30)
dt Dt
Vm Vm
donde B es una propiedad cualquiera que es continua en el espacio y en el tiempo.

Para concluir este apartado, es conveniente listar los pasos que se siguieron para obtener la
ecuación de continuidad:
1. Se escribió el principio de conservación en una ecuación con palabras.
2. Se tradujo la ecuación del punto anterior al lenguaje matemático.
3. Se buscó obtener una ecuación que estuviera expresada en la forma de una integral de volumen
igualada a cero.
4. Al concluir que el integrando de la ecuación del punto anterior debe ser cero, se obtuvo la
ecuación buscada.
Este procedimiento se utilizará en los siguientes apartados para deducir el resto de las ecuaciones
de conservación.

1.4.2. Conservación de cantidad de movimiento

Para deducir las ecuaciones correspondientes a este tipo de transporte, es conveniente seguir los
desarrollos de Whitaker (2009) y comenzar con la primera ley de Euler, la cual contiene a las tres
leyes de la mecánica Newtoniana y puede resumirse en la siguiente ecuación con palabras:
( ) ( )
La razón de cambio en el tiempo de la Las fuerzas que actúan
= (1.31)
cantidad de movimiento lineal de un cuerpo sobre el cuerpo

El cual es un balance de cantidad de movimiento lineal7 ; esta propiedad se define como el producto
de la masa por la velocidad, de esta forma,
( )
La razón de cambio en el tiempo de la d
Z
= ρv dV (1.32)
cantidad de movimiento lineal de un cuerpo dt
Vm
6
Como es el caso de algunos fluidos y sólidos a ciertas condiciones de temperatura y presión.
7
También conocida como momento.

12
Por otro lado, las fuerzas que actúan sobre un cuerpo pueden dividirse en fuerzas volumétricas y
superficiales. Denotando por b a la fuerza que actúa sobre todo el cuerpo por unidad de masa y
por t(n) a la fuerza superficial por unidad de área8 , se tiene entonces que
( ) Z
Las fuerzas que actúan
Z
= ρb dV + t(n) dA (1.33)
sobre el cuerpo
Vm Am

El vector de esfuerzos puede escribirse en términos del tensor de esfuerzos T, mediante el teorema
fundamental de Cauchy (la deducción de este teorema está disponible en varias fuentes, como en
el Apéndice D del trabajo de Whitaker, 2009),

t(n) = n · T (1.34)

donde n es el vector unitario normal a la superficie del cuerpo. Con estos desarrollos, se puede
expresar en lenguaje matemático la primera ley de Euler (Ec. 1.31) como sigue
d
Z Z Z
ρv dV = ρb dV + n · T dA (1.35)
dt
Vm Vm Am

Al igual que en el caso anterior, se busca agrupar a todos los términos de la ecuación anterior en
una sola integral de volumen. Para intercambiar diferenciación e integración en el término del lado
izquierdo de la ecuación anterior, se utiliza el teorema modificado de Reynolds,
Z  
d Dv ∂v
Z Z
ρv dV = ρ dV = ρ + v · ∇v dV (1.36)
dt Dt ∂t
Vm Vm Vm

Por otro lado, para convertir la integral de superficie en el último término de la Ec. (1.35) a una
integral de volumen, se recurre al teorema de la divergencia,
Z Z
n · T dA = ∇ · T dV (1.37)
Am Vm

Sustituyendo las ecs. (1.36) y (1.37) en la Ec. (1.35), resulta lo siguiente

Z    
∂v
ρ + v · ∇v − ρb − ∇ · T dV = 0 (1.38)
∂t
Vm

De donde se extrae finalmente, la ecuación de conservación de cantidad de movimiento, que se

conoce como primera ecuación de Cauchy:
 
∂v
ρ + v · ∇v = ρb + ∇ · T (1.39)
∂t
Para avanzar en nuestros desarrollos, es conveniente descomponer al tensor total de esfuerzos en
sus componentes normal y tangencial:
T = −pI + τ (1.40)
8
También conocido como vector de esfuerzos (Whitaker, 1992).

13
donde p, I y τ representan la presión termodinámica, el tensor identidad y el tensor esfuerzos
viscosos. De esta forma, para casos en los que la única fuerzas volumétrica sea la gravedad (b = g),
tenemos:
∂v
ρ + ρv · ∇v = −∇p + ρg + ∇ · τ (1.41)
∂t
la cual se conoce como ecuación de movimiento de Cauchy. Aquı́ tomamos en cuenta la siguiente
propiedad del tensor identidad:
a·I=I·a=a (1.42)
donde a es un vector cualquiera. Para flujos Newtonianos el tensor de esfuerzos viscosos se representa
como:  
h
T
i 2
τ = µ ∇v + (∇v) + I κ − µ ∇ · v (1.43)
3
donde µ y κ representan las viscosidades dinámica (o de corte) y de bulto, respectivamente.
Aplicando el operador de divergencia a la ecuación anterior tenemos:
  
h i 2
∇ · τ = µ ∇ · ∇v + ∇ · (∇v)T + ∇ κ − µ ∇ · v (1.44)
3

donde ∇ · ∇v = ∇2 v y ∇ · (∇v)T = ∇(∇ · v). Por lo tanto, para condiciones de flujo incompresible,
la ecuación de conservación de movimiento se transforma en la bien conocida ecuación de Navier-
Stokes:
∂v
ρ + ρv · ∇v = −∇p + ρg + µ∇2 v (1.45)
∂t
notemos que hemos considerado que la viscosidad es constante.

14
 2 | Ecuación de energı́a

Contenido
2.1. Ecuación de energı́a mecánica . . . . . . . . . . . . . . . . . . . . . . 15
2.2. Primera ley de la termodinámica . . . . . . . . . . . . . . . . . . . . 17
2.3. Ecuación de energı́a total . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.4. Ecuación de energı́a térmica . . . . . . . . . . . . . . . . . . . . . . . 20
2.5. Ley de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.6. Mecanismos de transferencia de calor . . . . . . . . . . . . . . . . . 23
2.6.1. Conducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.2. Convección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.6.3. Radiación . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.7. Condiciones de frontera . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.8. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27

2.1. Ecuación de energı́a mecánica

Para los desarrollos que siguen es necesario deducir la ecuación de energı́a mecánica, la cual
es el resultado de hacer el producto interno de la ecuación de movimiento con la velocidad del
fluido. Comencemos con la ecuación de movimiento de Cauchy, al hacer el producto mencionado
anteriormente, tenemos
 
∂v
ρ · v + (v · ∇v) · v = − (∇p) · v + ρg · v + (∇ · τ ) · v (2.1)
∂t

Para avanzar en nuestros desarrollos, tomemos en cuenta las siguientes identidades, las cuales
aplican a flujos incompresibles:

∂v 1 ∂v 2
·v= (2.2a)
∂t 2 ∂t
1
∇v · v = ∇v 2 (2.2b)
2
∇p · v = ∇ · (pv) (2.2c)

15
De esta forma, la Ec. (2.1) puede expresarse como:
ρ ∂v 2
 
2
+ v · ∇v = −∇ · (pv) + ρg · v + (∇ · τ ) · v (2.3)
2 ∂t
Más aún, tomando en cuenta la definición de la derivada material, podemos expresar a la Ec. (2.3)
como
ρ Dv 2
= −∇ · (pv) + ρg · v + (∇ · τ ) · v (2.4)
2 Dt
Dirigiendo la atención al término viscoso, podemos usar la identidad:
(∇ · τ ) · v = ∇ · (τ · v) − τ : ∇v (2.5)
Por lo que la Ec. (2.4) es ahora
ρ Dv 2
= −∇ · (pv) + ρg · v + ∇ · (τ · v) − τ : ∇v (2.6)
2 Dt
Ejercicio 2.1.
Demuestre que la ecuación anterior puede expresarse como sigue:
ρ Dv 2
= ρg · v + ∇ · (T · v) − T : ∇v (2.7)
2 Dt

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

τ = µ ∇v + (∇v)T
 
(2.8)
Por lo tanto,
∇ · (τ · v) = µ∇ · (∇v · v) + µ∇ · (∇v)T · v


(2.9a)
O bien, tomando en cuenta que ∇v · v = 21 ∇v 2 , tenemos
µ 2 2
∇ v + µ∇ · (∇v)T · v


∇ · (τ · v) = (2.10)
2
El resultado anterior puede modificarse al notar que
∇ · (∇v)T · v = ∇ · (v · ∇v) = ∇v : ∇v + v · ∇ (∇ · v) = ∇v : ∇v


(2.11)
En la última igualdad se tomó en cuenta la suposición de flujo incompresible. De esta forma, la Ec.
(2.6) toma la siguiente forma:
ρ Dv 2 1
= −∇ · (pv) + ρg · v + µ∇2 v 2 + µ∇v : ∇v − τ : ∇v (2.12)
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 (2.13)
Por lo que la forma final de la ecuación de energı́a mecánica es:
ρ ∂v 2
 
1
+ v · ∇v 2 = −∇ · (pv) + ρg · v + µ∇2 v 2 − µ∇vT : ∇v (2.14)
2 ∂t 2

16
Ejercicio 2.2.
Deduzca el resultado anterior a partir del producto punto de la ecuación de Navier-Stokes con el
vector de velocidad.

2.2. Primera ley de la termodinámica

Ahora que hemos repasado la deducción de las ecuaciones de conservación de masa, de cantidad
de movimiento y de energı́a mecánica, podemos dirigir la atención hacia la conservación de energı́a.
El principio fundamental de conservación de energı́a es la primera ley de la termodinámica y puede
enunciarse como sigue (Crawford, 1963)
( ) ( )
El cambio de energı́a interna Suma del calor y trabajo que actúan
= (2.15)
cinética y potencial de un sistema sobre la superficie del sistema

Esto es:
∆U + ∆Ec + ∆Ep = Q + W (2.16)

Donde
t=t
Z 2
∆ψ = ψ(t2 ) − ψ(t1 ) = dψ; ψ = U, Ec , Ep (2.17a)
t=t1

t=t
Z 2
Q= Q̇dt (2.17b)
t=t1
t=t
Z 2
W = Ẇ dt (2.17c)
t=t1

Por lo tanto, la primera ley de la termodinámica también puede escribirse como una sola integral
temporal:
t=t
Z 2h   i
dU + dEc + dEp − Q̇ + Ẇ dt = 0 (2.18)
t=t1

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

dU dEc dEp
+ + = Q̇ + Ẇ (2.19)
dt dt dt

En adelante nos referiremos a esta expresión como 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 nuestro contexto se refiere a un cuerpo o volumen

17
material. De esta forma, tenemos que:
Z
U= ρe dV (2.20a)
Vm
v2
Z
Ec = ρ dV (2.20b)
2
Vm
Z
Ep = ρφ dV (2.20c)
Vm

Donde e y φ denotan la energı́a interna por unidad de masa y la energı́a potencial por unidad de
masa, respectivamente.
Por otro lado, 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.21)
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 es: Z Z
Ẇ = tn · vdA = n · T · vdA (2.22)
Am Am

Usando estas definiciones, la forma instantánea de la primera ley de la termodinámica es ahora:

Z  
d 1
Z
ρe + ρv 2 + ρφ dV = n · (−q − qR + T · v) dA (2.23)
dt 2
Vm Am

Dirigiendo la atención al último término del el lado izquierdo del lado derecho de la ecuación
anterior tenemos, gracias al teorema modificado de Reynolds, que:
Z  
d Dφ ∂φ
Z Z
ρφ dV = ρ dV = ρ + v · ∇φ dV (2.24)
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 = − ρg · v dV (2.25)
dt
Vm Vm Vm

En la última igualdad hicimos uso de la definición del vector de gravedad:

g = −∇φ (2.26)

18
De esta forma, la Ec. (2.23) puede escribirse como:
Z  
d 1
Z Z
ρe + ρv 2 dV = n · (−q − qR + T · v) dA + ρg · v dV (2.27)
dt 2
Vm Am Vm

La cual es consistente con la ecuación con palabras que puede encontrarse en textos de fenómenos
de transporte (Bird y col., 2007; 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.28)
Notemos que en la expresión anterior denotamos 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 la sección 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 + ρg · v dV (2.29)
al cuerpo por trabajo
Am Vm
| {z } | {z }
trabajo superficial trabajo volumétrico

2.3. Ecuación de energı́a total

Para deducir la ecuación local e instantánea de conservación de energı́a es conveniente expresar la
Ec. (2.27) como una integral de volumen igualada a cero. Para ello, se utiliza el teorema modificado
de Reynolds en el término del lado izquierdo, 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.30a)
Por otro lado, utilizando el teorema de la divergencia en el primer término en el lado derecho
de la Ec. (2.27), da como resultado
Z Z
n · (−q − qR + T · v) dA = ∇ · (−q − qR + T · v) dV (2.30b)
Am Vm

El resultado de sustituir las ecs. (2.30) en la Ec. (2.27) es la siguiente expresión,

Z    
D 1 2
ρ e + v + ∇ · (q + qR − T · v) − ρg · v dV = 0 (2.31)
Dt 2
Vm

19
De donde se extrae la ecuación de energı́a total (Bird y col., 2007)
 
D 1 2
ρ e + v = −∇ · (q + qR ) + ∇ · (T · v) + ρg · v (2.32)
Dt 2
A pesar de la estructura que posee, el significado fı́sico de la ecuación permanece intacto, pues 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, mientras que el primer término del lado derecho es la tasa a la cual se suministra
calor y los últimos dos términos se refieren al trabajo realizado sobre el sistema.

2.4. Ecuación de energı́a térmica

En la sección anterior dedujimos la forma local e instantánea de la primera ley de la
termodinámica, la cual denotamos como la ecuación de energı́a total. Dicha ecuación combina tanto
las contribuciones térmicas como mecánicas de transporte de energı́a. Anteriormente dedujimos la
ecuación de energı́a mecánica, la cual, en términos del tensor total de esfuerzos es:

ρ Dv 2
= ρg · v + ∇ · (T · v) − T : ∇v (2.33)
2 Dt
Para aislar las contribuciones exclusivas por transporte de energı́a térmica, podemos restar la Ec.
(2.33) a la Ec. (2.32) para obtener la siguiente expresión
De
ρ = −∇ · (q + qR ) + T : ∇v (2.34)
Dt
la cual se conoce como la ecuación de energı́a térmica.
Al igual que en la sección anterior, para el caso de flujo de Stokes, se puede descomponer al
tensor de esfuerzos como T = −pI + τ para expresar a la Ec. (2.34) como
De
ρ = −∇ · q − p∇ · v − τ : ∇v − ∇ · qR (2.35)
Dt
Con el fin de expresar al lado izquierdo de la Ec. (2.35) 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.36)
ρ
Aplicando el operador de derivada material en ambos lados de la ecuación anterior da como
resultado
Dh De 1 Dp p Dρ
= + − 2 (2.37)
Dt Dt ρ Dt ρ Dt
Despejando de la ecuación anterior el término De
Dt y, tomando en cuenta la forma de la ecuación de
continuidad dada por
Dρ
= −ρ∇ · v (2.38)
Dt
se obtiene el siguiente resultado
De Dh 1 Dp p
= − − ∇·v (2.39)
Dt Dt ρ Dt ρ

20
Sustituyendo la expresión anterior en la Ec. (2.35), da lugar a

Dh Dp
ρ = −∇ · q − τ : ∇v + − ∇ · qR (2.40)
Dt Dt
Por último, de cualquier texto de termodinámica (ver, por ejemplo Callen, 1985) se sabe que la
entalpı́a es función de la temperatura y la presión, y dicha dependencia puede expresarse como
sigue
Dh DT 1 Dp
= Cp + (1 − T β) (2.41)
Dt Dt ρ Dt
donde Cp y β son la capacidad calorı́fica a presión constante y el coeficiente de expansión térmica,
respectivamente. Sustituyendo la Ec. (2.41) en la Ec. (2.40) permite llegar a la siguiente expresión
 
∂T
ρCp + v · ∇T = −∇ · q + Φ (2.42)
∂t

Para hacer más breve esta expresión se utilizó la siguiente definición

Dp
Φ = −τ : ∇v + T β − ∇ · qR (2.43)
Dt
Como en el caso de la conservación de masa para una especie quı́mica, note que si Φ > 0, existe
generación de calor en el sistema y si Φ < 0 el sistema tiene pérdidas de calor.

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, consideremos
el sistema mostrado en la Figura 2.1a) 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 100 a 1500 W/m2 . Se calculó la temperatura en la superficie de la lámina superior para cada
valor de Q̇/A manteniendo constate la separación de las placas en 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.44)
A L
donde k es la 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
fue el primero en notar esta correlación en el siglo XIX. Para este caso, la ley de Fourier puede
expresarse como:
∂T
qz = −k (2.45)
∂z

21
consider this problem in subsequent paragraphs; however, we will avoid a discussion of the
radiant flux, q R , until Chapter ??.

1.4.3.3. Fourier’s law.

We begin our simplification of Eq. (1.4-43) by considering the heat flux q and the
system shown in Fig. 1.4-2, which consists of two parallel plates separated by a distance L.
The lower plate is maintained at some temperature T0 by a cool stream of fluid, while the
upper plate is maintained at a higher temperature T1 by an electric heater. The heater
output is measured in order to determine the energy per unit time Q& which flows from the
top plate, through the test material, to the cold bottom plate.

a)of the thermal conductivity

Figure 1.4-2 Experimental determination
.
If we vary the heat output Q , and40the cell spacing, L, while measuring the temperature
difference, ΔT = T − T we obtain a series of values for the heat flux Q&/ A , as a function
1 o

of L and ΔT. Here A is the area

Q̇ over30which the heat flow takes place. By experiment it has
been found that Q&/ A is a linear
A function of ΔT / L as indicated in Fig. 1.4-3. The data
20

10
400 600 800 1,000 1,200 1,400
∆T /L
b)
26
Figura 2.1: a) Esquema del dispositivo experimental para determinr la conductividad térmica; b)
Resultados experimentales del flux de calor vs ∆T /L.

22
El signo negativo se agrega para que k siempre sea positivo ya ∂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.
En forma vectorial, la ley de Fourier puede escribirse como:

q = −k∇T (2.46)

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.47)

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

correlaciones.
Utilizando la forma de la ley de Fourier dada en la Ec. (2.46), la ecuación de energı́a térmica
se expresa como:  
∂T
ρCp + v · ∇T = ∇ · (k∇T ) + Φ (2.48)
∂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. A continuación describiremos brevemente cada uno
de ellos.

2.6.1. Conducción
Este modo de transferencia de calor se basa en el contacto directo entre las partı́culas. Por
lo tanto, la conducción de calor es imposible en el vacı́o. En el experimento que describimos
anteriormente se puede despreciar la radiación debido a que los cambios de temperatura no son
demasiado altos. Cuando este no es el caso, debe considerarse la radiación. Siempre que se obedezca
la ley de Fourier podemos afirmar que se da este modo de transferencia de calor.

2.6.2. Convección
Este modo de transporte de energı́a es el resultado del movimiento macroscópico. En los cursos
de termodinámica se ha encontrado este tipo de transporte en términos que incluyen ṁcp ∆T ,
donde ṁ = ρvA es el flujo másico. En el ejemplo de nuestro experimento, el calor se transporta
por convección a nivel microscópico por el movimiento de las moléculas de aire, pero a nivel
macroscópico el transporte se considera que se da por conducción. Para hacer esta distinción más
clara, clasificaremos el transporte no radiante de la siguiente forma:

El transporte de energı́a que es proporcional al flux de masa, ρv (convección).

El transporte de energı́a que es proporcional al gradiente de temperatura, ∇T (conducción).

23
2.6.3. Radiación
Es el transporte de energı́a que se da en ausencia de materia. El transporte de energı́a radiante
se da a la velocidad de la luz por medio de fotones. Para fines ingenieriles, nos interesa 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 Steffan-Boltzmann:
qR · n = q R = εσT 4 (2.49)
donde ε 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. cabe recordar que q R 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 q R [=]W (2.50)
Considere ahora un cuerpo que se encuentra dentro de una cavidad al vacı́o y la temperatura del
cuerpo se incrementa de T2 a T1 mediante una resistencia eléctrica. Un balance de energı́a requiere
que:
( ) ( )
la tasa a la cual se transfiere energı́a la tasa a la cual se emite energı́a
=
desde el cuerpo hacia las paredes radiante desde el cuerpo
( )
la tasa a la cual se absorbe energı́a
− (2.51)
radiante por el cuerpo

donde:
( )
la tasa a la cual se emite energı́a
= εA1 σT14 (2.52)
radiante desde el cuerpo
( )
la tasa a la cual se absorbe energı́a
= αA1 σT24 (2.53)
radiante por el cuerpo

donde α 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. De esta forma:

Q̇12 = A1 σ εT14 − αT24

(2.54)

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

Q̇12 = εA1 σ T14 − T24

(2.55)

2.7. Condiciones de frontera

Las condiciones de frontera se pueden presentar en 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

24
 = 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
placas que se encuentran a 25o C en un extremo (x = 0) y a 80o C en el otro (x = L) (ver Fig.
En x = L, T = 80 °C
2.2). Esta situación se expresa de la siguiente forma:
Breviario sobre modelado Francisco J. Valdés Parada
2. Condiciones tipo Neumann: SeEn
especifica
x = 0, el Tvalor
= 25 de◦ la
C derivada de la variable
Septiembre, 2012
dependiente en alguna posición. Por ejemplo, si la placa
◦ del ejemplo anterior se
localiza en medio de otras dosEn x =queL,se encuentran
placas T = 80 C a 25°C en un extremo (x
encuentra
donde aisladael(es
se especifica decir,
valor da su densidaddependiente,
la variable de flujo de calor es cero)
en este caso,en
la un extremo Notemos
temperatura.
= 0) y a 80°C en el otro (x = L) (Figura 1.6-1). Esta situación se expresa de la
que el transporte de calor en las interfaces se da por conducción.
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

encuentra aislada (es decir, su densidad de flujo de calor es cero) en un extremo

Figura 1.6-1: Esquema de una placa de largo L en contacto con otras dos placas
Figura 2.2: Esquema de una placa de largo L en contacto con otras dos placas que se encuentran
que se encuentran
a temperaturas a temperaturas
de 25o C y 80 de 25°C y 80°C, respectivamente.
o C, 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 Q (ver Fig. 2.3). En este caso, las
condiciones de frontera son:
dT
En x = 0, k =0
dxcontacto con otras dos placas
Figura 1.6-1: Esquema de una placa de largo L en
dT
que se encuentran a temperaturas
En de
x=25°C
L,y 80°C,
−k respectivamente.
=Q
dx
Figura 1.6-2: Esquema 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.

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

Figura 2.3: Esquema de una placa de largo L en contacto en un extremo con una placa aislante y
en el otro conuna placa
una aislante
fuente y en elconstante
de calor otro con una fuente de calor constante.

25

42
 dT
En x = 0, −k = 0
dx

dT
En x = L, −k = Q
dx

Figura 1.6-3: 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 a 80o C y
en el otro con
unaun flujoa 80°C
placa de aire a Tel∞otro
y en . con un flujo de aire a T∞ .

3. Condiciones de frontera mixtas (o tipo Robin o de Cauchy): Se especifica el valor de

una combinación lineal de los valores de una función y su derivada. En su forma general
unidimensional se expresa como sigue:

df
En x = a, α f (a) + β =γ (2.56)
dx
42 x=a

donde α, β y γ son coeficientes conocidos que podrı́an (o no) ser funciones de la variable
dependiente f . Note que, si β = 0, la Ec. (2.56) 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 (Bird

y col., 2007), la cual se utiliza para acoplar el transporte convectivo y difusivo en una frontera.
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 constante, T∞ y en el otro está en contacto
con una placa a temperatura constante (ver Fig. 2.4). 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 situación se expresa como sigue:

dT
En x = L, n · q = −k = h ( T |x=L − T∞ ) (2.57)
dx x=L | {z }
| {z } convección
conducción

donde n es el vector normal a la superficie, h [=] W/m2 K es el coeficiente de transferencia

de calor y es el encargado de capturar la información esencial de la transferencia de calor que
ocurre cerca de la frontera. Por supuesto, la Ec. (2.57) puede expresarse en la misma forma
que la Ec. (2.56) con una simple manipulación algebraica.

1
En realidad 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.

26
2.8. Resumen
En este capı́tulo dedujimos tres ecuaciones de transporte de energı́a. La más completa es la
ecuación de energı́a total:
 
D 1 2
ρ e + v = −∇ · (q + qR ) + ∇ · (T · v) + ρg · v (2.58)
Dt 2

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

ρ Dv 2
= ρg · v + ∇ · (T · v) − T : ∇v (2.59)
2 Dt
y de energı́a térmica
De
= −∇ · (q + qR ) + T : ∇v
ρ (2.60)
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.61)
| {z∂t} | {z } | {z } |{z}
Convección Fuente o sumidero
Conducción
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.62)
| {z∂t} | {z } |{z}
Fuente o sumidero
Conducción
Acumulación

La cual en estado estacionario se reduce a:

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

En el próximo capı́tulo estudiaremos problemas de aplicación de esta forma simplificada de la

ecuación de energı́a térmica.

27
28
 3 | Transferencia de calor en estado
estacionario

Contenido
3.1. Conducción de calor en una placa . . . . . . . . . . . . . . . . . . . 29
3.1.1. Temperaturas superficiales especificadas . . . . . . . . . . . . . . . . . 29
3.1.2. Transferencia de calor considerando convección externa . . . . . . . . 31
3.1.3. Placa con generación de calor uniforme . . . . . . . . . . . . . . . . . 33
3.2. Conducción unidimensional en estado estacionario en coordenadas
cilı́ndricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
3.3. Conducción de calor en coordenadas esféricas . . . . . . . . . . . . 38
3.4. Una introducción al escalamiento: La aleta de enfriamiento . . . . 40

3.1. Conducción de calor en una placa

3.1.1. Temperaturas superficiales especificadas
Considere la transferencia de calor en estado estacionario en una placa de sección transversal
rectanguar, la cual 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 T1 < T0 . Las otras
caras se mantienen aisladas. Deduzca el perfil de temperatura dentro de la placa.
En este caso, la ecuación de energı́a térmica es:

∂2T ∂2T ∂2T

0= + + (3.1a)
∂x2 ∂y 2 ∂z 2
La cual está sujeta a las siguientes condiciones de frontera:

En x = 0, T = T0 (3.1b)
En x = L, T = T1 (3.1c)
∂T
En y = 0, H =0 (3.1d)
∂y
∂T
En z = 0, W =0 (3.1e)
∂z

29
Integrando la ecuación diferencial desde y = 0 hasta y = H y desde z = 0 hasta z = W obtenemos,
tras considerar las condiciones de frontera:
y=H
Z z=W
∂2T
Z
0= dzdy (3.2)
∂x2
y=0 z=0

En este momento, es conveniente definir la temperatura promedio en las direcciones y y z:

y=H
Z z=W
1
Z
T (x) = T (x, y, z)dzdy (3.3)
HW
y=0 z=0

Dado que los lı́mites de integración en la Ec. (3.2) no dependen de la dirección x, es posible
intercambiar diferenciación e integración libremente para obtener:
d2 T
0= (3.4a)
dx2
La cual está sujeta a las siguientes condiciones de frontera:
En x = 0, T = T0 (3.4b)
En x = L, T = T1 (3.4c)
Aquı́ supusimos que T0 y T1 son constantes en las superficies localizadas en x = 0 y x = L.
La solución de este problema de valor a la frontera es
x
T = T0 + (T1 − T0 ) (3.5)
L
La evaluación de la solución se muestra en la Figura 3.1. Donde se puede apreciar que la pendiente
de la resta aumenta a medida que aumenta la diferencia de temperaturas entre los dos extremos.
Más aún, si definimos el promedio de la temperatura como:
Z y=H
x=L Z z=W x=L
1 1
Z Z
hT i = T (x, y, z)dzdydx = T dx (3.6)
LHW L
x=0 y=0 z=0 x=0

Obtenemos que, en este caso, su valor es la media aritmética:

T1 + T0
hT i = (3.7)
2
Por último, el flux de calor en la dirección x es:
∂T
qx = −k (3.8)
∂x
O bien, al promediar en las direcciones y y z:
dT T0 − T1 fuerza impulsora
q x = −k = = (3.9)
dx L/k resistencia
Donde es evidente que, a mayor longitud (o menor conductividad térmica), habrá mayores
resistencias a la transferencia de calor.

30
 Heat and Mass Transfer Fundamentals
Fall, 2009

Figura 3.1: Evaluación del2.4.perfil

Figure de temperatura
Temperature en with
profiles in a slab unaspecified
placa con temperaturas fijas en sus
temperatures.
extremos.

Ejercicio 3.1.
Eqs. (2.3-6) and (2.3-7) can be used to evaluate the constants of integration, A and B, and
Reconsidere the
este problema
temperature para
field el caso
is given by en que en x = L en lugar de especificar la temperatura, se
especifica el valor del flux de calor, q1 .
x
T = T0 − (T0 − T1 ) (2.3-9)
L
3.1.2. Transferencia de calor considerando convección externa
Most often the engineer is interested in this result because it can be used to compute the
Considereheat
la flux or the rate
situación at which en
mostrada energy is transferred
la Figura across
3.2, en dondethe flat
un plate.
tanqueThis is givencontiene
agitado by un fluido
a temperatura constante T` , la pared del tanque tiene un espesor ` y la temperatura ambiente es
∂T ⎛ k ⎞tanque y la tasa de transferencia de calor
Ta . Determine el perfil de temperatura qx = en− lak pared ⎜ ⎟ (T0 − T1 )
= del (2.3-10)
∂x ⎝L⎠
hacia el exterior.
Siguiendoandlos desarrollos
we see that the heatdel
fluxproblema anterior, in
is constant everywhere la theecuación diferencial
slab. Shorter queit isgobierna la
is better. and
temperaturadirectly
promediada en los ejes y y z en la pared del tanque es:
proportional to the thermal conductivity and the temperature difference at the walls
and inversely proportional to the slab thickness. It is helpful to think of Eq. (2.3-10) in the
form
d2 T
0= (3.10a)
dx2
driving force
= frontera:
heat flux de
La cual está sujeta a las siguientes condiciones (2.3-11)
resistance
dT
En x = 0, n · q = k = hliq (T x=0 − T` ) (3.10b)
dx x=0
dT
17
En x = `, n · q = −k = hair (T x=` − Ta ) (3.10c)
dx x=0

La solución de este problema es:

h  i
k k` + hair
1 1
T` + hliq Ta − (T` − Ta ) x (T` − Ta )

k

T =   = T` −   x+ (3.11)
1
k hliq + k` + hair
1 1
k hliq + k` + hair
1 hliq

31
 Mechanisms of Energy Transport 15

free convection

torced convection

'L X

Figura 3.2: Esquema del perfil de temperatura alrededor de las paredes de un tanque agitado y el
ambiente.
Fig. 1.3.7b Temperature profiles.

Here hliq is the liquid -side film heat transfer coefficient, k is the thermal conductivity of the wall material,
L is the thickness of the wall, and hair is the air-side film heat transfer coefficient. The value of knowing
the overall coefficient U should be obvious, for if it is known that we can calculate the rate at which energy
is transferred from the hot liquid to the surrounding air knowing only the temperature of the hot liquid Te
and the temperature of the ambient air, Ta. Sorne representative values for film heat transfer coefficients
are given in Table 1.3-1.
32
It should be noted at this point how very important it is to l">e able to sketch temperature profiles such as
those shown in Fig. 1.3.7, and to have sorne idea how these profiles will change as the air flow rate is
changed, or the thermal conductivity of the w~ll is changed, etc. The necessity for having a good "feel" for
the temperature field would perhaps appear more reasonable if we stopped to think about sorne of the
problems that have been successfully sol ved in earlier courses. Y o u should remember from a course in
El flux de calor en la pared exterior es entonces:

dT T` − Ta fuerza impulsora
q x = −k = 1 ` 1
= (3.12)
dx x=` hliq + k + hair
resistencia

Entonces, si denotamos la superficie externa del tanque como A, la tasa de transferencia de calor
desde el tanque hacia el ambiente es:

Q̇ = Aq x = U A (T` − Ta ) (3.13)

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.14)
hliq k hair
Esta expresión se asemeja a la resistencia equivalente de circuitos eléctricos en serie.
Ejercicio 3.2.
En algunas ocasiones, el flujo de calor emitido es demasiado para ser aceptable a nivel industrial,
por lo que se agrega una capa de material aislante a la pared del tanque. Considerando que esta capa
de aislante tiene un espesor `2 y una conductividad térmica k2 < k1 , demuestre que el coeficiente
global de transferencia de calor es ahora:

1 −1
 
1 ` `2
U= + + + (3.15)
hliq k k2 hair

3.1.3. 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:
d2 T Φ
0= 2
+ (3.16a)
dx k
En x = 0, T = T0 (3.16b)
En x = L, T = T1 (3.16c)

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

x Φ̄ L2  x x
T̄ = T0 − (T0 − T1 ) + 1− (3.17)
L k 2 L L
O bien, definiendo las variables adimensionales:

T − T1 x ΦL2
Θ= ; X= ; Ω= (3.18)
T0 − T1 L 2k(T0 − T1 )
tenemos,
Θ = [1 + ΩX] (1 − X) = 1 + (Ω − 1) X − ΩX 2 (3.19)

33
 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.

a) b)

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

LaInterestingly,
temperatura all the values
promedio en laof dimensionless heat flux, Π, converge to a value of 1 at
the es:
placa
the middle of the plate. This observation Z is confirmed
x=L X=1 by taking X = 0.5 in Eq. (2.3-32).
1
Z
Note that when Ω >0, the temperature
hT i = profiles
T dx =exhibit
T dXa maximum and when Ω < (3.20) 0, a
L
minimum appears in the temperature profiles.
x=0 The
X=0position of the maximum or minimum
O temperature
bien, can be found by taking the derivative of Θ with respect to X in Eq. (2.3-29) and
finding the value of X that makes the resulting
X=1 expression to be equal to zero; the resulting
1 Ω
Z
critical value is hΘi = ΘdX = + (3.21)
2 6
X=0
1 se calcula
1
El flux de calor es ahora una función de
X crit = y−
la posición mediante: (2.3-33)
2 2Ω
dT k dΘ
q x = −k = − (T0 − T1 ) (3.22)
dx L dX
Para escribir esta expresión totalmente en función de variables adimensionales definimos:
qxL
Π= (3.23)
k(T0 − T1 )
Por lo tanto, 21
dΘ
Π=− = 1 − Ω(1 − 2X) (3.24)
dX
La evaluación de los perfiles de temperatura y del flux de calor adimensionales se muestra en
la Figura 3.3 para diferentes valores de Ω. Como puede notarse, todos los valores del flux de calor

34
convergen a 1 en la mitad de la placa. Más aún, para valores de Ω 6= 0, los perfiles incluyen un
punto máximo (o mı́nimo) 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.25)
2Ω 2 2Ω
Este resultado puede verificarse en la Figura 3.3b) al buscar los valores para los cuales Π = 0.
Sustituyendo este valor de X en la Ec. (3.19), se puede predecir el valor máximo (o mı́nimo) de la
temperatura adimensional como:
(1 + Ω)2
Θmı́n / máx = (3.26)
4Ω
Ejercicio 3.3.
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
Φ = Φ0 1 − (3.27)
L
donde Φ0 es el valor máximo de la fuente y es una constante.

3.2. Conducción unidimensional en estado estacionario en coordenadas

cilı́ndricas
Considere un tubo de dos paredes compuestas como el mostrado en la Figura 3.4, el cual en
su compartimiento interior, lleva un fluido a temperatura constante T0 , la cual es mayor que la
temperatura del ambiente, 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 2 + (3.28)
r ∂r ∂r r ∂θ ∂z 2
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.29)
∂z
y que el transporte de calor es axisimétrico:
∂T ∂T
= (3.30)
∂θ θ=0 ∂θ θ=2π

Al igual que en el caso de la placa, definimos una temperatura promedio:

θ=2π
Z z=H
1
Z
T (r) = T (r, θ, z)dθdz (3.31)
2πH
θ=0 z=0

35
 (2.4-1)

If we as sume that the process under consideration is axially symmetric, the temperat~e is independent of (}

tHere we assume that the surrounding ftuids are ftowing at a sufficiently rapid rate so that the temperature is uniform in the fluid
except in a thin region near the pipe.

(a)

To

fo r, '•
r
(b)
Fig. 2.4.1 Radial heat conduction in a composite pipe.
Figura 3.4: a) Esquema de la sección transversal de un cilindro de paredes compuestas y b) Esbozo
del perfil de temperaturas.

36
Aplicando el operador de promediado a la ecuación diferencial y tomando en cuenta las condiciones
de frontera resulta que:
 
1 d dT 1
0= r ; r0 < r < r1 (3.32a)
r dr dr
 
1 d dT 2
0= r ; r1 < r < r2 (3.32b)
r dr dr
Las condiciones de frontera son:
dT 1
En r = r0 , k1 = h0 (T 1 − T0 ) (3.32c)
dr
En r = r1 , T1 = T2 (3.32d)
dT 1 dT 2
En r = r1 , − k1 = −k2 (3.32e)
dr dr
dT 2
En r = r2 , − k2 = h1 (T 2 − Ta ) (3.32f)
dr
Para resolver este problema es conveniente definir las siguientes variables y parámetros adimensionales:
T 1 − T0 T 2 − T0 r h0 r0 h1 r0 k2
Θ1 = ; Θ2 = ; ξ= ; Bi0 = ; Bi1 = ; κ= (3.33)
Ta − T0 Ta − T0 r0 k1 k2 k1
Por tanto, la forma adimensional del problema es:
 
1 d dΘ1
0= ξ ; 1 < ξ < ξ1 (3.34a)
ξ dξ dξ
 
1 d dΘ2
0= ξ ; ξ 1 < r < ξ2 (3.34b)
ξ dξ dξ
dΘ1
En ξ = 1, = Bi0 Θ1 (3.34c)
dξ
En ξ = ξ1 , Θ1 = Θ2 (3.34d)
dΘ1 dΘ2
En ξ = ξ1 , − = −κ (3.34e)
dξ dξ
dΘ2
En ξ = ξ2 , − = Bi1 (Θ2 − 1) (3.34f)
dξ
La solución de este problema de valor a la frontera es:
κBi1 ξ2 (Bi0 ln ξ + 1)
Θ1 = h  i (3.35a)
(Bi0 ln ξ1 + 1) κBi1 ξ2 − Bi0 Bi1 ξ2 ln ξξ21 − 1
Bi0 Bi1 ξ2 ln ξξ1 + κBi1 ξ2 (Bi0 ln ξ1 + 1)
Θ2 = h  i (3.35b)
(Bi0 ln ξ1 + 1) κBi1 ξ2 − Bi0 Bi1 ξ2 ln ξξ21 − 1

Para determinar la tasa de transferencia de calor en la superficie externa, se usará la siguiente

definición:
dT 2 r2 dΘ2
Q̇ = 2πr2 L q r |r2 = −2πr2 Lk2 = 2π Lk2 (T0 − Ta ) (3.36)
r2 dr r2 r0 dξ ξ2

37
O bien,
2πξ2 Lk2 (T0 − Ta )Bi0 Bi1 2πr2 L (T0 − Ta )
Q̇ =h  i = ξ2
(3.37)
r2 (Bi0 ln ξ1 + 1) κBi1 ξ2 − Bi0 Bi1 ξ2 ln ξξ12 − 1 h0 + r2
k1 ln ξ1 + r2
k2 ln ξξ21 + 1
h1

La cual puede expresarse en su forma más compacta como sigue

Q̇ = U2 A2 (T0 − Ta ) (3.38)
r2

donde el coeficiente global de transferencia de calor basado en el área externa es:

1 −1
 
ξ2 r2 r2 ξ2
U2 = + ln ξ1 + ln + (3.39)
h0 k1 k2 ξ1 h1
Ejercicio 3.4.
Determine la tasa de transferencia de calor basada en la superficie interna del tubo, es decir, Q̇ .
r1

3.3. 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= 2 r + 2 sin θ + 2 2 (3.40a)
r ∂r ∂r r sin θ ∂θ ∂θ r sin θ ∂φ2
La cual está sujeta a las siguientes condiciones de frontera:

En r = r0 , T = T0 (3.40b)
∂T
En r = r1 , −k = h(T − Ta ) (3.40c)
∂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.40d)
∂θ θ=0 ∂θ θ=π
∂T ∂T
= (3.40e)
∂φ φ=0 ∂φ φ=2π

Para eliminar la dependencia de la temperatura con las coordenadas angular y azimutal, definimos
el siguiente operador de promediado:
Z φ=2π
θ=π
1
Z
T = T sin θdθdφ (3.41)
4π
θ=0 φ=0

38
Aplicando este operador a la ecuación diferencial considerando las condiciones de periodicidad
tenemos:  
1 d 2 dT
0= 2 r (3.42a)
r dr dr
La cual está sujeta a las condiciones de frontera:
En r = r0 , T = T0 (3.42b)
dT
En r = r1 , −k
= h(T − Ta ) (3.42c)
dr
Antes de resolver el problema, es conveniente escribirlo en términos de las siguientes variables y
parámetros adimensionales:
T − Ta r r0 hr1
Θ= ; ξ= ; ξ0 = ; Bi = (3.43)
T0 − Ta r1 r1 k
De esta forma, el problema adimensional es:
 
1 d dΘ
0= ξ2 (3.44a)
ξ 2 dξ dξ

En ξ = ξ0 , Θ=1 (3.44b)
dΘ
En ξ = 1, − = BiΘ (3.44c)
dξ
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
ξ
Θ=   (3.45)
1
1 + Bi −1
ξ0
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), tenemos:
Bi → 0, Θ=1 (3.46)
Es decir, la temperatura de la coraza se mantendrá a T0 y no disipará el calor. Por otro lado, si la
coraza no ofrece resistencias a la transferencia de calor Bi  1, tenemos:
1
−1
ξ
Bi  1, Θ= (3.47)
1
−1
ξ0
El cual es el perfil de temperatura para el caso en el que en r = r1 , T = Ta (o bien, en ξ = 1,
Θ = 0).
Por último, el flux de calor es en este caso
dT̄ k dΘ h (T0 − Ta ) 1
q̄r = −k = − (T0 − Ta ) =   2 (3.48)
dr r1 dξ 1 + Bi ξ10 − 1 ξ

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

39
 Figura 3.5: Ejemplos de aletas de enfriamiento rectangulares y cilı́ndricas.

3.4. 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.5
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
Fig. 3.6). Las dimensiones son tales que se satisfacen las siguientes desigualdades
b<LW (3.49)
bajo estas condiciones, la ecuación de energı́a térmica se reduce a
∂2T ∂2T ∂2T
+ + =0 (3.50)
∂x2 ∂y 2 ∂z 2

Figura 3.6: Esquema de una pared con aletas de enfriamiento con sus dimensiones y ubicación del
sistema coordenado para estudiar una sola aleta.

40
Para simplificar esta ecuación es conveniente llevar a cabo análisis de órdenes de magnitud de los
términos que componen a esta ecuación. Esto es:
∂2T
 
T
2
=O (3.51a)
∂x L2
∂2T
 
T
=O (3.51b)
∂y 2 W2
∂2T
 
T
2
=O 2 (3.51c)
∂z b
Dado que L > b pero son del mismo orden de magnitud, tenemos que
∂2T ∂2T ∂2T
 , (3.52)
∂y 2 ∂x2 ∂z 2
Bajo estas condiciones no es necesario considerar la dependencia de la temperatura en y, y la
ecuación diferencial se reduce a:
∂2T ∂2T
+ =0 (3.53)
∂x2 ∂z 2
La cual está sujeta a las siguientes condiciones de frontera

En x = 0,T = Tp (3.54a)
∂T
En x = L, − k = h(T − Ta ) (3.54b)
∂x
∂T
En z = −b/2, k = h(T − Ta ) (3.54c)
∂z
∂T
En z = b/2, − k = h(T − Ta ) (3.54d)
∂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 simplificar el problema, se define la temperatura promedio
z=b/2
1
Z
hT i = T dz (3.55)
b
z=−b/2

Antes de continuar, es conveniente notar que existe simetrı́a en el transporte de calor en este
problema. Es decir, el perfil de temperatura desde z = −b/2 hasta z = 0 es el mismo que de
z = 0 hasta z = b/2. Entonces, ya no es necesaria la condición de frontera en z = −b/2, y puede
reemplazarse por la siguiente condición de simetrı́a:
∂T
En z = 0, =0 (3.56)
∂z
Lo anterior motiva a modificar la definición de la temperatura promedio por:
z=b/2
2
Z
hT i = T dz (3.57)
b
z=0

41
Aplicando este operador a la Ec. (3.53), resulta
d2 hT i 2 ∂T 2 ∂T
+ − =0 (3.58)
dx2 b ∂z z=b/2 b ∂z z=0

Tomando en cuenta las condiciones de frontera en z = 0 y en z = b/2, permite expresar la ecuación

anterior como
d2 hT i 2h  
− T | z=b/2 − Ta =0 (3.59)
dx2 bk
Si la conductividad de la aleta, es suficientemente grande respecto al producto hb/2, se puede
aceptar la siguiente suposición
hT i ∼ T |z=b/2 (3.60)
De esta forma, se puede escribir la Ec. (3.59) como sigue,
d2 Θ
− m2 Θ = 0 (3.61)
dX 2
donde, se usaron las siguientes definiciones
hT i − Ta x 2hL2
Θ= ; X= ; m2 = (3.62)
Tp − Ta L bk
La Ec. (3.61) está sujeta a las siguientes condiciones de frontera

En X = 0, Θ=1 (3.63a)
dΘ
En X = 1, − = BiΘ (3.63b)
dX
En la última ecuación, se definió el número de Biot como
hL
Bi = (3.64)
k
La solución de este problema es
Bi sinh [m (1 − X)] + m cosh [m (1 − X)]
Θ= (3.65)
Bi sinh (m) + m cosh (m)
Para obtener este resultado se usaron las siguientes identidades:
sinh [m (1 − X)] = sinh (m) cosh (mX) − cosh (m) sinh (mX)
(3.66)
cosh [m (1 − X)] = cosh (m) cosh (mX) − sinh (m) sinh (mX)
Con este resultado podemos calcular la tasa de transferencia de calor que recibe la aleta como:
∂T
Q̇f = W bqx = −kW b (3.67)
∂x x=0

O bien, promediando
Bi
dhT i Wb dΘ mkW b + tanh (m)
hQ̇f i = −kW b = −k (Tp − Ta ) = (Tp − Ta ) m (3.68)
dx x=0 L dX X=0 L Bi
tanh (m) + 1
m

42
 La eficiencia de la aleta, ηf se define como la razón del calor disipado por la superficie de la
aleta entre el calor que se disiparı́a si la superficie de la aleta se mantuviera a Tp . Esta tasa de
transferencia de calor máxima es:

Q̇máx = 2LW h(Tp − Ta ) (3.69)

Entonces, la eficiencia de la aleta es

Bi tanh (m)
2
+
ηf = m m (3.70)
Bi
tanh (m) + 1
m
Aquı́ tomamos en cuenta el hecho que: mkW b/L = 2LW h/m.
Más aún, para una placa que tenga n aletas de enfriamiento de espesor b, separadas a una
distancia d − b (ver Fig. 3.7), se tiene lo siguiente:

Af = nbW ; Au = n(d − b)W ; A = ndW (3.71)

Por lo que la eficiencia global superficial, de acuerdo a la Ec. (3.70), está dada por

Bi tanh (m)
Au + Af ηf A − Af (1 − ηf ) Af b b m2 + m
η= = =1− (1 − ηf ) = 1 − + (3.72)
A A A d d Bi
tanh (m) + 1
m
Aquı́ tomamos en cuenta que Au = A − Af . Para el caso en el que Bi → 0 , el resultado anterior
se reduce a
b b tanh(m)
η =1− + (3.73)
d d m
Ejercicio 3.5.
Muestre cómo se modificarı́a el procedimiento mostrado aquı́ si no se impone la suposición de
b < L  W.

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

43
44
 4 | Transferencia de calor en más de
una dirección

Contenido
4.1. Coordenadas rectangulares . . . . . . . . . . . . . . . . . . . . . . . . 45
4.1.1. Transporte bidimensional . . . . . . . . . . . . . . . . . . . . . . . . . 45
4.2. Coordenadas cilı́ndricas . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2.1. Transporte axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.2.2. Transporte no axisimétrico . . . . . . . . . . . . . . . . . . . . . . . . 50
4.3. Problema de aplicación . . . . . . . . . . . . . . . . . . . . . . . . . . 51
4.4. Aleta de enfriamiento revisada . . . . . . . . . . . . . . . . . . . . . 53
4.5. Conducción de calor transitoria . . . . . . . . . . . . . . . . . . . . . 55
4.6. Conducción y convección en estado transitorio . . . . . . . . . . . . 56

En esta sección consideraremos situaciones más realistas que las analizadas en el capı́tulo
anterior. En especı́fico, analizaremos problemas en los que el transporte de calor se da por
conducción en más de una dirección en estados estacionario y transitorio.

4.1. Coordenadas rectangulares

4.1.1. Transporte bidimensional

Considere una viga de sección transversal rectangular de base b, alto L y profundidad W .

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= + (4.1a)
∂x2 ∂y 2

45
La cual requiere de cuatro condiciones de frontera, las cuales son:

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

Para resolver este problema es conveniente adoptar los siguientes cambios de variables
T − T0 x y b
Θ= ; X= ; Y = ; X0 = (4.2)
T1 − T0 L L L
De esta forma, el problema puede escribirse en su forma adimensional como:

∂2Θ ∂2Θ
0= + (4.3a)
∂X 2 ∂Y 2
En X = 0, X0 , Θ = 0 (4.3b)
En Y = 0, Θ = 0 (4.3c)
En Y = 1, Θ = 1 (4.3d)

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 4.3 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 ) (4.4)

Los resultados de tomar esta condición de frontera se muestran en la Figura 4.4. 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.

4.2. Coordenadas cilı́ndricas

4.2.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 + (4.5a)
r ∂r ∂r ∂z 2

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

Figura 4.2: Perfiles de temperatura adimensionales en la sección transversal de una viga.

47
Figura 4.3: Dependencia de la temperatura promedio y el promedio de la componente normal del
flux de calor en las superficies de entrada, salida y laterales del sistema con el parámetro X0

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

∂T
En r = 0, =0 (4.5b)
∂r
∂T
En r = R, − k = h(T − Ta ) (4.5c)
∂r
En z = 0, T = T0 (4.5d)
∂T
En z = L, =0 (4.5e)
∂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 = (4.6)
T0 − Ta R R R k

48
Figura 4.4: Perfiles adimensionales de temperatura y de la componente y del flux de calor a diferentes
posiciones. En todos los casos X0 =0.5

49
De esta forma, el problema de valor a la frontera es ahora:
∂2Θ
 
1 ∂ ∂Θ
0= ξ + 2 (4.7a)
ξ ∂ξ ∂ξ ∂ζ
∂Θ
En ξ = 0, =0 (4.7b)
∂ξ
∂Θ
En ξ = 1, − = BiΘ (4.7c)
∂ξ
En ζ = 0, Θ = 1 (4.7d)
∂Θ
En ζ = λ, =0 (4.7e)
∂ζ
En la Figura 4.5 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.

4.2.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 esta forma, tenemos que:
∂T
En r = R, 0 ≤ θ ≤ π, −k = h(T − Ta ) (4.8a)
∂r
∂T Q̇
En r = R, π ≤ θ ≤ 2π, −k = (4.8b)
∂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π, − = = (4.9)
∂ξ k(T0 − Ta ) A A

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

50
Figura 4.6: 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.

La dependencia de los perfiles de temperatura con la fuente de calor se muestra en la Figura 4.7
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 Ω.

4.3. 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

51
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.

52
4.4. 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
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 (4.10a)

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

En Ω1 , Tp = T0 (4.10b)

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

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

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 = (4.11)
L T0 − Ta kp kp
De esta forma, el problema adimensional es:

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

Note que κp = 1. En la Figura 4.8 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.

Ejercicio 4.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.

Ejercicio 4.2.
Estudie el mismo problema con aletas de enfriamiento en coordenadas cilı́ndricas.

53
Figura 4.7: 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 4.8: 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.

54
4.5. Conducción de calor transitoria
Considere la transferencia de calor en una placa que inicialmente se encuentraba a una
temperatura 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 (4.13a)
∂t
La cual está sujeta a las siguientes condiciones de frontera:
∂T
En x = 0, k = h(T − Ta ) (4.13b)
∂x
El resto de las superficies están aisaldas:

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

Por último, la condición inicial es:

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

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

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 (4.16a)
dt∗
Cuando t∗ = 0, hΘi = 1 (4.16b)

Cuya solución es:

∗
hΘi = e−Bit (4.17)
En la Figura 4.9 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.

55
 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 4.9: 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).

4.6. 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.
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 − β (4.18a)
2
ρβ vβ · ∇vβ = −∇pβ + µβ ∇ vβ , en la fase − β (4.18b)

Las condiciones de frontera son:

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

56
 Mientras que, para la transferencia de calor, la ecuación de energı́a térmica en cada fase es:
∂Tσ
ρσ cpσ = kσ ∇2 Tσ , en la fase-σ (4.19a)
 ∂t 
∂Tβ
ρβ cpβ + vβ · ∇Tβ = kβ ∇2 Tβ , en la fase-β (4.19b)
∂t
Las condiciones de frontera son:
1
En z = 0, Tβ = (T0 + T1 ) (4.19c)
2
∂Tβ
En z = H, =0 (4.19d)
∂z
En r = R, Tβ = Tσ (4.19e)
∂Tβ ∂Tσ
En r = R, −kβ = −kσ (4.19f)
∂r ∂r
∂Tβ
En r = R2 , −kβ =0 (4.19g)
∂r
Por último, las condiciones iniciales son:
Cuando t = 0, Tσ = T0 (4.19h)
Cuando t = 0, Tβ = T1 (4.19i)
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σβ =
vref β pin − pout T0 − T1 R R αβ kβ
ρβ vref R vref R
Re = ,Pe = (4.20)
µβ αβ
De esta forma, el problema en forma adimensional es:
∇∗ · v∗β = 0, en la fase − β (4.21a)
Rev∗β · ∇∗ v∗β = −∇∗ p∗β + ∇∗2 v∗β , en la fase − β (4.21b)
Las condiciones de frontera son:
En z ∗ = 0, p∗β = 1 (4.21c)
∗
En z = 15, p∗β = 0 (4.21d)
∗
En r = 1, 5, vz∗ = 0 (4.21e)
Note que para escribir estas ecuaciones se tomó:
(pin − pout ) R
vref = (4.22)
µβ
La forma adimensional del problema de transferencia de calor es:
∂Θσ
= ασβ ∇2 Θσ , en la fase-σ (4.23a)
∂t∗
∂Θβ
+ P ev∗β · ∇+ Θβ = ∇∗2 Θβ , en la fase-β (4.23b)
∂t∗

57
Figura 4.10: Dinámica de la temperatura promedio en la esfera para diferentes números de Péclet
tomando Re=1000.

Las condiciones de frontera son:

1
En z ∗ = 0, Θβ = (4.23c)
2
∂Θβ
En z ∗ = H, =0 (4.23d)
∂z ∗
En r∗ = 1, Θβ = Θσ (4.23e)
∂Θβ ∂Θσ
En r∗ = 1, − ∗ = −kσβ + (4.23f)
∂r ∂r
∂Θβ
En r∗ = 5, =0 (4.23g)
∂r∗
Por último, las condiciones iniciales son:

Cuando t∗ = 0, Θσ = 1 (4.23h)
∗
Cuando t = 0, Θβ = 0 (4.23i)

El resultado de la dinámica de la temperatura promedio de la esfera se muestra en la Figura

4.10.

Ejercicio 4.3.
Resuelva de nuevo el problema considerando resistencias a la transferencia de calor y considere sólo
P e = 1000 y lleva a cabo un barrido paramétrico del número de Biot.

58
 5 | Capa lı́mite

Contenido
5.1. Calentamiento súbito de una placa semi-infinita . . . . . . . . . . . 59
5.1.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . 59
5.1.2. Solución exacta por el método de similaridad . . . . . . . . . . . . . . 60
5.1.3. Solución aproximada par el método de von Karman-Pohlhausen . . . . 63
5.2. Ecuaciones de capa lı́mite de energı́a y cantidad de movimiento
para flujo laminar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
5.2.1. Planteamiento y simplificación . . . . . . . . . . . . . . . . . . . . . . 64
5.2.2. Solución aproximada de las ecuaciones de capa lı́mite . . . . . . . . . 67

5.1. Calentamiento súbito de una placa semi-infinita

5.1.1. Planteamiento del problema
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 . 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
ρcp = k∇2 T (5.1)
∂t
Para reducir la dimensión del problema definimos a la temperatura promedio como:

y=H
Z z=H
1
Z
T = T dydz (5.2)
HW
y=0 z=0

Dado que todas las caras están aisladas excepto una, la ecuación diferencial se reduce a:

∂T ∂2T
=α 2 (5.3a)
∂t ∂x

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

En x = 0, T = T1 (5.3b)
En x → ∞, T = T0 (5.3c)
Cuando t = 0, T = T0 (5.3d)

Para resolver este problema es conveniente definir a la temperatura adimensional como:

T − T0
Θ= (5.4)
T1 − T0
De esta forma, la versión adimensional del problema es:

∂Θ ∂2Θ
=α 2 (5.5a)
∂t ∂x
En x = 0, Θ = 1 (5.5b)
En x → ∞, Θ = 0 (5.5c)
Cuando t = 0, Θ = 0 (5.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 la siguiente sección examinaremos 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. Llevaremos a cabo la solución de manera exacta y posteriormente
de manera aproximada y compararemos las predicciones de ambos tipos de solución. Con estas
herramientas estaremos preparados para estudiar el problema de capa lı́mite entre un fluido y un
obstáculo.

5.1.2. Solución exacta por el método de similaridad

Una forma de transformar la ecuación diferencial parcial a una ecuación ordinaria es definiendo
una nueva variable independiente que contenga a las variables independientes del problema original.
Esto es:
Θ(x, y) = Θ[η(x, y)] (5.6)
Donde η se conoce como variable de similaridad y en este momento no conocemos su relación con
x y t. Para atender este problema, proponemos la siguiente expresión:

η = axb tc (5.7)

Usando la regla de la cadena tenemos:

∂Θ dΘ ∂η dΘ
= = acxb tc−1 (5.8)
∂t dη ∂t dη
∂Θ dΘ ∂η dΘ
= = abxb−1 tc (5.9)
∂x dη ∂x dη
2 2
 
∂ Θ ∂ b−1 c dΘ 2 2 2(b−1) 2c d Θ dΘ
2
= abx t = a b x t 2
+ ab(b − 1)xb−2 tc (5.10)
∂x ∂x dη dη dη

60
Sustituyendo estos resultados en la ecuación diferencial obtenemos:
dΘ d2 Θ dΘ
acxb tc−1 = αa2 b2 x2(b−1) t2c 2 + αab(b − 1)xb−2 tc (5.11)
dη dη dη
De donde notamos que al tomar b = 1 se elimina el último término del lado derecho y la ecuación
anterior se reduce a:
dΘ d2 Θ
acxtc−1 = αa2 t2c 2 (5.12)
dη dη
O bien
d2 Θ cx dΘ
2
− =0 (5.13)
dη αatc+1 dη
Para eliminar a x usamos la definición de η, lo que da lugar a:
d2 Θ cη dΘ
− =0 (5.14)
dη 2 αa2 t2c+1 dη
De donde es evidente que debe tomarse c = −1/2 para eliminar la dependencia con el tiempo. De
esta forma, la ecuación es ahora:
d2 Θ η dΘ
2
+ =0 (5.15)
dη 2αa2 dη
Para eliminar el factor 2α y el número 2 adicional que resultará en el denominador al integrar, se
toma a2 = 1/(4α). De esta forma, la ecuación diferencial queda como:

d2 Θ dΘ
2
+ 2η =0 (5.16a)
dη dη
Y la variable de similaridad es:
x
η=√ (5.16b)
4αt
Las condiciones de frontera son:

En η = 0, Θ=1 (5.16c)
En η → ∞, Θ=0 (5.16d)

Para resolver este problema, es conveniente hacer el cambio de variable:

dΘ
P = (5.17)
dη
De esta forma la ecuación diferencial se reduce a:
dP
+ 2ηP = 0 (5.18)
dη
cuya solución es
dΘ
P = = c1 exp(−η 2 ) (5.19)
dη

61
Integrando semidefinidamente desde ξ = 0 hasta ξ = η da como resultado:
ξ=η
Z
Θ = 1 + c1 exp(−ξ 2 )dξ (5.20)
ξ=0

Usando la condición de frontera restante permite obtener que

1 2
c1 = − ξ=∞ = −√ (5.21)
R π
exp(−ξ 2 )dξ
ξ=0

Entonces la solución del problema es:

ξ=η
2
Z
Θ=1− √ exp(−ξ 2 )dξ (5.22)
π
ξ=0

Donde el último término corresponde a la definición de la función error, por lo que:

Θ = 1 − erf(η) = erfc(η) (5.23)

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 Θ = 0.01.
Usando tablas de la función error1 puede determinarse que el valor de η que corresponde a Θ = 0.01
es:
δ
η = √ = 1.83 (5.24)
2 αt
Es decir, √
δT = 3.66 αt (5.25)
Para concluir el problema, calculemos el flux en x = 0:

∂T dΘ ∂η k(T1 − T0 ) dΘ
qx |x=0 = −k = −k(T1 − T0 ) =− √ (5.26)
∂x x=0 dη ∂x x=0 4αt dη η=0

Tomando en cuenta las ecuaciones (5.19) y (5.21), el resultado anterior toma la forma:

k
qx |x=0 = √ (T1 − T0 ) (5.27)
παt
Al comparar este resultado con la ley de enfriamiento de Newton:

qx |x=0 = h(T1 − T0 ) (5.28)

permite concluir que:

k
h= √ (5.29)
παt
1
Ver por ejemplo: https://www.geophysik.uni-muenchen.de/∼malservisi/GlobaleGeophysik2/erf tables.pdf

62
Este resultado es importante pues permite calcular los valores del coeficiente interfacial de
transferencia 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`ref `ref
Nu = =√ (5.30)
k παt
donde `ref es una longitud de referencia. A continuación explicaremos un método alternativo para
resolver de manera aproximada este problema.

5.1.3. Solución aproximada par el método de von Karman-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
Θ = a0 + a1 + a2 2 + a3 3 (5.31)
δT δT δT

Este problema se resolverá únicamente en el dominio definido por la capa lı́mite. Por tanto, la
solución debe satisfacer:

En x = 0, Θ=1 (5.32a)
En x = δ, Θ≈0 (5.32b)

Con sólo dos condiciones de frontera no es posible calcular todas las constantes del polinomio. Para
atender este problema, imponemos que, después de la capa lı́mite, no hay transferencia de calor,
por tanto:
∂Θ
En x = δ, =0 (5.32c)
∂x
Por último, dado que en x = 0 la temperatura es constante, la ecuación diferencial permite deducir
que
∂2Θ
En x = 0, =0 (5.32d)
∂x2
Considerando estas condiciones de frontera resulta que a0 = 1, a1 = −3/2, a2 = 0, a3 = 1/2 y
entonces:
3 x 1 x3
Θ=1− + (5.33)
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
∂Θ ∂Θ 3α
dx = −α = (5.34)
∂t ∂x x=0 2δT
x=0

donde
x3
 
∂Θ 3 x dδ
= − (5.35)
∂t 2 δT2 δT4 dt

63
Por tanto
x=δ
Z T
∂Θ 3 dδ
dx = (5.36)
∂t 8 dt
x=0
Entonces, la ecuación diferencial para el espesor de la capa lı́mite es:
dδT
δT = 4α (5.37)
dt
la cual está sujeta a la condición inicial:
Cuando t = 0, δT = 0 (5.38)
La solución de este problema es: √ √
δT = 2 2αt = 2.8 αt (5.39)
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 ) (5.40)
∂x x=0 5.6 αt
lo cual lleva a concluir que:
0.536k
h= √ , modelo aproximado (5.41)
αt
Este resultado es de hecho muy similar al obtenido anteriormente:
0.564k
h= √ , modelo exacto (5.42)
αt

5.2. Ecuaciones de capa lı́mite de energı́a y cantidad de movimiento

para flujo laminar
5.2.1. Planteamiento y simplificación
Considere el caso mostrado en la Figura 5.1 en donde un flujo Newtoniano, laminar,
incompresible 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,
cantidad de movimiento y energı́a térmica son:
∂vx ∂vy ∂vz
+ + =0 (5.43a)
∂x ∂y ∂z
 2
∂ vx ∂ 2 vx ∂ 2 vx
  
∂vx ∂vx ∂vx ∂p
ρ vx + vy + vz =− +µ + + (5.43b)
∂x ∂y ∂z ∂x ∂x2 ∂y 2 ∂z 2
 2
∂ 2 vy ∂ 2 vy
  
∂vy ∂vy ∂vy ∂p ∂ vy
ρ vx + vy + vz =− − ρg+ µ + + (5.43c)
∂x ∂y ∂z ∂y ∂x2 ∂y 2 ∂z 2
 2
∂ 2 vz ∂ 2 vz
  
∂vz ∂vz ∂vz ∂p ∂ vz
ρ vx + vy + vz =− +µ + + (5.43d)
∂x ∂y ∂z ∂z ∂x2 ∂y 2 ∂z 2
 2
∂2T ∂2T

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

64
Figura 5.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.

65
Antes de continuar con las condiciones de frontera, es conveniente explorar la posibilidad de reducir
el problema. Para ello, consideremos que el largo, L aunque es mayor que el ancho, h de la placa,
es mucho menor que la profundidad de la placa, W . Esto es:
h<LW (5.44)
Más aún, supongamos 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, podemos
plantear los siguientes estimados de orden de magnitud:
∂v u 
∞
=O (5.45a)
∂x L
∂v  u∞ 
=O (5.45b)
∂z W
∂2v u 
∞
= O (5.45c)
∂x2 L2
∂2v u 
∞
= O (5.45d)
∂z 2 W 2

∂2T
 
T
2
=O (5.45e)
∂x L2
∂2T
 
T
=O (5.45f)
∂z 2 W2
Por tanto, dado que L  W , podemos 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
a:
∂vx ∂vy
+ =0 (5.46a)
∂x ∂y
 2
∂ vx ∂ 2 vx
  
∂vx ∂vx ∂p
ρ vx + vy =− +µ + (5.46b)
∂x ∂y ∂x ∂x2 ∂y 2
 2
∂ 2 vy
  
∂vy ∂vy ∂p ∂ vy
ρ vx + vy =− − ρg+ µ + (5.46c)
∂x ∂y ∂y ∂x2 ∂y 2
 2
∂2T

∂T ∂T ∂ T
vx + vy =α + (5.46d)
∂x ∂y ∂x2 ∂y 2
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, tenemos que las ecuaciones
anteriores se reducen a:
∂vx ∂vy
+ =0 (5.47a)
∂x ∂y
∂ 2 vx
 
∂vx ∂vx ∂p
ρ vx + vy =− +µ 2 (5.47b)
∂x ∂y ∂x ∂y
∂ 2 vy
 
∂vy ∂vy ∂p
ρ vx + vy =− − ρg + µ 2 (5.47c)
∂x ∂y ∂y ∂y
∂T ∂T ∂2T
vx + vy =α 2 (5.47d)
∂x ∂y ∂y

66
En este momento es oportuno corregir nuestro 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 tenemos
que:  
δH
vy = O u∞ (5.48)
L
Por lo tanto, a partir de las ecuaciones de cantidad de movimiento tenemos que:
 2 
∂p u u∞
=O ρ ∞ + µ 2 (5.49a)
∂x L δ
  2 H   
∂p δH u∞ u∞ δH ∂p
=O ρ +µ 2 = O (5.49b)
∂y L L δH L ∂x
Los resultados anteriores pueden usarse en la siguiente expresión:
∂p ∂p
dp = dx + dy (5.50)
∂x ∂y
para deducir que, en base a la restricción de escala δH  L, se tiene que:
dp ∂p
≈ (5.51)
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 podemos concluir que la presión
no cambia de valor en toda la capa lı́mite (esto es, dp/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 (5.52a)
∂x ∂y
∂vx ∂vx ∂ 2 vx
vx + vy =ν (5.52b)
∂x ∂y ∂y 2
∂T ∂T ∂2T
vx + vy =α 2 (5.52c)
∂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∞ (5.53a)
En y = 0, vx = vy = 0, T = T0 (5.53b)
En y = δH , vx = u∞ , T = T∞ (5.53c)

5.2.2. Solución aproximada de las ecuaciones de capa lı́mite

Campo de velocidad
Siguiendo el método de von Karman-Pohlhausen, proponemos que el perfil de velocidad es:
y y2 y3
vx = a + b +c 2 +d 3 (5.54a)
δH δH δH

67
donde δH es una función de x. La velocidad vx está sujeta a las siguientes condiciones de frontera

En y = 0, vx = 0 (5.54b)
En y = δH , vx = u∞ (5.54c)

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, es razonable imponer que
∂vx
En y = δH , =0 (5.54d)
∂y
Por último, dado que en y = 0, vx = vy = 0, la ecuación diferencial indica que:

∂ 2 vx
En y = 0, =0 (5.54e)
∂y 2
Usando estas condiciones de frontera podemos deducir que:
vx 3 y 1 y3
= − 3 (5.55)
u∞ 2 δH 2 δH
Para determinar la componente y de la velocidad usamos la ecuación de continuidad, la cual al ser
integrada semidefinidamente desde y = 0, da como resultado:
ζ=y
∂vx
Z
vy − vy |y=0 = − dζ (5.56)
∂x
ζ=0

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∞ y 2 1 y 4 dδH
 
vy =
4 2 − 2 δ4
δH dx
(5.57)
H
Para determinar el espesor de la capa lı́mite integramos 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 = −ν =− (5.58)
∂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ν
= (5.59)
280 dx 2δH u∞
La cual, puede integrarse, sujeta a la condición de frontera

En x = 0, δH = 0 (5.60)

para obtener p
δH = 4.64 νx/u∞ (5.61)

68
Campo de temperatura
Para la temperatura procedemos de manera similar y proponemos el siguiente polinomio:

y y2 y3
T =a+b +c 2 +d 3 (5.62a)
δT δT δT

Para estimar las cuatro constantes, usamos las siguientes condiciones de frontera:

∂2T
En y = 0, T = T0 ; =0 (5.62b)
∂y 2
∂T
En y = δT , T = T∞ ; =0 (5.62c)
∂y

Usando estas condiciones de frontera, la solución queda como:

1 y3
 
3 y
T = T0 + (T∞ − T0 ) − (5.63)
2 δT 2 δT3

Esta solución, junto con las otras obtenidas anteriormente para las componentes x y y del vector
velocidad se sustituirán 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 ) (5.64)
∂x ∂y ∂y y=0 2δT
y=0

para obtener:
y=δ
Z T 
y3 y 3 dδT 3 y2 1 y4 y 2 dδH
     
3y y 1 2α
− 3 − 2 + 4 + 2 − 4 − 3 dy = − (5.65)
δH δH δT δT dx 2 δH 2 δH δT δT dx δT u∞
y=0

O bien, definiendo
R = δT /δH (5.66)
obtenemos:
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
(5.67)

Suponiendo que R es constante podemos llegar al siguiente resultado:

R3
 
R dδT α
− δT = (5.68)
10 140 dx u∞

Esta ecuación puede resolverse sujeta a la condición de frontera:

En x = 0, δT = 0 (5.69)

69
De esta forma, obtenemos:

R3
 
R αx
− δT2 = (5.70)
20 280 u∞
O bien,
R5 20αx 20α 0.93
R3 − = 2 = 2 = (5.71)
14 u∞ δH (4.64) ν Pr
Notemos 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 (5.72)

O bien, r
νx
δT = 4.64 P r−1/3 (5.73)
u∞
Con estos resultados podemos calcular el flux interfacial de calor como:
r
∂T 3k u∞ 1/3
qy |y=0 = −k = (T0 − T∞ ) = 0.323k P r (T0 − T∞ ) (5.74)
∂y y=0 2δT νx

La cual, al comprarla con la ley de enfriamiento de Newton:

qy |y=0 = hloc (T0 − T∞ ) (5.75)

da como resultado r
u∞ 1/3
hloc = 0.323k Pr (5.76)
νx
O bien, al promediar sobre toda la placa:
x=L
1 k
Z
hhloc i = hloc dx = 0.646 Re1/2 P r1/3 (5.77)
L L
x=0

donde Re = u∞ L/ν. Definiendo el número de Nusselt como:

hhloc iL
Nu = (5.78)
k
Entonces, el resultado final de este análisis es:

N u = 0.646Re1/2 P r1/3 (5.79)

El cual está en excelente concordancia con el resultado exacto que es: N u = 0.664Re1/2 P r1/3 .

70
enteras que pueden funcionar mediante este tipo de energía.

Los colectores de energía solar (Figura 2.4-3) funcionan, casi en su totalidad, bajo

condiciones transitorias en particular en los periodos de arranque durante las primeras

horas de la mañana. El objetivo es determinar la dependencia temporal de la

temperatura de la superficie de un colector solar a lo largo de un año en el Valle de

México.
6 | Balances macroscópicos
El sistema consiste en una celda solar, dado que no hay convección ni generación de

calor al interior de ella, la ecuación de energía térmica [Ec. (1.5.5)] se reduce a

Contenido ∂T
6.1. Transferencia de calor ρen = k{∇T 2
C un
p panel fotovoltaico . . . . . . . . . . . 72 (2.4.40)
6.2. 14 2 4∂3.t .
Tanque agitado homogéneo . transporte
.difusivo
. . . . . . . . . . . . . . . . . . . 73
acumulación
6.2.1. Planteamiento del problema . . . . . . . . . . . . . . . . . . . . . . . . 73
6.2.2. Promediado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
Evidentemente
6.3. la temperatura
Tanque T es
agitado no función de. la
homogéneo . . posición
. . . . . . y. el
. . tiempo,
. . . . . .por
. . lo
. que
75 para
6.4. Transferencia de calor en un pozo geotérmico . . . . . . . . . . . . 77
simplificar 6.5.
el modelo se buscará
Transferencia desarrollar
de calor unafijo
en un lecho expresión
. . . . . .en
. .el. cual
. . . .la. temperatura
. . . 78 sólo
6.6. Transferencia de calor en intercambiadores de calor de doble tubo 80
dependa del tiempo. Para ello se define la temperatura promedio en el volumen como,

1
V V∫
T = T dV (2.4.41)

Figura
Figura 2.4-3: Colector
6.1: Esquema delsolar
panelinstalado en instalado
fotovoltaico la azoteaendella edificio B de
azotea del la UAM-Iztapalapa.
edificio B de la UAM-
Iztapalapa.

71
84
6.1. Transferencia de calor en un panel fotovoltaico
Considere un solo panel fotovoltaico como el que se muestra en la Figura 6.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
ρcp = ∇ · (k∇T ) (6.1)
∂t
La cual está sujeta a las siguientes condiciones de frontera:

nnc · q = 0, en Anc (6.2)

nc · q = − nc · q r + h(T − Tamb ) , en Ac (6.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 la superficies no
expuestas se supone que no hay transferencia de calor.
Deseamos encontrar un modelo que describa la temperatura promedio en el panel como función
del tiempo. Para ello, definimos el siguiente operador de promediado volumétrico
1
Z
hT i = T dV (6.4)
V
V

donde V denota el dominio ocupado por el panel fotovoltaico (el cual no cambia con el tiempo) y
V su volumen. Aplicando este operador de promediado a la Ec. (6.1) da lugar a:
dhT i
ρcp = h∇ · (k∇T )i (6.5)
dt
Usando el teorema de la divergencia en el término conductivo obtenemos
1 1
Z Z
h∇ · (k∇T )i = nc · k∇T dA + nnc · k∇T dA (6.6)
V V
Ac Anc

la cual al usar las condiciones de frontera toma la forma siguiente:

1 1
Z Z
h∇ · (k∇T )i = nc · qr dA − h(T − Tamb ) dA (6.7)
V V
Ac Ac

Para expresar este resultado en una forma más conveniente definimos la potencia promedio incidente
al panel solar como:
1
Z
S= nc · qr dA (6.8)
Ac
Ac

Además, en la superficie del panel descomponemos la temperatura en su valor promedio y

desviaciones por lo que:
1 1 1
Z Z Z
h(T − Tamb ) dA = h(hT i − Tamb ) dA + hT̃ dA (6.9)
V V V
Ac Ac Ac

72
 Figura 6.2: Esquema de un reactor tipo tanque agitado

Suponiendo que:
1
Z
T̃ dA  hT i (6.10)
Ac
Ac
y que el coeficiente interfacial de transferencia de calor, h, es constante, tenemos:
h∇ · (k∇T )i = ac S − hac (hT i − Tamb ) (6.11)
donde ac ≡ Ac /V . Bajo estas condiciones, el modelo macroscópico es:
dhT i
ρcp = ac S − hac (hT i − Tamb ) (6.12)
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.

Tarea

1. 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

6.2. Tanque agitado homogéneo

6.2.1. Planteamiento del problema
Considere un tanque agitado como el mostrado en la Figura 6.2 en donde ingresa una corriente
de fluido a una temperatura Ten y velocidad ven constantes. Suponga que dentro del tanque tiene

73
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.
Nuestro sistema es el lı́quido contenido en el tanque agitado. Para las condiciones que especifica
el problema, la ecuación de energı́a térmica se expresa como:
∂T
ρcp + ρcp v · ∇T = ∇ · (k∇T ) + RA ∆H (6.13)
∂t
O bien, suponiendo que el flujo es incompresible:
∂T
ρcp + ∇ · (ρcp vT − k∇T ) = RA ∆H (6.14a)
∂t
Las superficies (A ) con las que está en contacto el fluido son las siguientes: superficies de entrada
Ae , de salida As , paredes Ap 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, tenemos las siguientes condiciones de frontera:

En Ae , T = Ten , v = −ven , ne , ne · k∇T = 0 (6.14b)

En As , T = Tsal , v = vsal ns , ns · k∇T = 0 (6.14c)
En Aa , T = Ta , v = −va na , na · k∇T = 0 (6.14d)
En Ap , np · k∇T = 0, v = 0 (6.14e)

6.2.2. Promediado
Comenzaremos el proceso del cambio de escala definiendo el siguiente operador de promediado:
1
Z
hT i = T dV (6.15)
V
V

Aplicando el operador de promediado a la ecuación diferencial suponiendo que la cantidad de lı́quido

en el tanque es constante, tenemos:

dhT i 1
Z
ρcp + n · (ρcp vT − k∇T ) dA = hRA ∆Hi (6.16)
dt V
A

La integral de superficie puede descomponerse en cuatro integrales para cada superficie. Usando
las condiciones de frontera tenemos
dhT i 1 1 1 1
Z Z Z
− Ten ven dA + Tsal vsal dA − Ta va dA = hRA i∆H (6.17)
dt V V V ρcp
Ae Ae Aa

Aquı́ consideramos que la entalpı́a de reacción puede suponerse como constante. Para avanzar en
nuestros desarrollos haremos las siguientes suposiciones:

1. La temperatura de entrada es constante.

74
 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 QA 1
Z Z
− ven dA + vsal dA − = hRA i∆H (6.18)
dt V V ρcp ρcp
Ae Ae

Donde definimos:
ρcp
Z
Ta va dA = QA (6.19)
V
Aa

Por otro lado, de la ecuación de continuidad tenemos que:

Z
∇ · ρv dV = 0 (6.20)
V

O bien, usando el teorema de la divergencia y las condiciones de frontera:

Z Z Z
− ρven dA + ρvsal dA − ρva dA = 0 (6.21)
Ae As Aa

Sin embargo, dado que no hay transporte de masa hacia el agitador, tenemos que ρva = 0 y la
ecuación anterior se reduce a: Z Z
ven dA = vsal dA ≡ Q̇ (6.22)
Ae As

donde Q̇ es el flujo volumétrico; además tomamos en cuenta la suposición de flujo incompresible.

Tomando este resultado en cuenta, la ecuación macroscópica de energı́a térmica toma la forma:

dhT i Q̇
ρcp = ρcp (Ten − hT i) + QA + hRA i∆H (6.23)
dt V
Recordando la definición de QA y dado que ρva = 0, el modelo se reduce a:

dhT i ρcp
ρcp = (Ten − hT i) + hRA i∆H (6.24)
dt τ

donde τ = V /Q̇ es el tiempo de residencia.

6.3. 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.

75
 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 diferencial para el fluido en
el tanque en la escala microscópica es:
∂T
ρcp + ∇ · (ρcp vT − k∇T ) = 0 (6.25a)
∂t
En este caso las superficies con las que tiene contacto el fluido son: entradas (Ae ), salidas (As ),
paredes (Ap ), Agitador (Aa ), esferas (Aσ ) y resistencia eléctrica (Ar ). Las condiciones de frontera
en cad una de estas superficies son:

En Ae , T = Ten , v = −ven , ne , ne · k∇T = 0 (6.25b)

En As , T = Tsal , v = vsal ns , ns · k∇T = 0 (6.25c)
En Aa , T = Ta , v = −va na , na · k∇T = 0 (6.25d)
En Ap , np · k∇T = 0, v = 0 (6.25e)
En Aσ , nσ · (ρcp vT − k∇T ) = −nσ · kσ ∇Tσ (6.25f)
En Aσ , nσ · kσ ∇Tσ = hσ ( Tσ |r=R − hT i) (6.25g)
En Ar , v = 0, −nr · k∇T = hr (hT i − Tr ) (6.25h)

El operador de promediado sigue siendo el mismo definido en la sección anterior, por lo que
aplicándola a la ecuación diferencial resulta que:

dhT i 1
Z
ρcp + n · (ρcp vT − k∇T ) dA = 0 (6.26)
dt V
A

Utilizando las condiciones de frontera, la ecuación anterior se transforma en:

 
dhT i ρcp  hσ Aσ hr Ar
Z Z Z
ρcp − Ten ven dA − Tsal vsal dA + Tal va dA = ( Tσ |r=R − hT i) + (Tr − hT i)
dt V V V
Ae As Aa
(6.27)

Al igual que en el caso anterior, el término de agitación desaparece pues no hay transferencia de
masa al interior del agitador. Además, suponiendo que la temperatura de entrada es constante y
que la temperatura de salida es igual a la temperatura promedio en el tanque tenemos que:

dhT i ρcp hσ Aσ hr Ar
ρcp = (Ten − hT i) + ( Tσ |r=R − hT i) + (Tr − hT i) (6.28)
dt τ V V
Como puede notarse, esta ecuación está acoplada a la temperatura en las partı́culas, las cuales
hemos supuesto que son idénticas en tamaño y forma. La ecuación de energı́a térmica en cada
partı́cula es:
∂Tσ
ρσ cpσ = ∇ · (kσ ∇Tσ ) (6.29)
∂t

76
 Tarea
1. En este ejercicio se hizo la suposició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.

2. Reconsidere el ejercicio suponiendo ahora que la temperatura en la resistencia no es

constante, sino que cambia con la posición. ¿Cómo se modifica el modelo final?

6.4. 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β
ρβ cpβ + ∇ · (ρβ cpβ vβ Tβ − kβ ∇Tβ ) = 0, en la fase-β (6.30a)
∂t
∂Tγ
ργ cpγ + ∇ · (ργ cpγ vγ Tγ − kγ ∇Tγ ) = 0, en la fase-γ (6.30b)
∂t
Las cuales están acopladas en la interfase, mediante las siguientes

nβγ · (−kβ ∇Tβ + kγ ∇Tγ ) = nβγ · ρβ (vβ − w)∆Hvap , en la interfase β − γ (6.30c)

Tβ = Tγ = Tsat , en la interfase β − γ (6.30d)

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 β − ω (6.30e)

−nβω · kβ ∇Tβ = −nβω · kω ∇Tω (6.30f)

−nγω · kγ ∇Tγ = hγ (Tγω − Tω ), en la interfase γ − ω (6.30g)

nγω · kγ ∇Tγ = −nγω · kω ∇Tω (6.30h)

En el medio poroso la ecuación de energı́a térmica es:

∂Tω
ρω cpω = ∇ · (kω ∇Tω ) (6.30i)
∂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 (6.30j)

77
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ω,∞ (6.30k)

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 (6.31)
Vω
Vω

Aplicando este operador a la ecuación del medio poroso obtenemos:

dhTω i 1
Z
ρω cpω = nω · kω ∇Tω dA (6.32)
dt Vω
Aω

Usando las condiciones de frontera interfaciales tenemos:

dhTω i 1 1 1
Z Z Z
ρω cpω = hβ (Tβω − Tω ) dA + hγ (Tγω − Tω ) dA + ha (Ta − Tω ) dA
dt Vω Vω Vω
Aβω Aγω Aωa
(6.33)

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ω

ρω cpω = (hT iβω − hTω i) + (hT iγω − hTω i) + (hT iaω − hTω i) (6.34)
dt Vω Vω Vω
donde
1
hT iiω = Tiω dA, i = βγ, a (6.35)
Aiω

6.5. Transferencia de calor en un lecho fijo

Considere ahora la transferencia de calor en un tubo de radio R y largo L empacado con
partı́culas esféricas (fase-ω) idénticas de radio R por el cual circula un fluido (fase-β). Desarrolle
un modelo macroscópico que describa la dinámica de la transferencia de calor en el fluido a lo largo
del tubo. La ecuación de energı́a térmica para el fluido es:
∂Tβ
ρβ cpβ + ρβ cpβ ∇ · (vβ Tβ ) = kβ ∇2 Tβ (6.36a)
∂t
Cuyas condiciones de frontera son:

En Ap , np · kβ ∇Tβ = 0 (6.36b)

En Ae , −nβe · kβ ∇Tβ = −nβe · kω ∇Tω (6.36c)

En Ae , −ne,β · kω ∇Tω = he ( Tω |r=R − Tβ ) (6.36d)

78
La ecuación de energı́a térmica en las partı́culas es
∂Tω
ρω cpω = ∇ · (kω ∇Tω ) (6.36e)
∂t
En este caso promediaremos en la sección transversal del tubo que está ocupada por el fluido Aβ ,
de manera que nuestro operador de promediado es:
1
Z
hTβ i = Tβ dA (6.37)
Aβ
Aβ

En este punto es conveniente distinguir la diferenciación longitudinal de la transversal en la ecuación

diferencial. haremos esta distinción mediante subı́ndices. De esta forma, la ecuación diferencial es:
∂Tβ
ρβ cpβ + ρβ cpβ ∇L · (vβ Tβ ) + ρβ cpβ ∇T · (vβ Tβ ) = kβ ∇2L Tβ + kβ ∇2T Tβ (6.38)
∂t
Aplicando el operador a la ecuación del fluido obtenemos
∂hTβ i ρβ cpβ 1
Z Z
ρβ cpβ + ρβ cpβ ∇L · hvβ Tβ i + nT · (vβ Tβ )dσ = kβ ∇2L hTβ i + nT · kβ ∇T Tβ dσ
∂t Aβ Aβ
σT σT
(6.39)
El segundo término convectivo puede descartarse pues en todas las lı́neas de contacto del fluido con
el sólido la velocidad es nula. Utilizando las condiciones de frontera interfaciales obtenemos:
∂hTβ i he `eω
ρβ cpβ + ρβ cpβ ∇L · hvβ Tβ i = kβ ∇2L hTβ i + ( Tω |r=R − hTβ i) (6.40)
∂t Aβ
Note que en el último término hemos despreciado las desviaciones de la temperatura en la interfase
con el sólido. El término convectivo puede trabajarse al descomponer a la temperatura y a la
velocidad en sus promedios y desviaciones y reconociendo que el promedio de las desviaciones es
nulo. Bajo estas condiciones, tenemos
∂hTβ i he `eω
ρβ cpβ + ρβ cpβ ∇L · (hvβ ihTβ i) = kβ ∇2L hTβ i − ρβ cpβ ∇L · hṽβ T̃β i + ( Tω |r=R − hTβ i)
∂t Aβ
(6.41)
Usualmente se define la dispersión de calor como:
k ∗ ∇L hTβ i = kβ ∇L hTβ i − ρβ cpβ hṽβ T̃β i (6.42)
Por lo que, tomando en cuenta que la dirección longitudinal es la dirección z, la ecuación de energı́a
promedio para el fluido toma la forma
∂hTβ i ∂hTβ i ∂ 2 hTβ i he `eω
ρβ cpβ + ρβ cpβ hvβ,z i = k∗ + ( Tω |r=R − hTβ i) (6.43)
∂t ∂z ∂z 2 Aβ
para casos en los que la dispersión hidrodinámica sea despreciable respecto a la convección de calor
(esto es, para valores suficientemente altos del número de Péclet), la ecuación anterior se reduce a
la forma de flujo pistón:
∂hTβ i ∂hTβ i he `eω
ρβ cpβ + ρβ cpβ hvβ,z i = ( Tω |r=R − hTβ i) (6.44)
∂t ∂z Aβ

79
6.6. Transferencia de calor en intercambiadores de calor de doble
tubo
Considere un sistema como el dibujado en la Fig. 6.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 (6.45a)
∂z r ∂r ∂r ∂z 2

∂ 2 TII
 
∂TII kII ∂ ∂TII
± ρII CpII vzII = r + kII (6.45b)
∂z r ∂r ∂r ∂z 2
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 Ec. (6.45b), el signo + corresponde a una operación co-corriente
y el signo − a la operación en contracorriente.
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 − (6.46a)
ξi
| {z }
fI (ξ)

Figura 6.3: Intercambiador de calor de doble tubo con acercamiento en la sección transversal.

80
 " #
(ξe2 − 1) ln ξ/ξe + ξe2 − ξ 2 ln ξe


II
vzII = hvzII i (6.46b)
(1 + ξe2 ) ln ξe + 1 − ξe2
| {z }
fII (ξ)

donde se usaron las siguientes definiciones

r ri re
ξ= ; ξi = ; ξe = (6.47)
Ri Ri Ri
Ası́ como los siguientes operadores de promediado
r=r
Z i
I 2
hvzI i = 2 vzI rdr (6.48a)
ri
r=0
r=R
Z i
II 2
hvzII i = 2 vzII rdr (6.48b)
Ri − re2
r=re

La condición de frontera interfacial es:

 
∂TI ∂TII hTI iI − hTII iII  
− kI ri = −kII re =   = re Ue hTI iI − hTII iII (6.49)
∂r r=ri ∂r r=re 1
+ 1
ln re
+ 1
re he k ri ri hi

Dirigiendo la atención por el momento al fluido del tubo interno (región-I), se aplica el operador
de promediado definido en la Ec. (6.48a) para obtener

dhTI vzI iI d2 hTI iI

  
1 ∂ ∂TI
ρI CpI = kI r + kI (6.50)
dz r ∂r ∂r dz 2
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 (6.51)
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 Ec. (6.50) resulta,
tras multiplicar ambos lados por πri2

dhTI vzI iI d2 hTI iI

πri2 ρI CpI = −2πre Ue (hTI iI − hTII iII ) + πri2 kI (6.52)
dz dz 2
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 (6.53)

donde T̃I y ṽz,I representan las desviaciones espaciales de la temperatura y la velocidad,

respectivamente. Para obtener este resultado, se utilizó la descomposición espacial ψ = hψiI + ψ̃,

81
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 = hṽzI iI = 0 (6.54)

Por lo que,
hTI vzI iI = hTI iI hvzI iI + hT̃I ṽzI iI (6.55)
De esta forma, la Ec. (6.52) puede expresarse como sigue:

dhTI iI d2 hTI iI
πri2 hvzI iI ρI CpI = −2πro Uo (hTI iI − hTII iII ) + πri2 kI∗ (6.56)
dz dz 2
donde se definió el término de dispersión de calor como

d2 hTI iI d2 hTI iI dhT̃I ṽzI iI

kI∗ = k I − ρI CpI (6.57)
dz 2 dz 2 dz
Para simplificar más la Ec. (6.56) se hace la suposición que la transferencia de calor por convección
es mucho mayor que la dispersión de calor

d2 hTI iI dhTI iI
kI∗  hv zI iI
ρI C pI (6.58)
dz 2 dz
Bajo estas condiciones, se puede reducir la Ec. (6.56) a

dhTI iI
ṁI CpI = −2πre Ue (hTI iI − hTII iII ) (6.59)
dz
donde se usó la definición del flujo másico

ṁI = πri2 hvzI iI ρI (6.60)

Siguiendo pasos similares, se obtiene, para la región-II:

dhTII iII
± ṁII CpII = 2πre Ue (hTI iI − hTII iII ) (6.61)
dz
donde
ṁII = π(Ri2 − re2 )hvzII iII ρII (6.62)
En resumen, las ecuaciones promedio para cada región son:

dhTI iI 2πre Ue
=− (hTI iI − hTII iII ) (6.63a)
dz ṁI CpI
dhTII iII 2πre Ue
= (hTI iI − hTII iII ) (6.63b)
dz ±ṁII CpII

Restando la Ec. (6.63b) a la Ec. (6.63a) 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 1 1
(hTI iI − hTII iII ) = −2πre Ue + (hTI iI − hTII iII ) (6.64)
dz ṁI CpI ±ṁII CpII

82
Para expresar la Ec. (6.64) 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 ) (6.65)
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 = (6.66a)
(T2I − T1I )
Q̇
∓ṁII CpII = (6.66b)
(T2II − T1II )
Sustituyendo estas expresiones en la Ec. (6.64) resulta
1 d 2πre Ue  I
(hTII iII − hTI iI ) = − (T2 − T1I ) − (T2II − T1II )

(6.67)
(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 (6.68a)
∆T1 = T1II − T1I (6.68b)
∆T2 = T2II − T2I (6.68c)
Entonces, la Ec. (6.67) se escribe como
1 d∆T 2πre Ue
=− (∆T1 − ∆T2 ) (6.69)
∆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 ) (6.70)
∆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 (6.71)
aquı́ se usó la definición de la media logarı́tmica de temperaturas:
∆T2 − ∆T1 ∆T1 − ∆T2
∆Tln = ∆T2
= (6.72)
ln ∆T1 ln ∆T1
∆T2
Note de la Ec. (6.66) que, para el caso en que ṁI CP I = ṁII CP II 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 Ec. (6.72) cuando ∆T1 →
∆T2 y usar la regla de L’Hôpital para obtener
lı́m ∆Tln = ∆T1 = ∆T2 (6.73)
∆T1 →∆T2
y entonces
Q̇ = Ue Ae (T II − T I ) (6.74)
donde T II − TI = ∆T1 = ∆T2 .

83
84
 Bibliografı́a

Arfken, G., Weber, H., Harris, F. (2012). Mathematical Methods for Physicists. Academic Press.
séptima edición.

Bear, J., Cheng, A. (2010). Modeling Groundwater Flow and Contaminant Transport. Springer.

Bird, R., Stewart, W., Lightfoot, E. (2007). Transport Phenomena. Wiley. segunda edición.

Callen, H. (1985). Thermodynamics and an Introduction to Thermostatistics. Wiley. segunda

edición.

Crawford, F. (1963). Heat, Thermodynamics, and Statistical Physics. Rupert Hart-Davis, London,
Harcourt, Brace & World, Inc.

Herrera, I., Pinder, G. (2012). Mathematical Modeling in Science and Engineering. Wiley.

Pinder, G., Gray, W. (2008). Essentials of Multiphase Flow in Porous Media. Wiley.

Slattery, J. (1999). Advanced Transport Phenomena. Cambridge University Press.

Slattery, J., Oh, E., Fu, K. (2004). Extension of continuum mechanics to the nanoscale. Chemical
Engineering Science 59, 4621–4635.

Whitaker, S. (1983). Fundamental principles of Heat Transfer. Krieger.

Whitaker, S. (1992). Introduction to Fluid Mechanics. Krieger.

Whitaker, S. (1999). The Method of Volume Averaging. Kluwer.

Whitaker, S. (2009). Chemical engineering education: making connections at interfaces. Revista

Mexicana de Ingenierı́a Quı́mica 8, 1–33.

Wood, B. (2009). The role of scaling laws in upscaling. Advances in Water Resources 32, 723–736.

85

View publication stats