UNIVERSIDAD GABRIEL RENE MORENO

FACULTAD DE CIENCIAS EXACTAS Y TECNOLOGIA.

DEPARTAMENTO DE INGENIERIA DE ALIMENTOS

FENOMENOS DE TRANSPORTE

NOTAS DE CÁTEDRA: “UNIDAD TEMÁTICA 2”

ƒ ANALISIS ENVOLVENTE EN ESTADO ESTACIONARIO
ƒ EC. DIFERENCIALES PARA FLUJO DE FLUIDOS ISOTERMICOS
ƒ ANALISIS DIMENSIONAL

2021
Fenómenos de Transporte – Unidad Temática 2
2021

Página 2 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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.

FORMULACION
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 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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

APLICACIONES
Campo de utilización, estudio de limitaciones, ejemplos para diseño.

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)

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

Página 4 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

ƒ 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
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 nula con respecto al sólido en la
interfase, de acuerdo a la ya mencionada condición de no-deslizamiento.
ƒ Interfase líquido-gas: se considera que no hay resistencia en la fase líquida por parte de
la fase gaseosa, por lo que el gradiente de velocidad y el esfuerzo cortante en el líquido
es nulo.
ƒ Interfase líquido-líquido: para líquidos inmiscibles de distinta densidad, hay continuidad
en ambos lados de la interfase, siendo la velocidad idéntica 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.

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

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 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

2.B. ECUACIONES DIFERENCIALES PARA FLUJO DE FLUIDOS ISOTERMICOS

INTRODUCCION
El concepto de medio continuo explica el comportamiento de la materia como si se
encontrara llenando por completo el espacio que ocupa, no tomando en cuenta el
comportamiento individual de las moléculas. Al adoptar este punto de vista, y dentro de las
limitaciones para el que es válido, se puede tener una descripción matemática con funciones
continuas, para los intervalos considerados de las coordenadas de espacio y tiempo. Las
cantidades físicas resultantes son independientes de cualquier sistema particular de
coordenadas que las describa, aunque pueden referirse a un sistema apropiado por
conveniencia, tales como cartesianas, cilíndricas, etcétera.

TENSORES
Son entidades matemáticas que tienen existencia independiente de cualquier sistema de
coordenadas, aun-que pueden ser especificados mediante componentes. Las leyes físicas se
expresan entonces por ecuaciones tensoriales, que son válidas en cualquier sistema
coordenado, es decir, que son invariantes. Los tensores pueden ser clasificados por su orden,
reflejando también el número de componentes que poseen en un espacio n-dimensional.
ƒ Tensores de Orden Cero o Escalares; se especifican por una componente en el
espacio tridimensional, o sea, se caracterizan por su magnitud o valor, ej.; masa,
temperatura, volumen, tiempo, densidad.
ƒ Tensores de Orden Uno o Vectores; tienen tres componentes coordenadas en el
espacio físico, se caracterizan por ser magnitudes orientadas (en el espacio de tres
dimensiones tienen módulo, dirección y sentido), ej.; velocidad, cantidad de movimiento,
fuerza, densidad de flujo de energía.
ƒ Tensores de Orden Dos o Diadas; son el producto indeterminado de dos vectores, se
caracterizan por ser magnitudes de nueve componentes en el espacio tridimensional,
ej.; densidad de flujo de cantidad de movimiento, esfuerzos, producto diádico de
velocidades.

CAMPO TENSORIAL
Asocia un determinado tensor (función de sus coordenadas espaciales y el tiempo) a un
vector posición que varía en una región particular del espacio en un intervalo particular del
tiempo. Este campo será continuo (o diferenciable) si las componentes del tensor son
funciones continuas (o diferenciables). Si las componentes son funciones de sus coordenadas
e independientes del tiempo, se llama estacionario.
ƒ Campo Escalar: Cuando en cada punto de la región la magnitud es un escalar, por ej.;

Página 6 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

la densidad de un fluido. Se puede representar por;

ρ = ρ ( xi ,θ )
ƒ Campo Vectorial: Cuando en cada punto de la región existe un vector, por ej.; la
velocidad de las partículas de un fluido. Se representa por;

vi = vi ( x ,θ )
ƒ Campo Tensorial de 2do. Orden: Cuando en cada punto de la región existe una diada,
por ej.; el tensor deformación. Representado por;

