Unidad Temática 2 Rev2020

UNIVERSIDAD TECNOLOGICA NACIONAL
FACULTAD REGIONAL ROSARIO

DEPARTAMENTO DE INGENIERIA QUIMICA
NOTAS DE CÁTEDRA
FENOMENOS DE TRANSPORTE
UNIDAD TEMÁTICA 2: BALANCES

ENVOLVENTES Y ECUACIONES
DIFERENCIALES PARA SISTEMAS
ISOTÉRMICOS DE FLUIDOS PUROS
Fenómenos de Transporte - Unidad Temática 2
Revisión Prof. Sandra Godoy: Abril 2020
“Las notas presentadas en esta sección son una recopilación de diversas notas y apuntes de cátedra
dispersas que fueron elaboradas o redactadas a lo largo de muchos años, y han sido readaptadas a la
modalidad actual de enseñanza de la asignatura acorde a los medios disponibles.
Versiones anteriores digitalizadas se pueden encontrar en el sitio
http://www.frro.utn.edu.ar/catedras/catedra.php?m=57
Finalmente, se advierte que estas notas son solo una guía para el abordaje del estudio de los
Fenómenos de Transporte, debiendo consultarse la bibliografía recomendada en cada tema para lograr
un conocimiento pleno de los mismos.”
Cátedra de Fenómenos de Transporte
Página 2 de 26
2.A. ANALISIS ENVOLVENTE EN ESTADO ESTACIONARIO
2.A.1 INTRODUCCION A LA RESOLUCION DE BALANCES

Los balances de cantidad de movimiento, energía y masa pueden formularse aplicando
las leyes de conservación que rigen el comportamiento de un sistema. Estos balances se
pueden enunciar para un determinado período de tiempo (evolución discontinua o cíclica), o
más comúnmente introduciendo el concepto de velocidad de evolución, con las cantidades
relacionadas con la unidad de tiempo.
Su formulación general es la siguiente;
(Velocidad)
(Velocidad) (Velocidad)
(Velocidad) TRANSPORTE
TRANSPORTE GENERACION
ACUMULACION “NETO” de
“NETO” de “NETA” en el
“NETA” en el ENTRADA a
= - SALIDA a través + VOLUMEN
VOLUMEN través
de la del
del de la
SUPERFICIE SISTEMA
SISTEMA SUPERFICIE
del SISTEMA
del SISTEMA
El objetivo es obtener una descripción precisa y conveniente usando expresiones

matemáticas tan rigurosas como sea posible, con un mínimo de parámetros desconocidos.
Estas expresiones obtenidas pueden ser ecuaciones diferenciales ordinarias, en derivadas
parciales, diferenciales finitas, algebraicas, etcétera.
FORMULACIÓN
La secuencia habitual para el planteamiento y resolución de un balance es la siguiente;
1.- Desarrollar un modelo geométrico y elegir un sistema de coordenadas.
2.- Identificar las Entradas y Salidas.
3.- Identificar la Generación (si es positivo se crea y si es negativo se consume).
4.- Verificar el Estado de Régimen (las evoluciones en estado estacionario tienen acumulación
nula).
5.- Definir las Condiciones Límite o de Frontera (condición de una variable en la superficie de
frontera del sistema. Es la constante arbitraria en la solución de una ecuación diferencial.
Existen tantas Condiciones Límite como el orden de la ecuación).
6.- Establecer la Condición Inicial o estado de las variables en el momento inicial (evoluciones
transitorias o no estacionarias).
7.- Definir los Parámetros y los Requisitos. Condiciones o restricciones a cumplir por el sistema,
como por ejemplo, la temperatura, presión, constancia de propiedades, tipo de flujo, etcétera.
Página 3 de 26
EJECUCIÓN
1.- Tomar una envoltura finita lo más sencilla posible acorde con el sistema coordenado
elegido.
2.- Plantear el balance.
3.- Obtener la expresión simbólica representativa.
4.- Resolver la ecuación para lograr conocer la distribución de las variables o “perfiles”.
5.- Simplificarlas para lograr relaciones más simples o “soluciones generales”.
6.- Identificar los valores de frontera, los valores promedio y cualquier otro de interés para la
situación.
2.A.2. ANALISIS ENVOLVENTE DE CANTIDAD DE MOVIMIENTO
Consiste en la aplicación de la ley de conservación de cantidad de movimiento o 2da. Ley

de Newton a un sistema de flujo en estado estacionario, para un fluido de propiedades
constantes que circula en régimen isotérmico.
Velocidad de Velocidad de Velocidad de Velocidad de

