Unidad Tematica 4

UNIVERSIDAD AUTONOMA GABRIEL RENE
MORENO
CARRERA DE INGENIERIA DE ALIMENTOS
TEMA 1 TRANSFERENCIA DE
CALOR POR CONDUCCION
TRANSPORTE EN INTERFASE
BALANCES MACROSCÓPICOS ISOTÉRMICOS
“La presente es una recopilación de diversas notas y apuntes de cátedra dispersas que
fueron elaboradas o redactadas en los últimos años. Se agradece especialmente la
colaboración de los alumnos cursantes en 2008: Argüelles,Juan; Lopez, Carlos; Lopez,
Yanina; Navarro, Mauro; Pirani, Paula; Tejerina, Néstor; Bodo, Germán; Corrente,
Pablo; Piccinelli, Ariel; Zanzo, Juan Carlos; Regner, Leandro para su compilación y
organización, que con la coordinación del Auxiliar Juan M. Dominguez permite disponer
de esta versión revisada. Se advierte que estas notas son solo una guía para el
estudio, debiendo consultarse la bibliografía recomendada en cada tema para lograr un
conocimiento pleno de los mismos.”
Página 2 de 24
Fenómenos de Transporte – Unidad Temática 4
Revisión: Octubre 2008
4.A. TRANSPORTE EN INTERFASE
Las soluciones analíticas a los problemas de flujo laminar en geometrías sencillas pueden
ser halladas, al igual que si se conocen las relaciones empíricas entre la densidad de flujo de
cantidad de movimiento y los gradientes de velocidad de tiempo ajustado para el flujo
turbulento.
Para flujo en conducciones se pretende obtener la relación entre caída de presión y el
caudal; para circulación alrededor de objetos sumergidos lo que interesa es la relación entre la
velocidad de aproximación del fluido y la fuerza resistente.
En muchos casos de ambos tipos no se conoce (ni se puede calcular) la distribución de
velocidades por lo que la solución analítica no es posible.
Para obtener la respuesta se utilizan datos experimentales junto con el análisis
dimensional construyendo gráficos y “correlaciones” que permiten estimar el comportamiento
del flujo en sistemas geométricamente semejantes.
4.A.1. FLUJO EN CONDUCTOS CILINDRICOS

Para el flujo estacionario de un fluido de densidad constante que circula por un conducto
rectilíneo de sección uniforme, sobre la superficie sólida ejerce una fuerza Fk que depende del
comportamiento cinético y que tiene la misma dirección que la velocidad media <v> en el
conducto.
Esta fuerza cinética se puede expresar arbitrariamente como,
Fk = A ⋅ Kv ⋅ f
Siendo;
Fk : la fuerza cinética ejercida sobre el sólido
A : un área característica mojada por el fluido
Kv : una energía cinética característica por unidad de volumen
f : un número adimensional, el “factor de fricción”
Para un conducto cilíndrico de diámetro D, en la longitud L, el factor de fricción se puede definir por
las siguientes,
A = π DL = 2π RL r p0
z
1 2
Kv = ρ vz <v>
2 L
1 2
Fk = π DL ρ v z f Fk
2
pL
D
Página 3 de 24
Aplicando un balance de fuerzas al fluido entre h0 y hL, considerando flujo totalmente
desarrollado (dirección “z”),
π D2 π D2
Fk = ( p0 − pL ) + ρ g (h0 − hL )
4 4
y para P = p + ρ gh
π D2
Fk = (P0 − PL )
4
Igualando las ecuaciones de Fk;
π D2 1 2
(P0 − PL ) = π DL ρ v z f
4 2
2
(P − PL ) D 1 2f ρ L v z
obteniendo, f = 0 2 ó ( P0 − PL ) =
ρ vz L 2 D
que permite calcular f de datos experimentales, siendo la expresión la ya conocida Ecuación

de FANNING.
CORRELACIONES PARA f
Para el flujo estacionario, incompresible, isotérmico y desarrollado de un fluido en una
tubería lisa, la fuerza ejercida sobre las paredes, tanto en régimen laminar como turbulento
será,
Fk = ∫τ
sup
rz pared dA
siendo,
∂v z
τ rz = −µ
pared
∂r r =R
dA = Rdϕ dz
2π L
A= ∫
sup
dA = ∫ ∫ Rdϕdz
0 0
por lo tanto,
L 2π
∂v z
Fk = ∫ ∫ ( − µ )Rdϕ dz
0 0
∂r r =R
y de la ecuación de definición de f ;
Página 4 de 24
L 2π
∂v z
∫ ∫ (−µ ∂r r =R
)Rdϕ dz
f = 0 0
1 2
π DL ρ v z
2
Definiendo las variables adimensionales,
z r v p − p0 Dρ v
z* = ; r * = ; v z * = z ; p* = ; Re = ;
D D vz ρ vz µ
y reemplazando,
L
D 2π
∂ vz vz *
∫ ∫ (−µ
0 0
)
∂Dr * r * = 1
Rdϕ dDz *
f = 2
1 2
π DL ρ v z
2
4 D 2π
∂v z *
µ v z RD ∫ ∫ (− ∂r *
)dϕ dz *
1
0 0 r *=
f = 2
1 2
π DL ρ v z
2
4 D 2π
1D 1 ∂v z *
f =
π L Re ∫ ∫ (−
0 0
∂r * r * = 1
)dϕ dz *
2
Como para flujo desarrollado el perfil de velocidades es independiente de la posición axial