τ ij = τ ij ( x , θ )
Si observamos una región del espacio llena de un fluido, que se encuentra en movimiento
respecto de una terna derecha de ejes cartesianos estacionarios, podemos decir que estamos
en presencia de un "campo es-calar" cuando nos referimos a la densidad del mismo y a un
"campo vectorial" cuando nos referimos a las velocidades de las partículas del mismo.

OPERADORES DIFERENCIALES
Algunos de los más importantes son los siguientes;
∂
ƒ Diferencial respecto a una coordenada:
∂xi

ƒ Gradiente de un escalar: grad φ ≅ ∇ φ

ƒ Divergencia de un vector: div.v = ∇ v
ƒ Rotacional de un vector: rot.v ≡ ∇.v
ƒ Laplaciana de un escalar: ∇ .∇ φ ≡ ∇ 2φ

DERIVADAS DE FUNCIONES
Podemos definir las siguientes;
ƒ Derivada PARCIAL o LOCAL: Es la variación de una condición o propiedad con el
tiempo para un punto fijo en el espacio. En un punto fijo (P), podemos analizar la
variación de la velocidad de las partículas de fluido que circulan por aquel con el
transcurso del tiempo, esto da lugar a la "Derivada parcial de la velocidad del fluido" en
el punto (P), es decir, en símbolos;

⎛ ∂v ⎞
⎜ ⎟
⎝ ∂θ ⎠P

Página 7 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

ƒ Derivada TOTAL: Variación de una condición o propiedad con el tiempo en las

inmediaciones de un punto del espacio tomado como referencia. Para expresar el cambio
de la velocidad en los alrededores con respecto al tiempo, obtenemos la "Derivada Total
de la velocidad del fluido" para un desplazamiento infinitesimal, o sea;

dv ∂v ∂v dx ∂v dy ∂v dz
= + + +
d θ ∂θ ∂x d θ ∂y d θ ∂z d θ
ƒ Derivada MATERIAL, SUSTANCIAL ó SIGUIENDO EL MOVIMIENTO: Variación de una
condición o propiedad con el tiempo a lo largo de la trayectoria o movimiento, a partir de
un punto del espacio tomado como referencia. Para conocer el cambio de velocidad de
una partícula de fluido a lo largo del movimiento, escribimos la "Derivada Sustancial de la
velocidad del fluido" respecto del tiempo, que representa la aceleración o rapidez de la
variación de velocidad de la partícula en su trayectoria.

Dv ∂v ∂v ∂v ∂v
= + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z
Las siguientes figuras representan cada una de las respectivas derivadas enunciadas.

Punto fijo Po

⎛ ∂P ⎞
Propiedad varia según tiempo; “derivada parcial”= ⎜ ⎟
⎝ ∂ θ ⎠ xo , y o , z o
Punto desplazado PoÎP1
Propiedad varía según posición y velocidad de desplazamiento;

⎛ ∂P ⎞ ∂P ∂P ∂P
DERIVADA TOTAL: dP = ⎜ ⎟ dθ + dx + dy + dz
⎝ ∂ θ ⎠ ∂ x ∂ y ∂ z

Página 8 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

Punto desplazándose PoÎP1

Velocidad desplazamiento punto y vector posición iguales;

DP ∂P ∂P dx ∂P dy ∂P dy ∂P dz
DERIVADA SUSTANCIAL: = + + + +
Dθ ∂θ ∂x dθ ∂y dθ ∂y dθ ∂z dθ
dx dy dz
Siendo, vx = ; vy = ; vz =
dθ dθ dθ
DP ∂P ∂P ∂P ∂P
= + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z
Las expresiones generales anteriores pueden representarse para cada una de sus
componentes en el espacio. Por ejemplo, las componentes cartesianas del vector aceleración
serán;
D vx ∂vx ∂vx ∂vx ∂vx
ƒ Componente en "x": ax = = + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z
Dvy ∂v y ∂v y ∂v y ∂v y
ƒ Componente en "y": ay = = + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z
Dvz ∂vz ∂vz ∂vz ∂vz
ƒ Componente en "z": az = = + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z
Dv ∂v ∂v ∂v ∂v
ƒ Forma General: a= = + vx + vy + vz
Dθ ∂θ ∂x ∂y ∂z

Página 9 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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

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

Página 10 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