ENTRADA de SALIDA de ENTRADA de SALIDA de Sumatoria de
CANTIDAD de CANTIDAD de CANTIDAD de CANTIDAD de las FUERZAS
MOVIMIENTO - MOVIMIENTO + MOVIMIENTO - MOVIMIENTO + EXTERNAS =0
por FLUJO por FLUJO por TRANSP. por TRANSP. (presión y
GLOBAL (ó GLOBAL (ó MOLEC. (ó MOLEC. (ó gravitación)
convectivo) convectivo propagación) propagación)
2.A.2.1 FORMULACIÓN DEL BALANCE

1.- Desarrollar la forma de la envoltura elemental y seleccionar el sistema de coordenadas
adecuado a la misma.
2.- Identificar las Entradas y Salidas. El transporte de cantidad de movimiento a través de la
superficie de la envoltura puede originarse por dos mecanismos;
 Transporte Global o Convectivo, debido al movimiento conjunto del fluido que permite el
ingreso o egreso de cantidad de movimiento al sistema.
 Transporte Molecular o Propagación, a causa de la densidad de flujo de cantidad de
movimiento, o esfuerzos cortantes, originada por la viscosidad del fluido.
3.- Identificar Generación. La cantidad de movimiento puede ser originada, en el elemento de
volumen, por dos tipos de fuerzas;
 Presión, que actúa en la superficie de la envoltura.
 Gravitación, que lo hace en todo el volumen.
4.- Verificar el Estado de Régimen. Para el estado estacionario no existe acumulación de
cantidad de movimiento, es decir, no habrá fuerza resultante actuando sobre el volumen de
Página 4 de 26
fluido, o lo que es lo mismo, no existirá aceleración.

5.- Definir las Condiciones Límite o de Frontera. Para los balances envolventes de cantidad de
movimiento se pueden identificar tres condiciones de frontera;
 Interfase sólido-fluido: la velocidad del fluido es igual a la velocidad con que se mueve la
superficie misma; es decir, se supone que el fluido está adherido a la superficie sólida
con la que se halla en contacto (la ya mencionada condición de no-deslizamiento).
 Interfase líquido-gas: la densidad de flujo de cantidad de movimiento, y por lo tanto el
gradiente de velocidad y el esfuerzo cortante en la fase líquida, es extremadamente
pequeño, y en la mayor parte de los cálculos puede suponerse nulo.
 Interfase líquido-líquido: hay continuidad en ambos lados de la interfase, siendo la
velocidad y la cantidad de movimiento idénticas y común en la interfase para cada
líquido.
6.- Condición Inicial: no existe por ser estado estacionario.
7.- Parámetros y Requisitos: temperatura y propiedades físicas constantes. Trayectoria de las
partículas de fluido rectilíneas y paralelas, régimen de flujo laminar.
2.A.2.2 EJECUCIÓN
1.- Tomar una envoltura finita lo más sencilla posible acorde con el sistema coordenado
elegido.
2.- Plantear el balance según enunciado por notación simbólica.
3.- Tender a cero la envoltura y obtener una ecuación diferencial del esfuerzo en función de la
posición.
4.- Integrar la ecuación para obtener el perfil de densidad de flujo de cantidad de movimiento
(esfuerzo cortante) en función de la posición.
5.- Introducir la ecuación de Newton de la viscosidad y obtener una ecuación diferencial de la
velocidad en función de la posición.
6.- Integrar la ecuación para obtener el perfil de velocidad en función de la posición.
7.- Obtener la velocidad media de flujo, la velocidad máxima, el caudal volumétrico, la velocidad
de flujo de masa, la diferencia de presiones, el esfuerzo cortante en la pared y la fuerza en la
pared ejercida por el fluido.
2.A.2.3 APLICACIONES
En las condiciones enunciadas es válido para flujo en conductos de sección circular,
anular, rectangular, para películas de fluido descendentes en planos o cilindros, canales
abiertos, etcétera.
Página 5 de 26
2.A.2.3.1 Caso De Estudio 1: Balance Envolvente de Cant. de Mov. para una

Película de Fluido Descendente en Plano Inclinado
Como primer ejemplo, consideremos una superficie plana inclinada, tal como se indica en
la siguiente figura:
Estas películas se han estudiado en relación con torres de pared mojada, experiencias
de evaporación y absorción de gases, líneas de baño de chocolate y aplicación de capas de
pintura a rollos de papel. Se supone que la viscosidad y densidad del fluido son constantes y se
considera una región de longitud L, suficientemente alejada de los extremos de la pared, de
forma que las perturbaciones de la entrada y la salida no están incluidas en L; es decir, que en
esta región el componente vz, de la velocidad es independiente de z.
Hipótesis:
• Flujo Laminar, Newtoniano
• Estado Estacionario
• Propiedades físicas constantes (µ,ρ)
Siguiendo los pasos planteados en la sección anterior:
1) Comenzamos aplicando un balance de cantidad de movimiento z sobre una envoltura finita
de espesor ∆x, limitado por los planos z = 0 y z = L, y que se extiende hasta una distancia W
en la dirección y.
Página 6 de 26
τ xz x
τ xz x +∆x
2) Plantear el balance según enunciado por notación simbólica
TRANSPORTE MOLECULAR FLUJO GLOBAL FUERZAS EXTERNAS

Como vz vale lo mismo para z=0 y z=L
3) Tender a cero la envoltura y obtener una ecuación diferencial del esfuerzo en función de la
posición.
Página 7 de 26
Dividiendo miembro a miembro por W L Δx, y aplicando límite
4) Integrar la ecuación para obtener el perfil de densidad de flujo de cantidad de movimiento

(esfuerzo cortante) en función de la posición.
Para determinar el valor de C1 podemos definir la condición de contorno en la interfase Líquido-

Gas (x=0)
Para , C1=0
5) Introducir la ecuación de Newton de la viscosidad y obtener una ecuación diferencial de la
velocidad en función de la posición.
6) Integrar la ecuación para obtener el perfil de velocidad en función de la posición.
Para determinar el valor de C2 podemos definir la condición de contorno en la interfase Sólido-

Líquido (x=δ -espesor de la película de fluido-)
Para
7) Obtener la velocidad máxima y la velocidad media de flujo, el caudal volumétrico, y la fuerza

en la pared ejercida por el fluido.
Página 8 de 26
• Velocidad Máxima
La velocidad máxima se produce cuando x = 0
• Velocidad Media
• Velocidad volumétrica de flujo o caudal volumétrico
• Fuerza en la pared ejercida por el fluido
Comp. Z de TODO el peso del fluido contenido en la película
8) Resultados finales. Representación gráfica de las distribuciones de velocidad y de densidad