(z) y por la geometría cilíndrica también lo es del ángulo (φ) el gradiente de velocidad
adimensional en la pared dependerá solamente del Re.
Esto significa que no se introducen nuevas variables, y que el factor de fricción (f) será,
f = f (Re)
En caso que el perfil de velocidades no sea completamente desarrollado, puede interesar
∂v z *
la distancia desde la embocadura, y en ese caso, = f ( z*) quedando,
∂r *
L
f = f (Re, )
D
No obstante, si la longitud de entrada es pequeña respecto a L, no se considera la

dependencia con L/D y se utiliza la forma aproximada de f = f (Re) .
Página 5 de 24
Flujo Laminar
La ecuación de Fanning se puede combinar con la de Hagen-Poiseuille cuando Re <
2100, obteniéndose que
16
f =
Re
Flujo Turbulento
En este caso se observa que la caída de presión depende del estado de la superficie de
interfase. Para superficies limpias y pulidas, como porcelana, vidrio o plásticos se observa que
para un caudal determinado la caída de presión necesaria es menor que para otros materiales
como fundiciones, aceros, cemento, etc.
Los primeros se llaman “tubos lisos”, y algunas expresiones son,
0.0791
f = Fórmula de Blasius (válida para 2100 < Re < 100.000)
Re0.25
También se han propuesto,

f = 0.046Re −0.2 Ecuación de Von Karman
1
= 4log10 Re f − 0.40 Ecuación de Nikuradse-Prandtl (3000 < Re < 100.000)
f
Para tubos “rugosos”, es decir, aquellos en que el estado de la pared opone una
resistencia adicional al flujo, se define un parámetro adimensional llamado “rugosidad relativa”
(ξ/D) que es la relación entre la altura de la protuberancia o irregularidad frente al diámetro del
tubo.
Algunas formas irregulares se muestran, y depende del material, el método de
construcción del tubo, la antigüedad en servicio, etc.
ξ ξ
ξ ξ
L ξ
En general, f = f (Re, , )
D D
Algunas expresiones clásicas son;

1 R
= 4.06 log + 3.36 Ecuación de Von Karman
f ξ
Página 6 de 24
1 D
= 4.00 log + 2.28 Ecuación de Nikuradse
f ξ
1
válida para > 0.01. Esta fórmula es válida para flujo completamente desarrollado,
ξ
Re f
D
donde f se hace constante.
Otra es;
1 D D /ξ
= 4.0log + 2.28 − 4.0log(1 + 4.67 )
f ξ Re f
que es la Ecuación de Colebrook válida para Re >3000, y que tiene en cuenta la variación de f
con Re, siendo una de las más aproximadas a los resultados experimentales.
En todos los casos, experimentalmente se ha graficado el valor de f en función de Re

para la región laminar y turbulenta, en coordenadas logarítmicas, obteniendo un gráfico muy
difundido que se conoce como “diagrama de Moody” (ver figura en hoja en hoja siguiente).
Página 7 de 24
4.A.2. FLUJO EN CONDUCTOS NO CIRCULARES

Se define un “radio hidráulico” como la razón entre la superficie de flujo real y el perímetro
mojado,
Área Flujo
Rh =
Perímetro Mojado
Como para un tubo de diámetro D,
π D2
Rh = 4 =D
πD 4
se define,
4 ⋅ Rh ⋅ ρ v z
Re =
µ
Página 8 de 24
4.A.2. FLUJO EN CONDUCTOS NO CIRCULARES
Se define un “radio hidráulico” como la razón entre la superficie de flujo real y el perímetro
mojado,
Área Flujo
Rh =
Perímetro Mojado
Como para un tubo de diámetro D,
π D2
Rh = 4 =D
πD 4
se define,
4 ⋅ Rh ⋅ ρ v z
Re =
µ
y la ecuación de Fanning se puede escribir,

