1

CONVECCIN
El modo de transIerencia de calor por conveccion se compone de dos mecanismos de transporte,
que son, la transIerencia de energia debido al movimiento aleatorio de las moleculas (diIusion
termica) y el movimiento global o macroscopico del Iluido. El movimiento del Iluido del Iluidos
se asocia con el hecho de que, en cualquier instante, grandes numeros de moleculas se mueven de
Iorma colectiva o como agregados.
La transIerencia de calor se clasiIica de acuerdo con la naturaleza del Ilujo. Se habla de
conveccion for:ada cuando el Ilujo es causado por medios externos, tales como: ventilador,
bomba o vientos atmosIericos. Por otra parte, en la conveccion natural (o libre) el Ilujo es
inducido por Iuerzas de empuje que surgen a partir de la diIerencia de densidad ocasionada por la
variacion de la temperatura en los Iluidos.
La transIerencia de calor por conveccion ocurre
en una region de interIase entre un Iluido en
movimiento y una superIicie solida, entre dos
liquidos inmisibles en movimiento relativo y
entre un gas y un liquido que tiene movimiento
relativo, que estan a diIerentes temperatura. En
la Figura 4.1, se muestra el perIil de velocidad
tipico presente en la region de la interIace,
conocida como capa limite, zona en la cual
ocurre el Ienomeno de conveccion.
Sin importar la naturaleza particular del
proceso de transIerencia de calor por
conveccion, el Ilujo de calor por unidad de
area, es proporcional a la diIerencia de temperatura entre la superIicie (
s
T ) y del Iluidos que la
rodea (

T ). La ecuacion que modela adecuadamente este Ienomeno se conoce como ey de
enfriamiento de Newton y se expresa de la siguiente Iorma:
) (

= T T h q
s
(4.1)
Donde q es el Ilujo de calor por unidad de area expresada en W/m
2
y h es el coeIiciente de
conveccion (tambien llamado coeIiciente de pelicula o Iilmico) expresado en W/m
2
K.
FUNDAMENTOS DE LA CONVECCIN
PLANTEAMIENTO DE UN PROBLEMA DE CONVECCIN
onsidere la condicion de Ilujo de la Figura 4.1, donde el Iluido con velocidad

u y temperatura

T Iluye sobre una superIicie de Iorma arbitraria y de area A. Se supone que

> T T
s
por lo que
produce un Ilujo de calor local hacia el Iluido, dado por la ecuacion de Newton. La transferencia
de calor total que pasa a traves de la superIicie A se determina integrando el Ilujo de calor por
unidad de area sobre el area, es decir:
Figura 4.1. Esquema de la capa limite
hidrodinamica

u , 0

, 0 =

=
u
u y o

u


2



= =
A
s
A
dA T T h Q qdA Q ) (

(4.2)
Se deIine el coeficiente de conveccion promedio h sobre toda la superIicie A, como:
) (

= T T A h Q
s

(4.3)
Luego al iguala las ecuaciones (4.2) y (4.3) se obtiene una relacion entre el coeIiciente de
conveccion local y promedio
A hd
A
h
A

=
1
(4.4)
Para el caso de un problema sobre una placa plana de longitug , la ecuacio (4.4) se reduce a:

=

hd

h
0
1
(4.5)
omo se podra notar el problema de transIerencia de calor por conveccion consite en determinar
el coeIiciente promedio de convecion h
CAPA LIMITE DE CONVECCIN
[Link]. CAPA LIMI1E HIDRODIMAMICA
onsiderece una placa plana, como la Figura
4.2, donde un Ilujo de Iluido incide sobre
esta. Las particula de Iluido al hacer contacto
con la placa adquiren una velocidad relativa
a la placa nula. Estas particula actuan por
eIecto de la viscosidad sobre las capas
contiguas retardando su movimiento, asi
sucesivamente, hasta que a una distancia y
o los eIectos se hacen insigniIicante, esto
como conscuencia de que la velocidad del Iluido se hace igual a la velocidad de corriente libre.
La cantidad o se le llama espesor de capa limite hidrodinamica y se deIine como:

, 0 =

=
u
u y o (4.6)
uando el Iluido choca con la placa se producen dos regiones perIectamente deIinidas, la primera
region, llamada capa limite, que se caracteriza por una capa muy delgada en los que los
gradientes de velocidad y esIuerzos corte son grandes. La segunda region, llamada Ilujo de
potencia, se encuentra Iuera de la capa limite y los gradientes de velocidad y esIuerzos cortantes
son insigniIicantes.
A demas de el espesor de la capa limite, o, existen dos parametros impotente que se deber
calcular, estos son: El coeIiciente de Iriccion y el esIuerzo cortante en la pared.
El coeIiciente de Iriccion se deIine por:
Figura 4.2. Capa limite hidrodinamica


3

2
2
1

=
u
C
s
f
p
t
(4.7)
Donde p es la densidad del Iluido,

u es la velocidad de corriente libre y

=
y
s
y
u
3 t (4.8)
Donde 3 es la viscosidad dinamica del Iluido.
[Link]. CAPA LIMI1E 1ERMICA
Asi como se produce la capa limite
hidrodinamica, tambien se produce una
capa limite termica, siempre y cuando
exista un gradiente de temperatura entre
la superIicie de la placa y el Iluido que
entra en contacto con ella (ver Figura
4.3). En la entran del Iluido a la placa,
se puede considerar que la distribucion
de temperatura

T es uniIorme, luego al
estar la particula de Iluido en contacto con la placa que esta a
s
T , tal que

> T T
s
, se transmite
energia termica a las particulas de Iluido. uando se alcanza el equilibrio termico las particulas
de Iluido en contacto con la placa alcanzan la temperatura de la placa, las cuales a su vez
intercambia energia con las particulas adyacentes, produciendo un gradiente de temperatura en
esta region. Esta region donde existe el gradiente de temperatura se le llama capa limite termica y
su espesor (espesor de la capa limite termica) se determine por:

, 0 =

=
T T
T T
t
s
s
y o (4.)
Donde
s
T es la temperatura en la superIicie de la placa ( ) 0 , ( = = y T T
s
) y

T es la temperatura
del Iluido en la corriente libre
Haciendo el cambio de variable, en termino de exceso de temperatura adimensional 7

=
T T
T T
s
s
7 (4.10)
El espesor de la capa limite termica sera

, 0 =
=
7
o y
t
(4.11)
EL Ilujo de calor por unidad de area que pasa al Iluido a traves de la placa es
Figura 4.3. Capa limite termica


4

=
=

= T T h
y
T
k q
s
y
f s
0
(4.12)
Se debe destacar que las particulas de Iluido que estan en contacto con la superIicie de la placa
reciben el calor mediante conduccion, ya que no existe movimiento relativo es esta particula.
Desde estas particular de Iluido en contacto con el Iluido, el calor pasa al resto del Iluido por
conveccion.
Luego el coeIiciente de conveccion local se puede calcular por la siguiente relacion:

0 0 = =

=
y
f
y
s
f
y
k h o
y
T
T T
k
h
7
(4.13)
Ahora bien, para determinar el coeIiciente de conveccion local se debe calcular el gradiente de
temperatura en la superIicie de la pared.
Otro parametro de suma importancia en la conveccion es el numero de Nusselt, que se deIine por:

0 =

= =
y
s f

y
T
T T

k
h
Nu (4.15)
Usando

=
T T
T T
s
s
7 y haciendo

y
y =

, el numero de Nusselt se puede determinar por

=
y

y
Nu

7
(4.16)
Aqui tambien se nota que para calcular el numero de Nusselt se requiere calcular el gradiente de
temperatura adimensional en la superIicie de la pared.
ECUACIONES FUNDAMENTALES PARA CONVECCIN
omo se notara en las secciones anteriores, el estudio de la conveccion implica la combinacion
de movimiento del Ilujo, mezcla y transIerencia de calor. El movimiento del Ilujo se rige por la
ecuacion de continuidad y ecuacion de momentum, la transIerencia de calor la describe la
primera y segunda ley de la termodinamica y la mezcla la ecuacion de conservacion de las
especies.
Ecuacin de Continuidad
La ecuacion de continuidad (o ley de conservacion de la masa) esta dada en el sistema de
coordenada rectangular por:
0 ) ( 0
) ( ) ( ) (
= \ +

'
t
o
:
w
y
v

u
t
6
p
p p p p p
(4.17)
Para un Iluido incompresible, que es la situacion mas comun en la capa limite, se tiene:


5
0 0 = \ =

' o
:
w
y
v

u
6
(4.18)
En la Figura 4.4, se muestra el sistema de coordenada cilindricas y esIerica. Luego la ecuacion de
continuidad para densidad constante en coordenadas cilindrica es:
0
1
=

+ +

:
v v
r r
v
r
v
: r r
7
7
(4.1)
En coordenadas esIericas la ecuacion de continuidad para densidad constante es
0
sen
1
sen
sen
1 1
2
=

v
v v r
r r
r
(4.20)
Ecuacin de Momentum (o Balance de Fuerzas)
La ecuacion de momentum para Iluidos Newtoniano Iue propuesta por Navier-Stokes y se
expresa de la siguiente manera:

'
+

'

'
+

'

'
+

'

:
y

:
w
y
w

w
:
P
:
w
w
y
w
v

w
u
t
w

:
v
y
v

v
y
P
:
v
w
y
v
v

v
u
t
v

:
u
y
u

P
:
u
w
y
u
v

u
u
t
u
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
1
1
1
;
p
;
p
;
p
(4.21)
En Iorma vectorial, la ecuacion de Navier-Stokes es:
' P
Dt
' D
+ \ + \ =
6
6
2
1
;
p
(4.22)
En el sistema de coordenada cilindricas la ecuacion de Navier-Stokes es:
Figura 4.4: Sistema de coordenada cilindrica y esferica

y
:
P(r,7, :)
:
7
r

y
:
P(r,,7)
7
r

'
+

'

'
+

'

+ +

'
+

'

:
: : : :
:
:
: :
r
:
r
:
r
r
r
r r r r r
r
:
r r
r
r

:
v v
r r
v
r r
v
:
P
:
v
v
v
r
v
r
v
v
t
v

:
v
r
v
r
v
r r
v
r
v
r r
v
P
r :
v
v
r
v v v
r
v
r
v
v
t
v

:
v v
r
v
r r
v
r
v
r r
v
r
P
:
v
v
r
v v
r
v
r
v
v
t
v
2
2
2
2
2 2
2
2
2
2 2
2
2 2 2
2
2
2
2 2
2
2 2 2
2
2
1 1
1
2 1 1
1
2 1 1
1
7
;
p 7
7
;
7 p 7
7 7
;
p 7
7
7
7 7 7 7 7
7 7 7 7 7 7
7
7 7
(4.23)
En el sistema de coordenada esIerica la ecuacion de Navier-Stokes es:

'
+

'

+ \ +

= +

'
+

'

+ \ +

= +

'
+

'

\ +

=
+

7
7
7
7 7

7

7

7
7

7
;
7 p
7

;
p
7
;
p

v
r
v
r r
v
v
P
r r
v v
r
v v
t
Dv

v
r r
v
v
r
v
P
r r
v
r
v v
t
Dv

v
r r
v v
r r
v
v
r
P
r
v v
Dt
Dv
r
r
r
r
r
r
r
r
2 2 2 2 2
2
2 2 2 2 2
2
2
2 2 2 2
2
2 2
sen
cos 2
sen
2
sen
sen
1
tg
sen
cos 2
sen
2
1
tg
sen
2
tg
2
2 2
1
(4.24)
Donde

7
7

=
sen r
v
r
v
r
v
t Dt
D
r


2
2
2 2 2
2
2
2
sen
1
sen
sen
1 1

'
+

'

+
'
+

'

= \
r r r
r
r r

Ecuacin de Energia (Primera Ley de la 1ermodinmica)
La primera ley de la termodinamica para un Iluido cuya conductividad termica se puede
considerar isotropica se expresa de la siguiente Iorma:


7
' ) ( q
Dt
DP
T T k
Dt
DT
c
P
+ 1 + + \ \ = 3 . p (4.25)
Donde . es el coeIiciente volumetrico de expansion termica, 3 la viscosidad dinamica y 31 es
disipacion viscosa, que se obtiene por la siguiente expresion.

2
2
2
2 2
2 2
3
2
2

'
+

'

'
+

'

'
+

'

+
'
+

'

+
'
+

'

'
+

'

'
+

'

= 1
:
w
y
v

u
:
w
y
v

u
:
u

w
y
w
:
v

v
y
u
(4.26)
Los primeros tres terminos del lado derecho de la ecuacion (4.26), al ser multiplicados por 3
representan los esIuerzos cortantes viscosos y el resto a los esIuerzos normales. olectivamente,
los terminos representan la velocidad a la que la energia cinetica se convierte en Iorma
irreversible a energia termica debido a los eIectos viscosos.
Para el caso de un as ideal,
T
1
= . por lo que la ecuacion se simpliIica a
' ) ( q
Dt
DP
T k
Dt
DT
c
P
+ 1 + + \ \ = 3 p (4.27)
y para un liquido incompresible ( 0 = . ) y la ecuacion sera:
' ) ( q T k
Dt
DT
c
P
+ 1 + \ \ = 3 p (4.28)
La mayoria de los problema de conveccion, las propiedades del Iluido se pueden considerar
constantes, no existe generacion de calor, se desprecia la disipacion viscosa y los eIectos de
compresibilidad del Iluido, por lo que la primera ley se puede expresar de la siguiente Iorma:
T
Dt
DT
2
1
\ =

(4.2)
o que es lo mismo

2
2
2
2
2
2
1
:
T
y
T

T
:
T
w
y
T
v

T
u
t
T

'
+

'

(4.30)
En coordenadas cilindrica esta ultima ecuacion queda de la siguiente Iorma

'
+

'

+
'
+

'

+
'
+

'

=
'
+

'

:
T
:
T
r r
T
r
r r :
T
v
T
r
v
r
T
v
t
T
: r
7 7 7
7
2
1 1 1
(4.31)
En coordenadas esIericas la ecuacion (4.30) se expresa de la siguiente Iorma


8

'
+

'

'
+

'

+
'
+

'

'
+

'

7 7

7
7

T
r
T
r r
T
r
r r
T
r
v T
r
v
r
T
v
t
T
r
2 2 2
2
2
sen
1
sen
sen
1 1
sen
1
(4.32)
uando se tiene un Iluido cuya viscosidad es elevada, tal es el caso de los Iluidos lubricantes, la
disipacion viscosa no se puede despreciar. Por tanto, la primera ley queda de la siguiente Iorma:
1 + \ \ = 3 p ) ( T k
Dt
DT
c
P
(4.33)
En coordenadas cilindricas la disipacion viscosa se expresa por


2
2 2 2
2 2 2
3
2 1
2
1 1
'
:
v
r
v v
r r
v
r
v
:
v
:
v v
r
v
r r
v
r
v
: r r
: r : r
6
\

'
+

'

+
'
+

'

+
'
+

'

'
+

'

+
'
+

'

+
'
+

'

= 1
7
7 7
7
7 7 7
(4.34)
Donde
k
: r
r
r

1

= \ 7
7

k v v r v '
: r

+ + = 7
7
6

En coordenadas esIerica

2
2 2
2
2 2 2
3
2
tg sen
1 1
2
sen
1
sen
1
sen
sen 1
'
r
v
r
v v
r r
v
v
r r
v
r
v
r
r
v
r
v
r r
v
r
v
r r
v
r
r
r r r
r r
6
\

'
+

'

+ +

'
+

'

+
'
+

'

'
+

'

'
+

'

'
+

'

'
+

'

'
+

'

'
+

'

= 1
7
7 7

(4.35)
Donde
7
7

sen
1

= \
r r
r
r

7
7

v v r v '
r
+ + =
6

La Segunda Ley de la 1ermodinmica
%odo proceso de transIormacion de energia debe cumplir con la segunda ley de termodinamica, la
cual se expresa de la siguiente Iorma:


9
0 ) (
2
2
K 1 + \
T
T
T
k 3
(4.36)
Donde la temperatura en la ecuacion (4.36) debe ser expresada en escala absoluta (temperatura
absoluta).
CONSIDERACIONES ESPECIALES PARA ANLISIS DE LA CAPA LIMITE
Las ecuaciones mostradas en la seccion 4.1.3 presentan una explicacion completa de los procesos
Iisicos que determinan la capa limite hidrodinamica y termica. Sin embargo, esta ecuaciones se
pueden simpliIicar al considerar algunos aspectos de la capa limite, como son:
- El Ilujo es estacionario
- El Iluido se puede considerar incompresible (densidad constante)
- Las propiedades son constantes
- El problema se puede trata en Iorma bidimensional
- El espesor de la capa limite hidrodinamica es muy pequeo por lo que se cumple que:
v u >> ,

u
y
u

>>

,
y
v
y
u

>>

v
y
u

>>

, 0 =

y
P

- El espesor de la capa limite termica tambien en muy pequeo por lo que se cumple que

T
y
T

>>

Luego aplicar estas consideraciones, y tomado en cuanta que la capa limite se Iorma a lo largo
del eje , se tiene que:
0 =

y
v

u
(Ecuacion de continuidad) (4.37)

y
u
d
dP
y
u
v

u
u +

+ =

2
2
1
;
p
(Ecuacion de momentum) (4.38)

2
2
2

'
+

'

y
u
c y
T
y
T
v

T
u
P
v
(Ecuacion de eneria) (4.3)
Estas tres ecuaciones se conocen como ecuaciones de la capa limite, cuyas condiciones de borde
requeridas para resolverlas son las siguientes:

=
=
=
=
s
T T
v
u
y
) 0 , (
0 ) 0 , (
0 ) 0 , (
0 Condiciones de pared


10

=
=
=

T T
v
u u
y
) , (
0 ) , (
) , (
Condiciones de corriente libre
Las ecuaciones (4.37), (4.38) y (4.3) se adimensionan, haciendo el siguiente cambio de variable:

s
s
T T
T T
'
P
P
U
v
v
U
u
u

y
y

= = = = = =

7
p
,

, , , ,
2

Donde es una longitud caracteristica de la superIicie y U es la velocidad contracorriente de la
superIicie.
Al aplicar los cambios de variables y despreciando las Iuerzas de cuerpo, se tiene que
0

y
v

u
(Ecuacion de continuidad) (4.40)

2
2

#e
1

y
u
d
P d
y
u
v

u
u

+ =

(Ecuacion de momentum) (4.41)

Pr #e
1

'
+

'

y
u
y y
v

7 7
(Ecuacion de eneria) (4.42)
Donde:

v
'

= #e es el numero de #eynolds

v
= Pr es el numero de Prandtl
Las condiciones de borde para resolver las ecuaciones (4.40), (4.41) y (4.42), son:

=
=
=
=
0 ) 0 , (
0 ) 0 , (
0 ) 0 , (
0

v
u
y
7
Condiciones de pared

=
=
=

1 ) ,