de flujo de cantidad de movimiento.
Página 9 de 26
Distribución de velocidades
Distribución de densidad de flujo

de cantidad de movimiento
Válido para N° Reynolds <[4-25]
2.A.2.3.2 Caso De Estudio 2: Balance Envolvente de Cant. de Mov. para el Flujo

de un Fluido a Través de un Cilindro
El flujo de fluidos en tubos circulares se encuentra con frecuencia en física, química,
biología e ingeniería. El flujo laminar de fluidos en tubos circulares puede analizarse mediante
el balance de cantidad de movimiento que se ha descrito en el caso anterior, con la diferencia
de que se introduce aquí el uso de coordenadas cilíndricas, que son las coordenadas naturales
para describir las posiciones en una tubería circular.
Hipótesis:
• Flujo Laminar, Newtoniano
• Estado Estacionario
• Propiedades físicas constantes (µ,ρ)
• Tubo «muy largo» de longitud L y radio R (no existen “efectos finales”).
Siguiendo los pasos planteados en la sección anterior:

1) Tomamos una envoltura cilíndrica de espesor Δr y longitud L
2) Comenzamos a plantear el balance de cantidad de movimiento en la dirección z
considerando las distintas contribuciones según notación simbólica.
Página 10 de 26
TRANSPORTE MOLECULAR FLUJO GLOBAL FUERZAS EXTERNAS
3) Tender a cero la envoltura y obtener una ecuación diferencial del esfuerzo en función de
la posición.
Dividiendo miembro a miembro por 2 π L Δr, y aplicando límite
Siendo:
4) Integrar la ecuación para obtener el perfil de densidad de flujo de cantidad de
movimiento (esfuerzo cortante) en función de la posición.
Página 11 de 26
Para determinar el valor de C1 podemos definir la condición de contorno en el centro del

tubo (r=0)
Para C1=0
5) Introducir la ecuación de Newton de la viscosidad y obtener una ecuación diferencial de

la velocidad en función de la posición.
dvz
τ rz = − µ
dr
6) Integrar la ecuación para obtener el perfil de velocidad en función de la posición.
Para determinar el valor de C2 podemos definir la condición de contorno en la interfase

Sólido-Líquido (r=R –radio del cilindro-)
Para =
7) Obtener la velocidad media de flujo, la velocidad máxima, el caudal volumétrico, la

velocidad de flujo de masa, el esfuerzo cortante en la pared y la fuerza en la pared
ejercida por el fluido.
• Velocidad Máxima
La velocidad máxima se produce cuando r = 0
• Velocidad Media
rdrdϕ
rdrdϕ
Página 12 de 26
• Velocidad volumétrica de flujo (Producto área x V media)
Ecuación de Hagen-Poiseuille
• Fuerza en la pared ejercida por el fluido
8) Resultados finales. Representación gráfica de las distribuciones de velocidad y de

densidad de flujo de cantidad de movimiento.
Estos resultados son válidos solamente para valores del número de Reynolds inferiores a
2100, para los que el flujo es laminar. Para este sistema se acostumbra definir el número de
Reynolds por Re = D<vz>ρ/µ, siendo D = 2R, el diámetro del tubo.
2.A.2.3.3 Caso De Estudio 3: Balance Envolvente de Cant. de Mov. para el Flujo

de un Fluido a Través de una Sección de Corona Circular
Ver este caso en sección 2.4, del libro: Fenómenos de Transporte, R.B.Bird, W.E,
Stewart, E.N. Lightfoot. Ed. Reverté.
Página 13 de 26
2.B. ECUACIONES DIFERENCIALES PARA FLUJO DE FLUIDOS ISOTERMICOS
2.B.1 INTRODUCCION
Al trabajar con un nuevo problema de flujo no siempre es necesario formular un balance
envolvente de cantidad de movimiento. Es más rápido, más fácil, y más seguro, partir de las
ecuaciones de conservación de la materia y la cantidad de movimiento, expresadas en la forma
general, y simplificarlas con el fin de adaptarlas al problema de que se trate. Estas dos
ecuaciones describen todos los problemas de flujo viscoso isotérmico de un fluido puro, y se
denominan a veces «ecuaciones de variación».
La «ecuación de continuidad», corresponde a la aplicación de la ley de la conservación
de la materia a un pequeño elemento de volumen situado en el seno de un fluido en
movimiento. La «ecuación de movimiento», es una generalización del balance de cantidad de
movimiento.
También se utiliza la ecuación de movimiento para deducir una expresión que describe la
interconversión de las distintas formas de la energía mecánica de un fluido en movimiento. Esta
ecuación es particularmente útil para describir la degradación de la energía mecánica en
energía calorífica, que acompaña a todos los procesos reales de flujo. También es la base del
importante balance macroscópico de energía mecánica, o ecuación de Bernoulli, que se
estudiará mas adelante.
Las ecuaciones de variación también se pueden presentar en función de variables
adimensionales. Al escribir en esta forma las ecuaciones, se reúnen los «factores de escala»
(es decir, tamaño, del sistema, velocidad media del fluido, y propiedades del mismo) en un
corto número de relaciones adimensionales que resultan útiles para caracterizar los sistemas
de flujo, de modo de utilizar estas relaciones adimensionales para efectuar estudios empíricos
con un modelo a escala reducida, que permite un análisis de los sistemas demasiado
complejos.
2.B.2 DERIVADAS DE FUNCIONES