2
(P0 − PL ) ⋅ 2 ⋅ Rh f ρ vz L
f = ó P0 − PL =
ρ vz L
2
2 ⋅ Rh
Métodos de Cálculo:
Para tuberías puede ser necesario determinar la caída de presión en un sistema, el
caudal a circular para un -∆P dado, o el diámetro necesario con un caudal y caída de presión
definidos.
Como f = f (Re) ; exceptuando el caso del cálculo de -∆P, en los otros habrá que usar
métodos de tanteo iterativos o por resolución gráfica (del diagrama de Moody).
En los casos de solución numérica (no gráfica) se recomienda la ecuación de Colebrook
como la más aproximada a los resultados experimentales en todo el campo de aplicación (flujo
no-desarrollado o totalmente conseguido el perfil).
4.A.3. FLUJO ALREDEDOR DE OBJETOS SUMERGIDOS

En este caso las líneas de corriente sufren una deformación alrededor del objeto, de tal
forma que el vector velocidad del fluido a gran distancia del mismo forma distintos ángulos con
distintos puntos de la superficie de éste.
El rozamiento entre fluidos y sólidos es originado pos dos tipos de fuerzas:
1- Fuerzas tangenciales originadas por el esfuerzo cortante de la superficie (fricción de
superficie).
2- Fuerzas normales a la superficie originadas por la presión (fricción de forma).
Página 9 de 24
Los componentes en la dirección del flujo originan el rozamiento que es la suma de la

fricción de superficie y la de forma, a diferencia, del flujo en conductos donde sólo existía
fricción de superficie, porque las fuerzas de presión en la superficie del conducto no tenían
componentes en la dirección del flujo (perpendiculares en todos los puntos).
Fuerzas sobre el sólido Componentes de la resistencia al flujo en dirección al movimiento
Para el caso de una esfera que cae en el seno de un fluido por acción gravitatoria la
fuerza resultante de caída es igual a la diferencia entre el peso y el empuje. A esta fuerza
resultante se le opone una fuerza resistente (Fk) originada por el rozamiento. A medida que la
velocidad de caída aumenta, ésta también lo hace de forma que en un instante dado la
resultante es nula y la velocidad se hace constante. Esta se conoce como “velocidad limite de
sedimentación”.
La deformación que origina el sólido en las líneas de corriente perturba el flujo a gran
distancia por lo cual la velocidad se considera cuando no hay otras superficies sólidas en las
cercanías y en general para el movimiento relativo sólido-fluido se llama “velocidad de
aproximación” ( v ∞ ).
Definiendo la fuerza resistente de igual forma que antes se puede escribir:
Fk = f ⋅ Acaract ⋅ K v
donde el área característica es la proyección del sólido sobre un plano perpendicular al
movimiento del fluido no perturbado, que en el caso de esferas es el área de un círculo, de
forma de escribir:
1
A = π R2 Kv = ρv ∞ 2
2
⎛1 ⎞
( ⎝2
)
Fk = π R 2 ⎜ ρv ∞ 2 ⎟ f
⎠
(1)
Página 10 de 24
Cuando se alcanza la velocidad límite de sedimentación tendremos:

4 4
Peso = π R 3 ρs g Empuje = π R3 ρg
3 3
Peso – Empuje = Fk
4
Fk = π R 3 ( ρs − ρ )g (2)
3
Igualando las ecuaciones (1) y (2) se obtiene el coeficiente de frotamiento, también

llamado coeficiente de resistencia (cf) o coeficiente de arrastre (cd):
4 Dg ρs − ρ
f = ( )
3 v ∞2 ρ
Este coeficiente representa las pérdidas por fricción de superficie y de forma siguiente la
misma ley general que el coeficiente de fricción de Fanning, es decir que básicamente
representa lo mismo, lo cual se comprueba representándolo en función del número de
Reynolds definido por la siguiente ecuación;
ρv ∞ D
NRe =
µ
Stokes, obtuvo una expresión de la fuerza resistente para velocidades muy bajas
representadas en la ecuación:
Fk = 6πµ Rv ∞
e igualando ésta con la ecuación de definición (1) se obtiene:

24
f = (3)
Re
Que se cumple exactamente para Reynolds menores de 0,1; para Reynolds mayores la
variación de f con el Reynolds se ha determinado experimentalmente destacándose de que no
hay discontinuidad en la curva hasta valores de Reynolds muy elevados aproximadamente
mayores a 105.
Hacemos notar que el coeficiente de arrastre de la ecuación (3) puede estar representado
por una constante numérica distinta, que depende de la definición de Fk que se adopte.
12
f =
Re
En la región que va de Reynolds 500 a 200.000 la fuerza resistente es aproximadamente

proporcional al cuadrado de la velocidad de aproximación obteniendo valores f ≅ 0.44 .
Página 11 de 24
La caída brusca en la cercanía de Reynolds 200.000 se debe a la formación de un

conjunto dinámico homogéneo entre la partícula y la estela que lo acompaña.
Para el caso de sólidos distintos de las esferas se define un coeficiente de corrección
denominado “esfericidad” que representa la relación:
Área sup. esfera de igual volumen