ƒ caras perpendiculares al eje z: ( ρ v ⎢ − ρ v ⎢ ) ∆x∆y

z z z z +∆z

El balance de materia quedará entonces;

∆ρ
∆θ
( ) ( )
( ∆x∆y∆z ) = ρvx ⎢x −ρvx ⎢x+∆x ∆z∆y + ρvy ⎢y −ρvy ⎢y+∆y ∆x∆z + ρvz ⎢− (
z ρvz ⎢z +∆z ∆x∆y )
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 ⎠

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.
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 iguala la velocidad
“neta” de entrada de densidad de flujo de materia dividido por el volumen.
Efectuando la diferenciación de los términos y reuniendo las derivadas de la densidad en
el primer miembro, se obtiene;

∂ρ ⎛ ∂ρ ∂ρ ∂ρ ⎞ ⎛ ∂v x ∂v y ∂v z ⎞
+ ⎜ vx + vy + vz ⎟ = − ρ ⎜ + + ⎟
∂θ ⎝ ∂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 una trayectoria de una partícula de fluido. Por lo tanto, la
expresión quedará;

Página 11 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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.
Estas ecuaciones son el resultado de la formulación de la conservación de la materia, y
puede efectuarse para un elemento de volumen de una forma arbitraria cualquiera.
Para ilustrarlo, vamos a aplicar el balance de masa a un elemento de volumen (V), que se
encuentra limitado por una superficie (S), como el de la figura siguiente.
La formulación del balance de materia es idéntica a la anterior, pero ahora la masa
acumulada en la unidad de tiempo será la integral extendida a todo el volumen;

∂
∂θ ∫v
ρ dv

Mientras que la masa que atraviesa la superficie

elemental puede expresarse por;
JJG
ρ v.N d s = ρ v c o s ϕ d s
Donde ϕ es el ángulo formado entre la velocidad
másica y el versor unitario dirigido “hacia fuera” en cada punto de la superficie.
La cantidad “neta” de materia que atraviesa toda la superficie del sistema viene dada
por el “flujo neto”, que se representa por;
JJG
∫(
s
)
ρ v.N ds

En consecuencia, el balance de materia indica que la suma de las expresiones para la

acumulación en el volumen y el flujo neto deben ser idénticamente nulas (iguales pero de
distinto signo), o sea;

∂ JJG
∂θ ∫
v
ρ dv + ∫(
s
ρ v.N )ds = 0

El teorema de la divergencia (o de Gauss), se puede escribir;

JG JJG JG
∫(
s
)
A. N ds = ∫(
v
)
∇ . A dv

Por lo tanto, aplicándolo a la situación planteada, se obtiene;

Página 12 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008
JJG
∫(
s
ρ v.N ) d s = ∫ (∇ . ρ v ) d v
v

que indica que el flujo neto de materia a través de la superficie del elemento de volumen
considerado tiene que ser igual a la velocidad con que se acumula en ese volumen,
quedando;

∂
∂θ ∫ ρ d v + ∫ (∇ .ρ v ) d v
v v
= 0

y para un volumen que no varía en el tiempo, se debe cumplir;

∂ρ
∫
v
∂θ
dv + ∫ (∇ .ρ v ) d v = 0
v

o lo que es lo mismo;

∂ρ
+ ( ∇. ρ v ) = 0
∂θ

Formas especiales importantes de la ecuación de continuidad son las siguientes:

ƒ 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;