Para interpretar el significado físico de los tres tipos de derivadas con respecto al tiempo
que se utilizan en el curso, utilizaremos como ejemplo la concentración de peces en un río.
Puesto que los peces se están moviendo, su concentración C será una función de la
posición (x, y, z) y del tiempo (t).
Así si estuviéramos parados sobre un puente y observáramos cómo varía la
concentración de peces exactamente debajo de nosotros con el tiempo, lo que estaríamos
registrando es cómo varía la concentración con el tiempo, para una posición fija en el espacio.
Esto es la “derivada parcial de c con respecto a t: dc/dθ, manteniendo constantes x, y, z”.
Página 14 de 26
Si en lugar de estar sobre el puente, vamos en una lancha a motor que se mueve en el
río en todas direcciones, unas veces en contra de la corriente, otras a través, y tal vez otras a
favor de la corriente. Al referir la variación de la concentración de peces con respecto al tiempo,
los números que resultan han de reflejar también el movimiento de la lancha. Esto queda
representado por la “derivada total con respecto al tiempo”, la que viene dada por:
dc dc dc dx dc dy dc dz
= + + +
dθ dθ dx dθ dy dθ dz dθ
en la que dx/dt, dy/dt y dz/dt son los componentes de la velocidad de la lancha.

Finalmente, supongamos que vamos en una canoa a la que no se comunica energía, sino que
simplemente flota. En este caso, la velocidad del observador es exactamente la misma que la
velocidad de la corriente v. Al referir la variación de la concentración de peces con respecto al
tiempo, los números dependen de la velocidad local de la corriente. Esta derivada es una clase
especial de derivada total con respecto al tiempo que se denomina “derivada substancial” o,
a veces (más lógicamente), «derivada siguiendo el movimiento». Está relacionada con la
derivada parcial con respecto al tiempo de la forma siguiente:
Dc dc dc dc dc
= + vx + vy + vz
Dθ dθ dx dy dz
en la que vx, vy, y vz, son los componentes de la velocidad local del fluido v.
2.B.3 LA ECUACION DE CONTINUIDAD

Consideramos un elemento estacionario de volumen, un cubo de lados Δx, Δy, Δz, a
través del cual circula un fluido cuya densidad y velocidad son funciones de la posición y el
tiempo.
Si aplicamos un balance de materia al mismo, se podrá calcular el flujo de masa por
unidad de tiempo que atraviesa cada cara de ese elemento de volumen, de acuerdo a la
siguiente expresión;
Velocidad de Velocidad de Velocidad de

Acumulación de = Entrada de - Salida de
Materia Materia Materia
Página 15 de 26
La velocidad de acumulación de materia será la rapidez de variación, respecto del

tiempo, de la masa en el interior del cubo, o sea;
∆ρ
( ∆x∆y∆z )
∆θ
Para cada par de caras paralelas, la velocidad de entrada y salida de materia podrá
expresarse;
 caras perpendiculares al eje x: ( ρ v  − ρ v  ) ∆z∆y

x x x x +∆x
 caras perpendiculares al eje y: ( ρ v  − ρ v  ) ∆x∆z

y y y y +∆y
 caras perpendiculares al eje z: ( ρ v  − ρ v  ) ∆x∆y

z z z z +∆z
El balance de materia quedará entonces;
∆ρ
∆θ
( ∆x∆y∆=
z) ( ρ v  − ρ v  ) ∆z∆y + ( ρ v  − ρ v  ) ∆x∆z + ( ρ v  − ρ v  ) ∆x∆y
x x x x +∆x y y y y +∆y z z z z +∆z
Dividiendo la expresión por el elemento de volumen (Δx Δy Δz), y aplicando el límite

cuando el volumen y el tiempo tienden a cero, se obtiene;
∂ρ  ∂ρ vx ∂ρ v y ∂ρ vz 
=
− + + 
∂θ  ∂x ∂y ∂z 
Página 16 de 26
Esta expresión es la Ecuación de Continuidad, que describe la variación de la densidad

para un elemento de volumen fijo en el espacio, como resultado de las variaciones del vector
velocidad másica (ρ v).
Utilizando notación vectorial-tensorial, también puede escribirse;
∂ρ
+ ( ∇.ρ v ) =0
∂θ
El vector velocidad másica corresponde a una densidad de flujo de materia, y también

significa la velocidad de concentración de cantidad de movimiento. El término ( ∇.ρ v ) se
conoce como divergencia de (ρv), y significa la velocidad neta con que disminuye la densidad
de flujo de materia por unidad de volumen. La ecuación establece que la velocidad de aumento
de densidad en el interior del elemento de volumen fijo en el espacio, es igual a la velocidad
“neta” de flujo de materia dividido por su volumen.
Efectuando la diferenciación de los términos y reuniendo las derivadas de la densidad en
el primer miembro, se obtiene;
∂ρ  ∂ρ ∂ρ ∂ρ   ∂v ∂v y ∂vz 
+  vx + vy + vz =−ρ  x + + 
∂θ  ∂x ∂y ∂z   ∂x ∂y ∂z 
Se observará que el primer miembro es la derivada substancial de la densidad, o sea, la

derivada con respecto al tiempo para un recorrido que sigue el movimiento del fluído. Por lo
tanto, la expresión quedará;
Dρ
+ ρ ( ∇.v ) =
0
Dθ
La anterior representa la variación de la densidad de una partícula de fluido tal como la

observa quien se mueve con el fluido.
2.B.3.1 Formas especiales importantes de la ecuación de continuidad

 Para estado estacionario, la densidad no varía en el tiempo en el punto fijo del