ϕ=
Área sup. partícula
Factor de fricción f
18,5
f ≅
Re3/5
Asíntota Ley de Stokes f =

24
Re
Ley de Stokes Ley Intermedia Ley de Newton
Número de Reynolds
Página 12 de 24
4.B. BALANCES MACROSCÓPICOS ISOTÉRMICOS
Dado que la resolución de las ecuaciones generales para los balances de materia
(ecuación de continuidad) y de cantidad de movimiento (ecuación de movimiento) son
comúnmente difíciles de alcanzar (excepto sistemas sencillos), en la práctica ingenieril se
prefiere utilizar estos balances a sistemas “macroscópicos”, con entrada y salida de fluidos, en
los cuales se omiten aquellos términos de las mismas ecuaciones que resultan despreciables
en un problema dado.
4.B.1. BALANCE DE MATERIA
Calor (Q)
Sección (2)
Sistema
Sección (1)
Trabajo
En el sistema de la figura existe una entrada de fluido en la seccion (1) y una salida en
(2); ademas puede existir transmisión de calor (tomada como positiva si ingresa) y además
puede realizar trabajo (tomado como positivo si produce trabajo externo), y además tendremos
en cuenta que;
a) En las secciones (1) y (2) la velocidad alisada de tiempo ajustado temporalmente es
paralela a las paredes del conducto.
b) Todas las propiedades fisicas (µ, ρ) no varian en la seccion transversal S1 y S2.
De aquí:
dmt
= ρ1 < v1 > S1 − ρ2 < v 2 > S2 (1)
dθ
Donde:
mt es la masa total del fluido contenida en el sistema (mt = ∫ ρdv )
v
< v i > : es el promedio de velocidades medias alisadas temporalmente.
Página 13 de 24
Teniendo en cuenta que:

w = ρ <v >S
∆w = w 2 − w1
dmt
= −∆w
dθ
Balance macroscópico general de materia (2)
Para estado estacionario, la masa de fluido no se acumula en el sistema, luego:

dmt
=0
dθ
∆w = 0
Balance de materia para estado estacionario (3)

w1 = w 2
Es menester aclarar que la velocidad tal como aparece en las ecuaciones anteriores
depende de la naturaleza del flujo, de forma que se debe evaluar el perfil de velocidad en las
secciones transversales y obtener la velocidad media por las relaciones de flujo potencial,
laminar o turbulento, o aproximarla.
4.B.2. BALANCE DE CANTIDAD DE MOVIMIENTO
El principio de la cantidad de movimiento (deducido de la 2° Ley de Newton) expresa que

“la suma de todas las fuerzas externas que actúa sobre una masa de fluido es igual a la
rapidez de variación de la cantidad de movimiento respecto al tiempo”.
uur
ur dM d ∫ ρ vdV
F= t
= v (4)
dθ dθ
Mt: Cantidad de movimiento total

F: Fuerzas externas
Y donde las fuerzas externas pueden ser:
1) Fuerzas normales a las fronteras (tales como fuerzas de presión).
Página 14 de 24
2) Fuerzas paralelas a las fronteras (tales como esfuerzos tangenciales).

3) Fuerzas de campo (tales como campos magnéticos o gravitacionales).
De esta forma:
uuur
dMT
= ρ1 < v12 > S1 − ρ 2 < v 22 > S2 + p1S1 − p2S2 − F + mT g (5)
dθ
Donde:
ρ i < v i2 > Si : es la velocidad de cantidad de movimiento debido al movimiento global

del fluido en la entrada y salida.
pi Si : son las fuerzas de presión en los extremos del sistema.
F : Fuerzas del fluido sobre el sólido (viscosidad, presión, rozamiento, son salidas o
disipaciones de cantidad de movimiento, de aquí negativas).
mT g : Acción gravitatoria, peso del fluido contenido en el sistema.
uuur uur ur ur
La anterior es una ecuación vectorial (vectores MT ,Si , F , g ), donde la dirección y sentido
de Si corresponden a la velocidad alisada en dicha sección. Introduciendo los vectores:

uur uur
w1 = ρ1 < v1 > S1 w 2 = ρ2 < v 2 > S2
La ecuación anterior quedara:
dMT ⎡< v2 > ⎤

= −∆ ⎢ w + pS ⎥ − F + mt g
dθ ⎣⎢ < v > ⎦⎥
Balance macroscópico general de la cantidad de movimiento (6)
Para el estado estacionario, la cantidad de movimiento no se acumula en el sistema y se

obtiene:
⎡< v2 > ⎤
F = −∆ ⎢ w + pS ⎥ + mt g
⎣⎢ < v > ⎦⎥
Balance de cantidad de movimiento para estado estacionario (7)
Generalmente, en el flujo turbulento (perfil de velocidad plano) se sustituye < v >