( ∇.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.

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

Página 13 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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

Página 14 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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 “velocidad másica” (ρ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 las densidades de flujo de cantidad de movimiento pueden considerarse como
esfuerzos (τij). La convección indica que el primer subíndice indica la dirección de la
propagación y el segundo la de la velocidad. Para este caso, el segundo 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 + ∆ 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: τ z x ∆ x ∆ y ⎢z

ƒ Salida por z+ ∆ z: τ zx ∆ x ∆ y ⎢z + ∆ z

Página 15 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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 le la variación de cantidad de movimiento en el tiempo y espacio.

∆ρvx ρvxvx ⎢−
x ρvxvx ⎢x+∆x
ρv v ⎢ −ρvyvx ⎢y+∆y
= + y x y
∆θ ∆x ∆y
ρ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 + ρ vx v y + ρ vz v y ⎟ − ⎜ τ xy + τ yy + τ zy ⎟ − + ρ gy
∂θ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂x ∂y ∂z ⎠ ∂y

ƒ Componente “z”:

Página 16 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

∂ρ vz ⎛ ∂ ∂ ∂ ⎞ ⎛ ∂ ∂ ∂ ⎞ ∂p
= − ⎜ ρ vx vz + ρ vx 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 τ j 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 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:

∂ρ 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.
Para la expresión del término en dirección “x”, si aplicamos la ecuación de continuidad al
primer miembro y al primer termino del segundo miembro, la misma puede reordenarse,
quedando:

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θ

Página 17 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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.
∂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,
como tampoco distingue entre los fluidos compresibles o incompresibles.
Aplicando las relaciones entre esfuerzo-deformaciones (ver apéndice), aquellas
expresiones pueden transformarse para obtener las distribuciones de velocidad del fluido,
obteniendo ecuaciones generales que describen el movimiento del fluido.

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,

∂ρ
= −∇p −[∇.τ ] + ρg
∂θ

Página 18 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

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 la representación de la segunda Ley de Newton a un fluido en
movimiento.
Sin embargo, de acuerdo a las relaciones obtenidas entre esfuerzos y deformaciones
(ver Apéndice Páginas 21 a 26), podemos decir 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,

∂ρ v x ⎛ ∂ρ v x v x ∂ ρ v y vx ⎞ ⎛ ∂τ xx ∂τ yx ∂τ zx
= −⎜ + + ∂ρ v v
z x
⎟−⎜ + +
⎞ ∂p
⎟− + ρ gx
∂θ ⎝ ∂x ∂y ∂z ⎠ ⎝ ∂ x ∂ y ∂z ⎠ ∂ x

Que, combinada con la ecuación de continuidad,

∂ρ ⎛ ∂ρ vx ∂ ρ vy ∂ρ vz ⎞
= −⎜ + + ⎟
∂θ ⎝ ∂x ∂y ∂z ⎠
Permite obtener la componente x de la aceleración por unidad de volumen de fluido,

Dvx ∂ρ ∂τ xx ∂τ yx ∂τ zx ∂p
ρ = − −( + + ) − + ρgx
Dθ ∂x ∂x ∂y ∂z ∂x
Introduciendo las expresiones del esfuerzo del Apéndice, quedará,

Dvx ∂ρ ∂ ⎡ ∂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.
La expresión general de la Ecuación del Movimiento, en notación vectorial-tensorial
quedará entones,

Página 19 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

Dv µ
= −∇p − µ∇2v + ∇(∇.v) + ρ g
Dθ 3
Para el caso de fluidos de viscosidad constante. Como la viscosidad es función de la
temperatura principalmente y en pequeña medida de la presión, la anterior se puede aplicar
aún en los casos de moderados cambios de presión.

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 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. Sin embargo, las soluciones
generales de las Ecuaciones de Navier-Stokes aún no han sido encontradas debido a la
naturaleza no lineal de las ecuaciones diferenciales de derivadas parciales de segundo orden.
No obstante, un buen número de soluciones particulares se pueden obtener, y se
obtienen, en ciertos sistemas de flujo más simples o sencillos. Para ello debe disponerse de la
descripción física del sistema y usando las tablas en el sistema coordenado adecuado,
descartando los términos de valor cero. Este procedimiento tiene que hacer uso de una cierta
intuición para poder aproximarse a la solución.

ECUACION 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θ
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.

Página 20 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

ECUACION DE LA ENERGIA MECANICA

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 ínterconversió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,

∂
(½ρv2 ) = -(∇½ρv2v) - (∇.р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 v) : Velocidad neta de entrada de energía cinética debida al fluido global.
2
ƒ

ƒ - (∇.рv ) : Velocidad de trabajo producido por la presión de los alrededores sobre el

elemento de volumen.
ƒ - ρ(-∇.v ) : Velocidad de conversión “reversible” en energía interna.

ƒ - (∇.[τ .v]) : Velocidad de trabajo producido por las fuerzas viscosas que actúan sobre el

Página 21 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

elemento de volumen.

ƒ - (−τ : ∇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
elemento de volumen.

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.

APENDICE: TRANSPORTE MOLECULAR Y ESFUERZOS

Para fuerzas normales y tangenciales que actúan sobre un poco en un fluido serán
completamente definidas por nueve componentes tal como se muestra en la figura para las 3
caras positivas de un elemento cúbico en el espacio, cuyo centro está en el origen de
coordenadas.

Para las 3 caras negativas actúan las mismas ternas de esfuerzos pero en sentido
negativo. Estos no son independientes, ya que se demuestra que los esfuerzos cuyos
subíndices están permutados (p.ej. τxy= τyx) son iguales.
Para demostrarlo, mostraremos para una cara situada en un plano z=0, los esfuerzos
tangenciales que actúan sobre la misma en las direcciones x e y.

Página 22 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

La velocidad “neta” de deformación angular del elemento se la diferencia entre las

velocidades angulares dx y dy,

∂v y⎛ ∂v ⎞ ∂v y ∂v x
veloc .neta .def .ang . = −⎜− x ⎟ = +
∂x ⎝ ∂y ⎠ ∂x ∂y
Anteriormente se obtuvo la expresión de Newton que relacionaba la densidad de flujo de
cantidad de movimiento con la deformación angular en un fluido con una sola componente de
viscosidad. Luego, para el caso de una cara del elemento de volumen, bidimensional, la
relación será,

⎛ ∂v x ∂v y ⎞
τ xy = τ yx = − µ ⎜ + ⎟
⎝ ∂ y ∂x ⎠
⎛ ∂v y ∂v z ⎞
τ yz = τ zy = −µ ⎜ + ⎟
⎝ ∂ z ∂y ⎠
⎛ ∂v z ∂v x ⎞
τ xz = τ zx = − µ ⎜ + ⎟
⎝ ∂ x ∂z ⎠
Para los esfuerzos normales se llega a demostrar que,

∂vx 2
τ xx = −2µ + µ ( ∇ .v )
∂x 3
∂vy 2
τ yy = −2µ + µ ( ∇ .v )
∂y 3
∂vz 2
τ zz = −2µ + µ ( ∇ .v )
∂z 3

Página 26 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

Demostración
Suponiendo una masa m que gira alrededor de un eje con cierta distancia r con una
velocidad angular w, y a la que se aplica una fuerza tangencial F. Si consideramos que no
actúan otras fuerzas, la F tiende a aumentar la velocidad angular en dw.

El trabajo ejecutado será:

dw = Fds = Frd ϕ
dw = Md ϕ
Siendo:
ƒ M: momento respecto al eje

La energía cinética inicial será:

Ec = ½mv 2 = ½mr 2 w 2
Siendo:

ƒ v = rw
Y el incremento,

d ( Ec ) = mr 2 wdw
Si no hay rozamiento,

dw = d ( Ec)
Md φ = mr 2 wdw
dw
M = mr 2 w
dφ
Y como la aceleración angular (α) será:

dw dφ dw dw
α =w = =
dφ dθ dφ dθ
M = mr2α
M = Iα
Siendo:
ƒ I=mr2 el momento de inercia respecto al eje.

Página 27 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

2.C. ANALISIS DIMENSIONAL

Muchos problemas planteados en ingeniería, tanto en el flujo de fluidos como en la

transmisión del calor o la transferencia de masa, no se pueden solucionar por métodos
exclusivamente analíticos. Tampoco es posible obtener expresiones generales únicamente por
vía experimental, ya que cuando una función depende de varias variables, para obtener una
relación completa se deberán fijar todos los factores menos uno, para determinar entonces el
efecto de este y así sucesivamente, resultando extremadamente laborioso relacionar todos los
datos en un resultado único que comprenda todas las alternativas.
Como resultado de la experimentación, se puede llegar a ecuaciones empíricas donde no
se necesita prestar atención a su coherencia dimensional, es decir, que se pueden aplicar a
ecuaciones cuyos términos sean aparentemente de diferentes dimensiones. En estas
ecuaciones dimensionalmente heterogéneas, solo se pueden utilizar las variables en las
unidades que están especificadas.
Entonces, por ejemplo, si se estudia un fenómeno que depende de cinco variables, por el
método empírico habría que realizar por lo menos 100.000 experimentos para obtener series
de 10 ensayos donde cada variable está determinada por las otras cuatro fijadas en distintos
valores. Además de la cantidad de datos a correlacionar, generalmente algunas de las
variables son difíciles de manipular y el método resulta bastante limitado de aplicar para
obtener resultados amplios y generales.
Un método intermedio consiste en la aplicación del Análisis Dimensional, que se basa en
la proposición de que si hay una solución analítica que relacione todas las variables, la
expresión resultante deberá ser dimensionalmente homogénea, es decir, todos sus términos
tienen las mismas dimensiones, mientras que las constantes numéricas que ellas puedan
contener son adimensionales, resultando posible obtener la forma en que aparecen las
variables. Una ecuación de esta índole puede aplicarse con cualquier sistema de unidades,
siempre que estas sean coherentes entre sí, o lo que es lo mismo, utilizar las mismas para
cada magnitud fundamental.
Queda entonces determinado que, si en una ecuación de este tipo se dividen todos los
términos por uno de ellos, estos se transformarán en relaciones adimensionales. Una
combinación de variables tal que sus dimensiones se anulan, reciben el nombre de grupo,
número o razón adimensional y su valor es siempre el mismo, sin importar el sistema de
unidades que se utilice mientras sea homogéneo.
Las expresiones obtenidas solo conducen a estas relaciones entre las variables, no se
conocen los valores de las constantes numéricas, por lo que no se obtiene el resultado final
completo de la ecuación. Sin embargo, es un poderoso y valioso método que orienta a la

Página 28 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

experimentación mínima indispensable para corroborar los resultados, ya que reduce

drásticamente la cantidad de ensayos que deberían realizarse.
El Análisis Dimensional se puede aplicar directamente a las denominadas Ecuaciones de
Cambio (Ecuación de Continuidad, Ecuación del Movimiento) o puede estudiarse a partir de un
teorema fundamental (Teorema de Buckingham) siguiendo un proceso deductivo-inductivo que
se apoya en resultados experimentales previos y posteriores relacionados con el problema en
cuestión.

ANÁLISIS DIMENSIONAL DE LAS ECUACIONES DE CAMBIO Ó VARIACIÓN

Partiremos de las Ecuaciones de Continuidad y del Movimiento deducidas previamente,
para mayor sencillez en esta introducción nos limitaremos a aplicarlo en fluidos cuya densidad
y viscosidad son constantes, aunque debe destacarse que el procedimiento es general y no
tiene limitaciones en cuanto a estas condiciones.
Para un sistema de flujo cualquiera, es posible identificar arbitrariamente una longitud
característica "D" y una velocidad característica "V" del mismo, Por ejemplo, "D" es para flujos
en conductos cilíndricos el diámetro del tubo circular y para flujo alrededor de esferas el
diámetro de la misma, mientras "V" es la velocidad media en el tubo y la velocidad de
aproximación, respectivamente. Si bien la elección es arbitraria, deberá ser cuidadosamente
especificada en cada caso que se trate.
A partir de esa elección, se podrán definir una serie de variables adimensionales y
operaciones diferenciales, que en todos los casos cumplirán las mismas condiciones de
transformación, según las siguientes relaciones donde los símbolos tienen el significado
habitual;
v* = v/V ; p* = (p- p0 )/ρV2 ; Θ* = Θ V/D ;

x* = x/D ; y* = y/D ; z* = z/D ;

∇* = D∇ = δ1 (∂ /∂x*) + δ2 (∂ /∂y*) + δ3 (∂ /∂z*) ;

∇*2 = D2∇ 2 = ∂2 /∂x*2 + ∂2 /∂y*2 + ∂2 /∂z*2

(D /DΘ*) = (D/V) (D /DΘ)

Recordando las expresiones de las Ecuaciones de Cambio, para fluidos newtonianos de

densidad y viscosidad constante;
Ecuación de Continuidad: (∇ · v) = 0

Página 29 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

Ecuación del Movimiento:

ρ(Dv/DΘ)= -∇p + µ ∇2v + ρg
e introduciendo las variables adimensionales definidas anteriormente, estas ecuaciones
quedarán de la siguiente forma;
Ec. de Continuidad con variables adimensionales: 1/D (∇* · v*V) = 0

Ec. del Movimiento con variables adimensionales):

ρ (V /D)(D v*V /DΘ*) = - 1/D (∇* · p* ρV2)+ (µ /D2 ) ∇*2(v*V) + ρg

Reordenándolas y simplificando términos, se obtiene;

Ec. de Continuidad (adimensional): (∇* · v*) = 0

Ec. del Movimiento (adimensional):

(Dv* /DΘ*) = - (∇* p*)+ [µ /DVρ] ∇*2v* + [g D/V2] (g/g)

Las anteriores son las Ecuaciones de Variación equivalentes a las de partida, pero
presentadas con la forma de sus variables adimensionales, describiendo el movimiento del
fluido de igual manera que las expresiones dimensionales. En la Ecuación del Movimiento aquí
presentada, por otra parte, aquellas variables que describen el tamaño total, la velocidad del
sistema y las propiedades físicas se agrupan en un coeficiente que también es adimensional y
que se puede denominar como el "factor de escala" de la expresión.
De acuerdo a los Principios de Semejanza, si en dos sistemas diferentes sus factores de
escala son iguales para ambos, ambos están descritos por idénticas ecuaciones diferenciales
adimensionales. Además, si las condiciones adimensionales inicial y límite son las mismas, lo
que implica la restricción de Semejanza Geométrica, los dos sistemas son matemáticamente
idénticos, o en otras palabras, la distribución de la velocidad y de la presión adimensionales
son las mismas para cada uno de ellos, cumpliendo entonces la condición de Semejanza
Dinámica.
Para el pasaje de escala en procesos que no se conocen muy bien, deberían mantenerse
las restricciones de Semejanza Geométrica y Dinámica para obtener resultados comparables.
En este caso particular, como los factores de escala obtenidos intervienen con mucha
frecuencia en los estudios de ingeniería, reciben los nombres particulares de Número de
Reynolds y Número de Froude, en honor a los celebres investigadores en el estudio de la

Página 30 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

mecánica de los fluidos. Cabe destacar que en las expresiones se obtienen las inversas de los
mismos, que relacionan las fuerzas inerciales, las fuerzas viscosas y las gravitatorias, según;

ρ v D Fuerzas Inerciales
Número de Reynolds: Re = ≡
µ Fuerzas Viscosas

Partiendo de la expresión de la Ecuación del Movimiento en función de las presiones

locales y la gravitación, se hubiera obtenido un segundo factor de escala, que es la inversa del
denominado Número de Froude, que relaciona las fuerzas inerciales con las gravitatorias,
según;

v2 Fuerzas Inerciales
Número de Froude: Fr = ≡
g D Fuerzas Gravitatorias

ANÁLISIS DIMENSIONAL POR EL TEOREMA DE BUCKINGHAM

Tal como se presentó anteriormente en este curso, ya en el año 1850 Stokes estableció
que en un sistema de flujo geométricamente semejante puede utilizarse el número de Reynolds
como criterio de semejanza dinámica. En 1914, Buckingham formuló el denominado teorema
de π, que es el fundamento de la descripción siguiente. Según Bridgman, los principios
fundamentales del método son los siguientes:
1) Todas las magnitudes físicas pueden expresarse como funciones potenciales de
un reducido número de magnitudes fundamentales.
2) Las ecuaciones que relacionan las magnitudes físicas son homogéneas desde
un punto de vista dimensional, es decir, todos los términos tienen las mismas dimensiones.
3) Teorema π de Buckingham; si una ecuación es dimensionalmente homogénea,
puede reducirse a una serie completa de números adimensionales, en los que figuran todas
las variables físicas que influyen en el fenómeno, las constantes dimensionales que puedan
corresponder al sistema de unidades elegido (si este fuera redundante) y las constantes
universales que puedan intervenir en el fenómeno que se trate.