espacio, y la expresión quedará;
( ∇.ρ v ) =
0
 Para fluidos incompresibles (densidad constante) se debe cumplir que;
Página 17 de 26
( ∇.v ) =
0
Aunque ningún fluido es totalmente incompresible, en general se admite que los líquidos
tienen densidad constante, sin introducir casi error. Aún para los gases, dentro de moderados
cambios de presión pueden suponerse con densidad constante.
Para que la última sea válida, sólo es necesario que la densidad sea constante para el
elemento de fluido que se mueve a lo largo de una “línea de corriente”, entendiendo por tal a
la curva imaginaria que conecta sucesivos puntos del espacio cuyos vectores velocidad son
tangentes a la misma, indicando la dirección del movimiento y la trayectoria del elemento
Dρ
( =0).
Dθ
2.B.4 LA ECUACION DEL MOVIMIENTO

Para un fluido viscoso, en flujo isotérmico, se puede aplicar el principio de conservación
de la cantidad de movimiento al mismo elemento de volumen utilizado anteriormente.
El balance de cantidad de movimiento, aplicado al elemento de lados Δx, Δy, Δz se
puede escribir;
Velocidad de Velocidad de Velocidad de Suma de
acumulación de entrada de _ salida de fuerzas que
= +
cantidad de cantidad de cantidad de actúan sobre el
movimiento movimiento movimiento sistema
Este balance tiene en cuenta que;

 La expresión representa la ecuación de un vector, con componentes para cada una de
las direcciones coordenadas (x, y, z).
 El fluido es puro e isotérmico.
 El comportamiento puede ser no-estacionario.
 El movimiento del fluido es arbitrario a través de la superficie del elemento de
volumen.
 La velocidad de entrada y salida de cantidad de movimiento se puede efectuar por dos
mecanismos; (a) Flujo global o convectivo, o sea, la contribución del movimiento
conjunto de las partículas del fluido que atraviesan una superficie, y (b) Transporte
molecular o propagación, debido a la contribución viscosa originada en los gradientes
de velocidad.
 La aceleración gravitatoria es un vector de tres componentes.
Para iniciar el balance solo tendremos en cuenta la componente x de cada uno de los
Página 18 de 26
términos de la ecuación, tal como lo muestra la figura siguiente. Los correspondientes a los
ejes y, z se pueden obtener por analogía.
1- El termino de flujo global o velocidad de flujo de cantidad de movimiento por

convección a través de las seis caras, puede representarse por el producto de la componente
de la densidad de flujo de masa (ρv) con la velocidad en dirección x (vx) multiplicadas por el
área transversal o superficie de la cara en cuestión. Por lo tanto, se puede escribir;
 Entrada en la cara x: ( ρ vx )vx ∆y∆z x
 Salida por x + ∆ x: ( ρ vx )vx ∆y∆z x + ∆x
 Entrada en la cara y: ( ρ v y )vx ∆x∆z y
 Salida por y+ ∆ y: ( ρ v y )vx ∆x∆z y +∆y
 Entrada en la cara z: ( ρ vz )vx ∆x∆y z
 Salida por z+ ∆ z: ( ρ vz )vx ∆x∆y z +∆z
2- Para la propagación, tenemos que puede formularse para todas las caras en forma
similar, donde como ya se vio en la Unidad 1, las densidades de flujo de cantidad de
movimiento pueden considerarse como esfuerzos (τij), donde el primer subíndice indica la
dirección de la propagación y el segundo la de la velocidad. Para este caso, el segundo
Página 19 de 26
subíndice será siempre ”x”. El término τxx es el esfuerzo normal que actúa sobre la cara x, y los
términos τyx y τzx son los esfuerzos tangenciales, quedando:
 Entrada en la cara x: τ xx ∆y ∆z x
 Salida por x+ ∆ x: τ xx ∆y ∆z x + ∆x
 Entrada en la cara y: τ yx ∆x∆z y

 Salida por y+ ∆ y: τ yx ∆x∆z y +∆y
 Entrada en la cara z: τ zx ∆x∆y z
 Salida por z+ ∆ z: τ zx ∆x∆y z +∆z
3- Las únicas fuerzas que se tienen en cuenta en este caso son las de presión del fluido y
la gravitación, en la dirección del moviendo considerado (dirección x), resultando;
 Para la presión en x: ( p x − p x +∆x )∆y∆z
 Para la gravitación en x: ρ g x (∆x∆y∆z )
4- la velocidad de acumulación de cantidad de movimiento (en la dirección “x”) en el

elemento de volumen será;
∆ρ vx
 (∆x∆y∆z )
velocidad de acumulación en x:
∆θ
Sustituyendo cada una de las anteriores en la expresión general correspondiente al
balance y dividimos por el elemento de volumen (Δx Δy Δz) obteniendo una expresión
incremental de la variación de cantidad de movimiento en el tiempo y espacio.
∆ρ vx ρ vx vx x − ρ vx vx x +∆x ρ v y vx y − ρ v y vx y +∆y
= +
∆θ ∆x ∆y
Página 20 de 26
ρ vz vx z − ρ vz vx z +∆z τ xx x − τ xx x +∆x τ yx x −τ yx y +∆y

+ + + +
∆z ∆x ∆y
τ zx x −τ zx z +∆z p x − p x +∆x
+ + + ρ gx
∆z ∆x
Aplicando límite cuando la norma de los incrementos (∆x, ∆y, ∆z , ∆θ ) tiende a cero, se
obtiene la componente “x” de la ecuación del movimiento.
 Componente “x”:
∂ρ vx  ∂ ∂ ∂   ∂ ∂ ∂  ∂p
− ρ vx vx + ρ v y vx + ρ vz vx  −  τ xx + τ yx + τ zx  − + ρ g x
=
∂θ  ∂x ∂y ∂z   ∂x ∂y ∂z  ∂x
De igual manera se pueden obtener las correspondientes a las coordenadas “y”, “z”:
 Componente “y”:
∂ρ v y ∂ ∂ ∂   ∂ ∂ ∂  ∂p
−  ρ vx v y +
= ρ v y v y + ρ vz v y  −  τ xy + τ yy + τ zy  − + ρgy
∂θ  ∂x ∂y ∂z   ∂x ∂y ∂z  ∂y
 Componente “z”:
∂ρ vz  ∂ ∂ ∂   ∂ ∂ ∂  ∂p
−  ρ vx vz +
= ρ v y vz + ρ vz vz  −  τ xz + τ yz + τ zz  − + ρ gz
∂θ  ∂x ∂y ∂z   ∂x ∂y ∂z  ∂z
Los términos se pueden identificar de la siguiente forma;

 Velocidad másica (ρv), un vector cuyos tres componentes espaciales son: ρvx, ρvy, ρvz
 Aceleración gravitatoria (g), un vector cuyos tres componentes son: gx, gy, gz.
∂p ∂p ∂p
 Gradiente de presión (∇p ) , un vector cuyos tres componentes son: , ,
∂x ∂y ∂z
 Tensor esfuerzo τij de segundo orden, correspondiente a la “propagación”, con nueve
componentes: τxx, τxy, τxz, τyx, τyy, τyz, τzx, τzy,, τzz
 Densidad de flujo convectivo de cantidad de movimiento (ρvv), producto diádico 1 de los
vectores velocidad por velocidad másica, de nueve componentes: ρvxvx, ρvxvy, ρvxvz,
ρvyvx, ρvyvy, ρvyvz, ρvzvx, ρvzvy, ρvzvz
En notación vectorial-tensorial, se puede escribir:
Página 21 de 26
∂ρ v
= − [∇.ρ vv ] − [∇.τ ] − ∇p + ρ g
∂θ
El término [∇.ρ vv ] representa la velocidad de perdida de cantidad de movimiento (un
vector) por unidad de volumen debido al flujo de fluido, mientras que [∇.τ ] representa la
perdida de cantidad de movimiento debido a la viscosidad. Estos términos no son divergencias

simples, por el origen tensorial del producto diádico y tensor esfuerzo.
Por definición de Derivada substancial:
Dvx ∂p  ∂ ∂ ∂ 
ρ − −  τ xx + τ yx + τ zx  + ρ g x
=
Dθ ∂x  ∂x ∂y ∂z 
Idéntica transformación se puede realizar para las direcciones “y” y “z”, obteniendo una
expresión en notación compacta que describe el movimiento de una partícula de fluido que se
mueve bajo el efecto de las fuerzas que actúan sobre ella, a saber;
Dv
ρ = −∇p − [ ∇.τ ] + ρ g
Dθ
La ecuación del movimiento expresada así establece una expresión de la segunda ley de
Newton, en la cual (masa x aceleración= suma de fuerzas). El balance de cantidad de
movimiento aplicado a un fluido es totalmente equivalente a la segunda ley de Newton, donde
las fuerzas que actúan son las originadas en la presión, la viscosidad, y la gravitación.
Los términos se pueden identificar de la siguiente forma:
∂ρ v
 : Velocidad de aumento de cantidad de movimiento por unidad de volumen.
∂θ
 − [∇.ρ vv ] : Velocidad de ganancia de cantidad de movimiento por convección
por unidad de volumen.
 −∇p : Fuerza de presión sobre el elemento por unidad de volumen.
 − [∇.τ ] : Velocidad de ganancia de cantidad de movimiento por transporte
viscoso, por unidad de volumen (equivale a fuerza viscosa sobre elemento por
unidad de volumen).
 ρ g : Fuerza de gravitación que actúa sobre el elemento de volumen.
1
Producto diádico (o diádas) son los tensores que resultan de multiplicar entre sí dos vectores
Página 22 de 26
∂v
 ρ : Masa por unidad de volumen, por aceleración.
∂θ
En cada caso, la expresión con la “derivada parcial” representa un balance aplicado a un
elemento de volumen fijo en el espacio, mientras que la correspondiente a la “derivada
sustancial” describe las variaciones del elemento que sigue el movimiento del fluido. Ambas
son validas para medios continuos, sin importar la naturaleza de la relación esfuerzo-
deformación, es decir, pueden aplicarse tanto a fluidos newtonianos como no-newtonianos.
Aplicando las relaciones correspondientes entre esfuerzo-deformaciones, aquellas
expresiones pueden transformarse para obtener las distribuciones de velocidad del fluido,
obteniendo ecuaciones generales que describen el movimiento del fluido.
2.B.5 ECUACIONES DE NAVIER-STOKES