2
por
<v >
esto < v > , ya que el error cometido oscila entre el 2% y el 5%.
Página 15 de 24
Este resultado es útil para calcular las fuerzas que actúan sobre piezas o aparatos tales
como paletas de turbinas y curvas de tuberías; también es útil, junto con otros balances
macroscópicos para describir el funcionamiento de eyectores de vapor de agua y otros
aparatos en los que hay mezcla de corrientes de fluidos.
4.B.3. BALANCE DE ENERGÍA MECÁNICA

Plantearemos este caso para sistemas en estado estacionario, es decir que, cumplan las
siguientes condiciones:
a) El fluido que ingresa tenga propiedades y velocidades invariables en el tiempo.
b) Igual para el egreso aunque estas puedan ser distintas en el ingreso.
c) Las propiedades físicas del fluido dentro del sistema son constantes en el tiempo.
d) La masa contenida es constante (entradas = salidas).
e) El calor suministrado y trabajo entregado por el sistema son constantes.
El principio de conservación de la energía indica que para la energía total que ingresa
debe ser igual a la energía total que sale; para el análisis se omite considerar efectos
eléctricos, magnéticos o químicos.
En el sistema tendremos en cuenta:
A. Energía transportada por el fluido.
B. Energía transferida entre el fluido y los alrededores del sistema.
A. Entre estas tenemos:

A.1.- Energía interna (U) que es una propiedad del fluido.
A.2.- Energía potencial externa ( Φ ) debida a la posición del fluido respecto a un plano
de referencia.
A.3.- Energía cinética externa (K) asociada al movimiento del fluido.
A.4.- La energía de presión (pV) transportada por el fluido al introducirse al sistema.
B. En estas incluimos:
B.1.- El calor absorbido por el sistema (Q’), que puede o no variar la temperatura del
fluido. Debe pasar las fronteras del sistema.
B.2.- El trabajo ( τ ) transferido del fluido al exterior por acción de algún equipo (por
ejemplo una turbina). Debe pasar las fronteras del sistema.
Por lo tanto, el balance de todas las energías incluidas en el sistema de flujo,

agrupando energía de entrada al sistema (subíndice 1) y las de salidas del sistema
(subíndice 2), de aplicación rigurosa en estado estacionario será:
Página 16 de 24
w1 [U1 + Φ1 + K1 + pV
1 1 + Q´] = w 2 [U 2 + Φ 2 + K 2 + p2V2 + τ ] (8)
También se puede escribir referido a la unidad de masa del sistema ( $ )

∧ ∧ ∧ ∧ ∧ ∧ ∧ ∧ ∧ ∧
U 1 + Φ1 + K 1 + p1 V1 + Q' = U 2 + Φ 2 + K 2 + p2 V 2 + τ (9)
O en forma de incrementos ( ∆ = salidas - entradas):
∧ ∧ ∧ ∧ ∧ ∧
∆ U + ∆ Φ+ ∆ K + ∆( pV ) = Q'− τ (10)
Para los fluidos reales, como resultado del flujo se genera una fricción que convierte la
energía mecánica irreversiblemente en energía calorífica, de tal forma que:
1) O bien, el calor real que atraviesa las fronteras del sistema será:
∧ ∧ ∧
Q = Q + Ev (11)
2) O bien, todo el trabajo realizado por el fluido no se transfiere íntegramente al exterior.
∧ ∧ ∧
τ = τ + Ev (12)
Cualquiera de estas expresiones anteriores representa lo mismo, “las pérdidas por
∧
fricción” ( Ev ) irreversibles y conduce al mismo resultado.
Por lo tanto, en la ecuación expresada, como incrementos:

∧ ∧ ∧
⎛ ∧⎞ ∧ ∧ ∧
∆ U + ∆ Φ+ ∆ K + ∆ ⎜ pv ⎟ = Q − Ev − τ (13)
⎝ ⎠
Donde:
∧
Q : cantidad de calor (por unidad de masas) que realmente atraviesa la frontera.
∧
τ : Trabajo (por unidad de masa) que realmente atraviesan las fronteras.
∧
Ev : Pérdidas por fricción (por unidad de masa) irreversibles dentro del sistema.
Del primer principio de la termodinámica:

∧
v2
∧ ∧ ∧ ∧ ∧
∆ U = Q− τ = Q− ∫ p ⋅ d v (14)
∧
v1
Y teniendo en cuenta:
Página 17 de 24
∧
v2 p2
⎛ ⎞ ∧ ∧ ∧
∆ ⎜ pv ⎟ =
⎝ ⎠
∫ pd v + ∫ v dp
∧ p1
(15)
v1
Se obtiene de la ecuación (14):