(
0 ) ,

) , (

v
U
u
u
y
7
Condiciones de corriente libre
El coeIiciente de conveccion se determina por


11

0

=
y
f
y
k
h
7
(4.43)
y el numero de Nusselt

0

= =
y
f

y k
h
Nu
7
(4.44)
Ecuaciones de la capa limite para flujo turbulento
En la Figura 4.5 se muestra el comportamiento en el
tiempo de la velocidad del Iluido (componente x de
la velocidad) en un punto. Se puede notar que la
velocidad es Iluctuante en cada punto del Iluido, al
igual que el resto de las propiedades del Iluido
(densidad, temperatura, presion, etc.), por lo que se
debe expresar dichas propiedades de la siguiente
manera:
' u u u + = (4.45)
Donde u es el valor promedio de la velocidad en
un intervalo de tiempo y ' u es la Iluctuacion de la velocidad en torno a u . El promedio de la
velocidad (o de otra propiedad) se determina por

=
T
udt
T
u
0
1
(4.46)
Donde T es el periodo donde la Iluctuacion de la velocidad se repite.
El resto de las propiedades se pueden expresar de la siguiente Iorma:
' , ' , ' P P P T T T v v v + = + = + = (4.47)
Antes de introducir estos conceptos a las ecuaciones de continuidad, momento y energia se
presenta las reglas generales que esta notacion presenta. Estas reglas son:

0 , 0 ,
' , ' '
0 ' ,
2 2 2
=
'
+

'

=
'
+

'

+ = + =
= + = +
t
u
t
u

u
u u u v u v u uv
u u v u v u
(4.48)
Finalmente las ecuaciones para Ilujo turbulento en la capa limite son:
0 =

y
v

u
(Ecuacion de continuidad) (4.4)
Figura 4.5: Comportamiento en el
tiempo de la velocidad u en un punto
para un flufo turbulento
u
t
u
u


12

y
v u
y
u
d
P d
y
u
v

u
u +

+ =

) ' ' ( 1
2
2
;
p
(Ecuacion de momentum) (4.50)

y
T v
y
T
y
T
v

T
u

) ' ' (
2
2
(Ecuacion de eneria) (4.51)