Hemos obtenido anteriormente, la Ecuación de Continuidad,
∂ρ
= −(∇.ρ v)
∂θ
Y la Ecuación del Movimiento en función de los esfuerzos cortantes,
Dv
ρ = −∇p − [∇.τ ] + ρ g
Dθ
En esta última se establece que un pequeño elemento de volumen que se mueve con el
fluido es acelerado por las fuerzas que actúan sobre él (fuerzas de presión, gravitación y
viscosidad). Estas ecuaciones son la representación de la segunda Ley de Newton a un fluido
en movimiento.
Sin embargo, se puede demostrar que tenemos seis componentes del esfuerzo
desconocidas y tres componentes de velocidad también desconocidas, con lo que existen
nueve incógnitas en las ecuaciones.
Se tratará de reducir el exceso de incógnitas por la introducción de ecuaciones
constitutivas que relacionan el esfuerzo con la velocidad, para obtener expresiones que
relacionen solamente velocidades, presiones y gravitación.
Para ello se desarrollara la ecuación cartesiana de cantidad de movimiento para una
dirección (el eje x), que puede expresarse:
∂ρ vx  ∂ρ vx vx ∂ ρ v y vx   ∂τ xx ∂τ yx ∂τ zx  ∂p
=
− + + ∂ρ v vz x
− + + − + ρ gx
∂θ  ∂x ∂y ∂z   ∂x ∂y ∂z  ∂x
Página 23 de 26
Que, por definición de derivada substancial permite obtener la componente x de la

aceleración por unidad de volumen de fluido:
Dvx ∂p ∂τ xx ∂τ yx ∂τ zx
ρ =
− −( + + ) + ρ gx
Dθ ∂x ∂x ∂y ∂z
Introduciendo las relaciones para fluidos Newtonianos, entre los distintos esfuerzos en
función de los gradientes de velocidad y las propiedades del fluido (se puede encontrar la
demostración en la literatura), obtendremos las ecuaciones generales de movimiento para un
fluido newtoniano que presenta variación de la densidad y la viscosidad
Dvx ∂p ∂  ∂vx 2  ∂  ∂vx ∂vy  ∂  ∂vx ∂vz 

ρ =
− +  2µ − µ (∇.v)  +  µ ( + + µ( + + ρ gx
Dθ ∂x ∂x  ∂x 3  ∂y  ∂y ∂x  ∂z  ∂z ∂x 
Ecuaciones similares se pueden obtener para las direcciones “y”, “z”.

Estas ecuaciones, junto a la ecuación de continuidad, la ecuación de estado p=p(ρ), la
variación de viscosidad con la densidad µ=µ(ρ), y las condiciones iniciales y límite, determinan
en forma completa la presión, densidad y los componentes de velocidad de un fluido en
movimiento isotérmico.
Para el caso de fluidos newtonianos (viscosidad constante) e incompresibles (densidad
constante), la ecuación de continuidad predice que,
(∇.v) = 0
Por lo que la anterior se puede simplificar para obtener la conocida Ecuación de Navier-
Stokes, que en coordenadas cartesianas y en notación vectorial-tensorial viene representada
por:
Dv
ρ = -р
∇ + µ∇ 2 v + ρg
Dθ
Las componentes cartesianas, así como las expresiones correspondientes a
coordenadas cilíndricas y esféricas se encuentran tabuladas.
En principio, con la Ecuación del Movimiento y la Ecuación de Continuidad se puede
obtener la solución de cualquier problema de flujo isotérmico.
2.B.6 ECUACIÓN DE EULER

Para el caso de considerar un fluido no-viscoso, el término de esfuerzos constantes es
nulo, y la Ecuación del Movimiento se reduce a:
Dv
ρ = -р
∇ + ρg
Dθ
Página 24 de 26
Que dice que los cambios de velocidad del fluido se deben exclusivamente a las fuerzas
de presión y gravitatorias. Si bien los fluidos reales tienen viscosidad, esta ecuación se puede
aplicar el caso de que los efectos viscosos son poco importantes.
2.B.7 ECUACIÓN DE LA ENERGÍA MECÁNICA

Dada la Ecuación de Movimiento, expresada de la forma:
Dv
ρ -р - [τ]
=∇ ∇.ρ + g
Dθ
Se forma el producto escalar de la velocidad con la misma. La ecuación resultante es una
ecuación escalar, y describe la velocidad de variación de la energía cinética por unidad de
masa (½ v2) para el elemento de fluido que se mueve con la corriente, quedando:
D (½v 2 )
ρ = -(р)-
v.∇ (v[τ]
∇. ) + ρ(v.g )
Dθ
La Ecuación se conoce con el nombre de Ecuación de la Energía Mecánica y representa
la interconversión de energía mecánica de un fluido en movimiento.
Desarrollando la ecuación en función de las derivadas parciales, se puede escribir:
 Para la derivada sustancial, según la ecuación de continuidad,
D (½v 2 ) ∂
ρ = (½ρv 2 ) + (∇½ρv 2 v )
Dθ ∂θ
 Para el término de presión, según las reglas de derivación.
(∇.рv=
) (v.∇р) + р(∇v )
 Para la contribución viscosa, de igual forma,
(∇.[τ.v ]) = (v.[∇.τ]) + ( τ : ∇v)

Quedando entonces la Ecuación de la Energía Mecánica en función de un elemento
estacionario de volumen por donde circula el fluido,
∂
(½ρv 2τ)= ]-(∇
р½τ ρ-р
v 2 v-) - (∇.- v ) - v-) (∇.[ .v ) ( : ∇v) + ρ(v.g )
( ∇.
∂θ
Donde el significado físico de cada término es el siguiente:
∂
 (½ρv 2 ) : Velocidad de incremento de energía cinética por unida de volumen.
∂θ
 - (∇½ρv 2 v) : Velocidad neta de entrada de energía cinética debida al fluido global.
Página 25 de 26
 - (∇.рv ) : Velocidad de trabajo producido por la presión de los alrededores sobre el

elemento de volumen.
 - p (-∇.v) : Velocidad de conversión “reversible” en energía interna.
 - (∇.[τ .v]) : Velocidad de trabajo producido por las fuerzas viscosas que actúan sobre el
 - ( −τ : ∇v) : Velocidad de conversión “irreversible” en energía interna (1)
 ρ(v.g ) : Velocidad de trabajo producido por la fuerza de gravedad que actúa sobre el
El término (1) será desarrollado posteriormente, pero indica que en los sistemas de flujo
hay una degradación de energía mecánica en energía calorífica, y que los procesos reales no
son reversibles.
Un sistema “isotérmico”, es por lo tanto, aquel en que el calor generado o absorbido no
da variación apreciable de temperatura.
Esta ecuación servirá posteriormente para la deducción del balance de energía mecánica
macroscópico o Ecuación de Bernoulli.
Página 26 de 26

Unidad Temática 2 Rev2020

Cargado por

Copyright:

Formatos disponibles

Unidad Temática 2 Rev2020

Cargado por

Información del documento

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 Temática 2 Rev2020

Cargado por

Copyright:

Formatos disponibles

UNIVERSIDAD TECNOLOGICA NACIONAL

FACULTAD REGIONAL ROSARIO

UNIDAD TEMÁTICA 2: BALANCES

Cátedra de Fenómenos de Transporte

2.A. ANALISIS ENVOLVENTE EN ESTADO ESTACIONARIO

2.A.1 INTRODUCCION A LA RESOLUCION DE BALANCES

El objetivo es obtener una descripción precisa y conveniente usando expresiones

2.A.2. ANALISIS ENVOLVENTE DE CANTIDAD DE MOVIMIENTO

Consiste en la aplicación de la ley de conservación de cantidad de movimiento o 2da. Ley

Velocidad de Velocidad de Velocidad de Velocidad de

2.A.2.1 FORMULACIÓN DEL BALANCE

fluido, o lo que es lo mismo, no existirá aceleración.

2.A.2.3.1 Caso De Estudio 1: Balance Envolvente de Cant. de Mov. para una

2) Plantear el balance según enunciado por notación simbólica

TRANSPORTE MOLECULAR FLUJO GLOBAL FUERZAS EXTERNAS

Dividiendo miembro a miembro por W L Δx, y aplicando límite

4) Integrar la ecuación para obtener el perfil de densidad de flujo de cantidad de movimiento