p2
∧ ∧ ∧
⎛ ∧⎞
∆ U = Q + ∫ v dp − ∆ ⎜ pv ⎟ (16)
p1 ⎝ ⎠
Reemplazando en la ecuación (13) y simplificando:
p2
∧ ∧ ∧ ∧ ∧
∆ Φ+ ∆ K + ∫ v dp = − Ev − τ (17)
p1
O también:
p2
∧ ∧
1 ∧ ∧
∆ Φ+ ∆ K + ∫ρ
p1
dp + Ev = − τ
Balance macroscópico de energía mecánica en estado estacionario (18)
Es de notar, que cualquiera de las ecuaciones 8, 9, 10, 13, 17 ó 18 representan la ley de

conservación de la energía, simplemente que estas últimas (17,18) se expresan en unidades
mecánicas.
Otras formas comunes del balance energía, introducen el término de entalpía:
∧ ∧ ∧ ∧ ∧
∆ Φ+ ∆ K + ∆ H = Q − τ (19)
∧
Donde ∆ H es la diferencia de entalpías, que como es un diferencial total depende sólo
de valores iniciales y finales del fluido, independientemente del camino seguido por el mismo.
También se puede expresar en función de la energía libre de Gibbs, tal como:
∧ ∧ ∧ ∧ ∧
∆ Φ+ ∆ K + ∆ G ibbs + Ev = − τ (20)
Donde para un sistema isotérmico:

p2
∧
1
∆ G ibbs = ∫ ρ dp
p1
Para fluidos ideales, sin fricción, que no intercambia trabajo mecánico con el medio
∧ ∧
ambiente, los términos τ y Ev son nulos, y queda el balance expresado como la ecuación de
Bernoulli:
Página 18 de 24
p2
∧ ∧
1
∆ Φ+ ∆ K + ∫ ρ dp = 0
p1
Ecuación de Bernoulli (21)
∧
El término ∆ K (energía cinética por unidad de masa), se representará por una expresión
del tipo:
∧
1 1 〈v 3 〉
(velocidad ) = ∆
2
∆K = ∆
2 2
〈v 〉
Como en la sección transversal de flujo la velocidad varía, en una sección transversal

(dS), tendremos un gasto:
dω = ρ i v i dSi donde v i ⊥ Si
vi 2
y cada unidad de masa de fluido que atraviesa (dS) tiene una energía cinética:
2
luego, el flujo de energía cinética por unidad de tiempo en el área (dSi) será:
vi 2 1
dK i' = ( ρ i v i dSi ) = ρ i v i3dSi
2 2
donde K’ es energía/tiempo.
Luego, en toda la sección transversal, (S) para un fluido incompresible la energía cinética
será:
1
K' = ρ ∫ v i3dSi
2 sup
Y el gasto:
ω = ρ ∫ v i dSi
sup
Luego, la energía cinética por unidad de masa:
1 1
∧ ρ ∫ v i3dSi
K' ∫ v i3dSi
= 2 = 2
s s
K=
ω ρ ∫ v i dSi v S
s
Esta relación aparece a veces bajo la forma de un factor de corrección de la energía

cinética ( α ), tal que:
Página 19 de 24
∫v
3
i dSi
α= s
3
v S
2
∧ α v
K= (22)
2
Para evaluar α debe conocerse la variación de velocidad local en la sección transversal.
Así, para flujo laminar, donde < v >= 0,5 ⋅v max.
τ yx
ds = 2π rdr v=
2Rµ
(R 2
−r2 )
τ yx R
< v >=
4µ
∧
Se obtiene: α = 2 y K lam =< v > 2
Para el flujo turbulento, donde el perfil de velocidad es casi plano;

∧
< v >2
α ≅1 K turbulento = (23)
2
^
El término ∆ Φ : la energía potencial por unidad de masa, se representa por:
∧
∆ Φ = g ∆h = g ( h2 − h1 ) (24)
Para los casos habituales de aceleración gravitatoria (g) constante, siendo h1 y h2 las
alturas sobre un mismo nivel de referencia.
El termino ∆ ω isotérmico: la energía de presión para el sistema isotérmico se evalúa

según sea la variación de la densidad con la presión. Para fluidos incompresibles, ( ρ =
constante)
p2
1 1
∫ ρ dp = ρ ( p
p1
2 − p1 ) (25)
Para gases ideales:

p2 p2
1 RT RT p2 p1 p2
∫ ρ dp = ∫ Mp dp =
p1 p1
M
ln = ln
p1 ρ1 p1
(26)
Página 20 de 24
4.B.4. LAS PERDIDAS POR FRICCIÓN