En forma muy general y según algunos autores, solo existirían cinco constantes
universales imprescindibles; (i) la velocidad de la luz, (ii) la constante de Planck, (iii) la
constante de gravitación, (iv) la constante dieléctrica y (v) la constante de Boltzmann. En
general en el campo de la Ingeniería Química y en especial en el de los Fenómenos de
Transporte, las anteriores no suelen intervenir habitualmente, aunque sí las magnitudes físicas

Página 31 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

fundamentales, tales como la masa (M), la longitud (L), el tiempo (θ) y la temperatura (T), entre
otras.
Así, una descripción completa de un fenómeno estará representada por la función;
Ŧ (N1a, N2b, N3c, . . . Nik) = 0
donde (N1, N2, N3, . . . Ni ) son las variables o constantes dimensionales y (a, b, c, . . .k) los
exponentes a los que deben estar elevadas para obtener la solución del problema.
El teorema π de Buckingham del tercer apartado es una consecuencia de los dos
anteriores, teniendo una demostración matemática que excede los alcances del curso y de
estas notas. Plantea que una serie de números adimensionales es completa cuando todos ellos
son independientes entre sí, de manera que cualquier otro que pueda formarse con las mismas
variables sería una combinación de dos o más números adimensionales de los contenidos en
la serie completa original. En forma restringida, la cantidad de números adimensionales π
independientes, que pueden formarse con N variables y constantes dimensionales es igual a su
diferencia con la cantidad de magnitudes fundamentales F del sistema de unidades elegido, es
decir;
π=N-F
En una forma menos restringida, la relación anterior se puede expresar como (π = N - j),
donde j será la característica de una matriz formada por los exponentes en que están elevadas
las magnitudes fundamentales en las ecuaciones dimensionales de las distintas variables y
constantes dimensionales con relación a un determinado sistema de unidades, o sea, el orden
del mayor determinante distinto de cero que puede formarse a partir de dicha matriz. En la
mayoría de los casos, j coincide con F.
En conclusión, siendo π1, π2, π3, . . . πn, una serie completa de relaciones
adimensionales que pueden formarse con las variables que determinan un cierto fenómeno, el
Teorema π de Buckingham se puede representar por,
ƒ (π1, π2, π3, . . . πn) = 0
donde ƒ es una función universal independiente del sistema de unidades elegido.