Para determinar el valor de C1 podemos definir la condición de contorno en la interfase Líquido-

6) Integrar la ecuación para obtener el perfil de velocidad en función de la posición.

Para determinar el valor de C2 podemos definir la condición de contorno en la interfase Sólido-

7) Obtener la velocidad máxima y la velocidad media de flujo, el caudal volumétrico, y la fuerza

• Velocidad volumétrica de flujo o caudal volumétrico

• Fuerza en la pared ejercida por el fluido

Comp. Z de TODO el peso del fluido contenido en la película

8) Resultados finales. Representación gráfica de las distribuciones de velocidad y de densidad

Distribución de densidad de flujo

Válido para N° Reynolds <[4-25]

2.A.2.3.2 Caso De Estudio 2: Balance Envolvente de Cant. de Mov. para el Flujo

Siguiendo los pasos planteados en la sección anterior:

TRANSPORTE MOLECULAR FLUJO GLOBAL FUERZAS EXTERNAS

Dividiendo miembro a miembro por 2 π L Δr, y aplicando límite

Para determinar el valor de C1 podemos definir la condición de contorno en el centro del

5) Introducir la ecuación de Newton de la viscosidad y obtener una ecuación diferencial de

6) Integrar la ecuación para obtener el perfil de velocidad en función de la posición.

Para determinar el valor de C2 podemos definir la condición de contorno en la interfase

7) Obtener la velocidad media de flujo, la velocidad máxima, el caudal volumétrico, la

• Velocidad volumétrica de flujo (Producto área x V media)

• Fuerza en la pared ejercida por el fluido

8) Resultados finales. Representación gráfica de las distribuciones de velocidad y de

2.A.2.3.3 Caso De Estudio 3: Balance Envolvente de Cant. de Mov. para el Flujo

2.B. ECUACIONES DIFERENCIALES PARA FLUJO DE FLUIDOS ISOTERMICOS

2.B.2 DERIVADAS DE FUNCIONES

en la que dx/dt, dy/dt y dz/dt son los componentes de la velocidad de la lancha.

2.B.3 LA ECUACION DE CONTINUIDAD

Velocidad de Velocidad de Velocidad de

La velocidad de acumulación de materia será la rapidez de variación, respecto del

 caras perpendiculares al eje x: ( ρ v  − ρ v  ) ∆z∆y

 caras perpendiculares al eje y: ( ρ v  − ρ v  ) ∆x∆z

 caras perpendiculares al eje z: ( ρ v  − ρ v  ) ∆x∆y

El balance de materia quedará entonces;

Dividiendo la expresión por el elemento de volumen (Δx Δy Δz), y aplicando el límite

Esta expresión es la Ecuación de Continuidad, que describe la variación de la densidad

El vector velocidad másica corresponde a una densidad de flujo de materia, y también

Se observará que el primer miembro es la derivada substancial de la densidad, o sea, la

La anterior representa la variación de la densidad de una partícula de fluido tal como la

2.B.3.1 Formas especiales importantes de la ecuación de continuidad

2.B.4 LA ECUACION DEL MOVIMIENTO

Este balance tiene en cuenta que;

1- El termino de flujo global o velocidad de flujo de cantidad de movimiento por

 Entrada en la cara x: ( ρ vx )vx ∆y∆z x