Las pérdidas por energía mecánica por unidad de masas circulante debido a la fricción se
podrían calcular si se conocen todos los términos restantes de la ecuación de la energía
mecánica (18). No obstante lo común es estimar dichas pérdidas de acuerdo a las condiciones
de flujo conocidas. El desarrollo analítico dependerá de que el fluido circule por un conducto o
el fluido circule alrededor de un objeto sumergido.
a) Flujo en conductos: Para un fluido incompresible ( ρ =constante) que circula por un
conducto de sección circular transversal constante (S) en régimen estacionario originado en
fuerzas de presión y gravitatorias la ecuación del movimiento y de la energía se pueden escribir
de la siguiente forma:
F = ( p1 − p2 )S + ρ gLS (Movimiento) (27)
∧ p1 − p2
Ev = + gL (Energía) (28)
ρ
Multiplicando la ecuación (28) por ρ S y restando miembro a miembro de la ecuación del

movimiento, se obtiene:
∧ ( p1 − p2 )
E v ρS = ρ S + gL ρ S
ρ
∧ ∧
F
E v ρS − F = 0 Ev = (29)
ρS
La fuerza F es la fuerza de fricción que existe solamente cuando hay flujo, por lo cual en
lo sucesivo la denominaremos fuerza cinética (Fk).
De los balances envolventes para el flujo en conducto se determinó que esta fuerza
dependía de las diferencias de presiones y de la sección transversal de flujo:
Fk = (Po − PL )π R 2 (30)
Definiremos ahora a dicha fuerza como el producto de tres términos: un coeficiente de

frotamiento adimensional (f), un área característica de presión (A) y una energía cinética por
volumen de fluido (K), tal que:
Fk = f .A.Kv (31)
El flujo en conducto cumple con las características ya definidas;
A = 2π .R.L (área de fricción)
Página 21 de 24
1
Kv = ρ .〈v〉 2 (energía cinética/volumen)
2
Por lo tanto, podemos relacionar la ecuación (31) que no es una ley de la mecánica de
los fluidos, sino una definición de f, con la ecuación (30):
1
f .(2π .R.L ). ρ .〈v 〉 2 = (P0 − PL ).π .R 2 (33)
2
Obteniendo entonces f, en función de de la diferencia de presión, de las dimensiones de

conducto, de la densidad y velocidad del fluido:
(P0 − PL ).D
f = (34)
2.L.ρ.〈v 〉 2
2.f .L.ρ .〈v 〉 2

P0 − PL = (35)
D
Esta última da la pérdida de presión debida a la fricción del fluido y se conoce como
ecuación de Fanning, que ya se obtuvo. Alguna bibliografía utiliza un factor mayor, es decir el
aquí definido es cuatro veces más chico que aquel, pero presenta el mismo fenómeno y por lo
tanto, este coeficiente de fricción será también función del Reynolds y se podrá representar en
un gráfico (de f en función del Re) en coordenadas logarítmicas.
Página 22 de 24
A fin de generalizar la ecuación de Fanning para conductos de sección transversal no

circular se define un coeficiente llamado “radio hidráulico”, que representa la relación superficie
de flujo perímetro mojado, de tal forma que para conductos circulares se obtiene la siguiente
relación:
S π .R 2 R D
RH = = = = (36)
P 2π .R 2 4
Por lo tanto, la ecuación de Fanning quedará:
f .L.ρ .〈v 〉 2
P0 − PL = (37)
2.RH
donde f queda definido para un número Re en función del radio hidráulico, tal que:
4RH .〈v 〉.ρ

NRe = (38)
µ
A fin de relacionar la ecuación (37) con la ecuación de las pérdidas por fricción (29) la
podemos expresar de la siguiente forma:
P0 − PL 1 2 L
= 〈v 〉 .f .( ) (39)
ρ 2 RH
De la cual se obtiene que de las pérdidas de energía mecánica por unidad de masa de
fluido se pueden calcular del producto de la energía cinética del fluido por la siguiente
expresión:
∧
1 2 L
Ev = 〈v 〉 .f .( ) (40)
2 RH
Teniendo en cuenta que el cociente de las fuerzas de fricción sobre la sección transversal
de flujo representan la pérdida por fricción. La ecuación (V.40) se puede expresar como el
producto de una energía cinética por unidad de masa, por un “parámetro de fricción”
representado por la siguiente ecuación:
∧
1 2
Ev = 〈v 〉 .ev (41)
2
Página 23 de 24
El parámetro de fricción depende si consideramos flujo en conductos rectos donde las

líneas de corriente son paralelas o si consideramos flujo en sistemas donde éstas se desvían
por la existencia de perturbaciones tales como curvas, derivaciones, válvulas, codos, o demás
accesorios. En el primer caso ev es fácilmente determinable si se conocen las dimensiones del
conducto, el número de Reynolds y la rugosidad relativa, de acuerdo a la expresión (41).
Para el caso de perturbaciones ev se podrá determinar en base a la geometría del