Como ocurre en muchos casos, si su resolución en términos de funciones conocidas

resulta muy difícil, se utilizan formas de representación gráfica o de aproximaciones.
Si bien existen dos métodos principales para resolver un sistema aplicando el teorema,
identificados como métodos de Buckingham y de Rayleigh, en nuestras aplicaciones
utilizaremos una forma simplificada propuesta por Kay (considerando que j ≡ F) , siguiendo
ordenadamente los siguientes pasos;
a) Identificar y establecer todas las variables (N) que condicionan el fenómeno en
estudio.

Página 32 de 33
 Fenómenos de Transporte – Unidad Temática 2
Revisión: Junio 2008

b) Formular el sistema de magnitudes fundamentales (F) que se desee. Esta

elección deberá contener no menos de todas las magnitudes que intervienen en el
fenómeno. Ocasionalmente, puede ser redundante conteniendo magnitudes derivadas de
las fundamentales.
c) Calcular la cantidad de grupos adimensionales (π = N - F) que describirán el
fenómeno.
d) Separar una cantidad π del total de N variables, teniendo en cuenta que cada
magnitud fundamental deberá aparecer por lo menos una vez entre las elegidas (en este
método simplificado, el proceso de selección es por tanteo, necesitando conocer
previamente algo de la naturaleza física del fenómeno). Estas variables serán conocidas
como "características", porque aparecerán solamente en uno de los números
adimensionales que se formen (el que la caracteriza).
e) Formular ecuaciones entre cada variable "característica" y todas las restantes
variables "no características". Deben quedar conformadas π ecuaciones, en cada una se
encuentra solo una variable característica, un grupo adimensional y todas las variables
no-características elevadas a un exponente particular.
f) Reemplazar en cada una las variables por sus dimensiones (no unidades) y
resolver entre cada variable "característica" y todas las restantes variables "no
características". Formar una matriz de los exponentes a que quedan afectadas las
magnitudes fundamentales en las ecuaciones dimensionales de las variables. Al resolver,
se debería obtener el valor del exponente que afecta a cada variable.
g) Introducir los resultados en el conjunto de expresiones del apartado "e". Se
obtendrán las expresiones de los números adimensionales (π1, π2, π3, . . . πn ) formado
por productos, cocientes o potencias de las variables que intervienen en la solución del
problema. No se obtienen las constantes numéricas adimensionales de las ecuaciones,
que deberán determinarse por experimentación.

En conclusión, se puede verificar una significativa reducción en la cantidad de ensayos a

realizar para obtener una correlación que describa el fenómeno. Así, en el ejemplo inicial de
cinco variables, si fueran tres las dimensiones fundamentales consideradas la descripción se
lograría con dos grupos adimensionales independientes, por lo que se deberían realizar
solamente 100 ensayos con las mismas premisas anteriores.

Página 33 de 33