sistema (ver problema de ensanchamiento brusco) o en casos más complicados por valores
experimentales, los cuales se representan por nomogramas de longitudes equivalentes, donde
la pérdida por fricción en el accesorio se equipara a la pérdida que ocasiona una determinada
longitud de tubería recta.
De aquí entonces las pérdidas por fricción en un sistema resultarán de la suma de dos
términos generales; uno de ellos la sumatoria de las pérdidas en las longitudes rectas y otro la
sumatoria de las longitudes equivalentes en los accesorios.
El balance de energía mecánica se puede expresar entonces de la siguiente forma:
p2
∧ ∧
1 ∧
∆ Φ+ ∆ K + ∫ ρ dp + τ + ∑ pérdidas por fricción en conductos +
p1
∑ pérdidas por fricción en accesorios =0
Página 24 de 24

Unidad Tematica 4

Cargado por

Copyright:

Formatos disponibles

Unidad Tematica 4

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Unidad Tematica 4

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD AUTONOMA GABRIEL RENE

CARRERA DE INGENIERIA DE ALIMENTOS

4.A. TRANSPORTE EN INTERFASE

4.A.1. FLUJO EN CONDUCTOS CILINDRICOS

que permite calcular f de datos experimentales, siendo la expresión la ya conocida Ecuación

Como para flujo desarrollado el perfil de velocidades es independiente de la posición axial

No obstante, si la longitud de entrada es pequeña respecto a L, no se considera la

También se han propuesto,

Algunas expresiones clásicas son;

En todos los casos, experimentalmente se ha graficado el valor de f en función de Re

4.A.2. FLUJO EN CONDUCTOS NO CIRCULARES

4.A.2. FLUJO EN CONDUCTOS NO CIRCULARES

y la ecuación de Fanning se puede escribir,

4.A.3. FLUJO ALREDEDOR DE OBJETOS SUMERGIDOS

Los componentes en la dirección del flujo originan el rozamiento que es la suma de la

Fuerzas sobre el sólido Componentes de la resistencia al flujo en dirección al movimiento

Definiendo la fuerza resistente de igual forma que antes se puede escribir:

Cuando se alcanza la velocidad límite de sedimentación tendremos:

Igualando las ecuaciones (1) y (2) se obtiene el coeficiente de frotamiento, también

e igualando ésta con la ecuación de definición (1) se obtiene:

En la región que va de Reynolds 500 a 200.000 la fuerza resistente es aproximadamente

La caída brusca en la cercanía de Reynolds 200.000 se debe a la formación de un

Área sup. esfera de igual volumen

Asíntota Ley de Stokes f =

Ley de Stokes Ley Intermedia Ley de Newton

4.B. BALANCES MACROSCÓPICOS ISOTÉRMICOS

4.B.1. BALANCE DE MATERIA

 < v i > : es el promedio de velocidades medias alisadas temporalmente.

Teniendo en cuenta que:

Balance macroscópico general de materia (2)

Para estado estacionario, la masa de fluido no se acumula en el sistema, luego:

Balance de materia para estado estacionario (3)

4.B.2. BALANCE DE CANTIDAD DE MOVIMIENTO

El principio de la cantidad de movimiento (deducido de la 2° Ley de Newton) expresa que

 Mt: Cantidad de movimiento total

Y donde las fuerzas externas pueden ser:

1) Fuerzas normales a las fronteras (tales como fuerzas de presión).

2) Fuerzas paralelas a las fronteras (tales como esfuerzos tangenciales).

 ρ i < v i2 > Si : es la velocidad de cantidad de movimiento debido al movimiento global

de Si corresponden a la velocidad alisada en dicha sección. Introduciendo los vectores:

La ecuación anterior quedara:

dMT ⎡< v2 > ⎤

Balance macroscópico general de la cantidad de movimiento (6)

Para el estado estacionario, la cantidad de movimiento no se acumula en el sistema y se

Balance de cantidad de movimiento para estado estacionario (7)

Generalmente, en el flujo turbulento (perfil de velocidad plano) se sustituye < v >

esto < v > , ya que el error cometido oscila entre el 2% y el 5%.

4.B.3. BALANCE DE ENERGÍA MECÁNICA

A. Entre estas tenemos:

Por lo tanto, el balance de todas las energías incluidas en el sistema de flujo,

También se puede escribir referido a la unidad de masa del sistema ( $ )

O en forma de incrementos ( ∆ = salidas - entradas):

Por lo tanto, en la ecuación expresada, como incrementos:

Del primer principio de la termodinámica:

Se obtiene de la ecuación (14):

Reemplazando en la ecuación (13) y simplificando:

Balance macroscópico de energía mecánica en estado estacionario (18)

Es de notar, que cualquiera de las ecuaciones 8, 9, 10, 13, 17 ó 18 representan la ley de

Donde para un sistema isotérmico:

< v i > : es el promedio de velocidades medias alisadas temporalmente.

Mt: Cantidad de movimiento total

ρ i < v i2 > Si : es la velocidad de cantidad de movimiento debido al movimiento global