Mecánica de Fluidos: Notas de Clase
Mecánica de Fluidos: Notas de Clase
net/publication/334451380
CITATIONS
0
1 author:
Francisco J. Valdes-Parada
Metropolitan Autonomous University
132 PUBLICATIONS 1,119 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Francisco J. Valdes-Parada on 14 July 2019.
2
Notas de clase: Francisco J. Valdés Parada
Índice
Página
1. Introducción 7
1.1. Generalidades sobre el modelado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.2. Clasificación de modelos matemáticos . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.3. Construcción de modelos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4. Un primer vistazo a la mecánica de fluidos . . . . . . . . . . . . . . . . . . . . . . . . 10
1.5. Vectores y notación indicial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.6. Producto punto o escalar . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.7. Producto cruz o vectorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.8. Producto diádico y tensores de segundo orden . . . . . . . . . . . . . . . . . . . . . . 16
1.9. Propiedades de los productos diádicos . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.10. Operadores diferenciales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.11. Hipótesis del continuo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.12. Tres tipos de derivadas temporales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
1.12.1. Derivadas temporales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
1.13. Teoremas integrales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.13.1. Teorema de la divergencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.13.2. Teorema general del transporte . . . . . . . . . . . . . . . . . . . . . . . . . . 30
1.14. Conservación de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.15. Transporte en la interfase . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35
1.15.1. Condición de salto de masa en una interfase . . . . . . . . . . . . . . . . . . . 36
1.16. Lı́neas de camino, lı́neas de corriente y lı́neas de marcha . . . . . . . . . . . . . . . . 38
1.16.1. Lı́neas de camino . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
1.16.2. Lı́neas de corriente . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.16.3. Lı́neas de marcha . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
3
Notas de clase: Francisco J. Valdés Parada
2.7.2. Fluidos Newtonianos generalizados . . . . . . . . . . . . . . . . . . . . . . . . 59
2.8. Tipos de flujo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
2.8.1. Flujo reptante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.8.2. Flujo potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.9. La ecuación de energı́a mecánica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
3. Hidrostática 67
3.1. Ecuación fundamental de la hidrostática . . . . . . . . . . . . . . . . . . . . . . . . . 67
3.2. Presión hidrostática en lı́quidos y gases . . . . . . . . . . . . . . . . . . . . . . . . . 68
3.3. Manometrı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.4. Fuerzas sobre superficies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75
4
Notas de clase: Francisco J. Valdés Parada
7. Flujo reptante 135
7.1. Flujo reptante alrededor de una esfera . . . . . . . . . . . . . . . . . . . . . . . . . . 135
10.Turbulencia 159
10.1. Preliminares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160
10.2. Suavizado en el tiempo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 160
10.3. Problema de cerradura . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 162
Referencias 190
5
Notas de clase: Francisco J. Valdés Parada
6
Notas de clase: Francisco J. Valdés Parada
1 | Introducción
Contenido
1.1. Generalidades sobre el modelado . . . . . . . . . . . . . . . . . . . . 7
1.2. Clasificación de modelos matemáticos . . . . . . . . . . . . . . . . . 8
1.3. Construcción de modelos . . . . . . . . . . . . . . . . . . . . . . . . 9
1.4. Un primer vistazo a la mecánica de fluidos . . . . . . . . . . . . . . 10
1.5. Vectores y notación indicial . . . . . . . . . . . . . . . . . . . . . . . 12
1.6. Producto punto o escalar . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.7. Producto cruz o vectorial . . . . . . . . . . . . . . . . . . . . . . . . 15
1.8. Producto diádico y tensores de segundo orden . . . . . . . . . . . . 16
1.9. Propiedades de los productos diádicos . . . . . . . . . . . . . . . . . 17
1.10. Operadores diferenciales . . . . . . . . . . . . . . . . . . . . . . . . . 20
1.11. Hipótesis del continuo . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
1.12. Tres tipos de derivadas temporales . . . . . . . . . . . . . . . . . . . 26
1.12.1. Derivadas temporales . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
1.13. Teoremas integrales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.13.1. Teorema de la divergencia . . . . . . . . . . . . . . . . . . . . . . . . . 28
1.13.2. Teorema general del transporte . . . . . . . . . . . . . . . . . . . . . . 30
1.14. Conservación de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.15. Transporte en la interfase . . . . . . . . . . . . . . . . . . . . . . . . 35
1.15.1. Condición de salto de masa en una interfase . . . . . . . . . . . . . . . 36
1.16. Lı́neas de camino, lı́neas de corriente y lı́neas de marcha . . . . . . 38
1.16.1. Lı́neas de camino . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
1.16.2. Lı́neas de corriente . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
1.16.3. Lı́neas de marcha . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
7
Notas de clase: Francisco J. Valdés Parada
rediseño. La operación consiste en utilizar un equipo ya existente que represente un cierto aspecto de
la realidad que se desea modelar bajo condiciones controladas. Esta actividad es posterior al diseño,
en la cual se deben construir los elementos necesarios para conformar el equipo experimental que,
al operarlo, representará algún fenómeno a estudiar. Por supuesto, los diseños son mejorables para
abarcar, por ejemplo, una gama más amplia de situaciones a estudiar. A la actividad de mejorar
un diseño, la llamaremos rediseño. De lo anterior puede deducirse que las actividades de operación,
diseño y rediseño no están aisladas pues para operar correctamente un equipo es necesario conocer
su diseño, lo cual puede llevar a encontrar puntos de mejora que conduzcan a un rediseño del equipo.
8
Notas de clase: Francisco J. Valdés Parada
y los modelos basados en principios fundamentales son modelos a priori de la experiencia.
1. ¿Cuáles atributos, que son cruciales en la caracterización del objeto de estudio, se pueden
incorporar al modelo?
2. ¿Cuáles son los principios y las relaciones sobre las que se sostendrá el modelo y qué tipo de
ecuaciones resultan de dichos fundamentos?
Las preguntas 1, 3 y 4 están relacionadas con la necesidad de entender con suficiente claridad el
sistema que se está estudiando. Lo cual puede requerir la construcción de un esquema o incluso de
un modelo fı́sico. Por su parte, las preguntas 2-5 están acopladas, pues la información que aportará
el modelo depende de las leyes que lo fundamentan, la geometrı́a, el papel que juega el tiempo
y las suposiciones adoptadas. Las respuestas a las preguntas 3 y 4 son cruciales pues a menudo
constituyen el punto de retorno cuando un modelo no describe satisfactoriamente a un fenómenos
de la realidad.
Además de estas cuestiones, Aris (1976) propuso una serie de Máximas para el modelado
matemático, las cuales listamos a continuación:
2. Escoja una notación que le resulte agradable, pero no se vuelva adicto a ella.
3. Vuelva las variables adimensionales, ..., pero no pierda de vista a aquellas cantidades que
tengan que ser modificadas más adelante en el problema, no olvide el origen fı́sico de cada
parte.
4. Use (siempre que sea posible) lı́mites a priori de origen fı́sico o matemático para mantener a
todas las variables en el mismo orden de magnitud.
6. Use métodos cortos pero eficaces, sin llevarlos más allá de su punto de utilidad.
8. Estudie casos lı́mite y determine cómo se enlazan con problemas más simples que pueden
resolverse explı́citamente.
9
Notas de clase: Francisco J. Valdés Parada
9. Use aproximaciones burdas, ..., reflexione sobre las analogı́as que sugieren, pero recuerde sus
limitaciones.
10. Rearregle el problema, ..., Esté preparado para trabajar con soluciones implı́citas.
Como puede notarse, Aris formula sugerencias que no se limitan al planteamiento del modelo, sino
también a su análisis y reformulación. Por supuesto, ésta es una visión particular del autor y por
ello la relevancia de la máxima número 13. Para una discusión más detallada de la filosofı́a y el
procesos de modelado matemático asociado a fenómenos de transporte se recomienda al lector el
texto de Valdés-Parada (2014).
10
Notas de clase: Francisco J. Valdés Parada
tangenciales (fuerzas viscosas). La teorı́a del flujo potencial aplica sin mayor problemas a flujos de
fluidos suficientemente veloces y lejos de superficies sólidas. Este último punto es importante, pues
dicha teorı́a, originalmente propuesta por Euler, gobernaba la hidráulica a finales del siglo XIX. El
problema de esta teorı́a es que, al no considerar las fuerzas de rozamiento es que predice que un
cuerpo puede moverse infinitamente bajo la acción de una fuerza inercial (paradoja de D’Alembert).
Esta problemática condujo a Navier y a Stokes a proponer la inclusión de un término que incluyera
dichas fuerzas, con lo cual se pasó de la hidráulica empı́rica a la hidrodinámica. Sin embargo, el
problema que presentan las ecuaciones de Navier-Stokes es que sólo pueden ser resueltas de manera
exacta para condiciones de flujo no inercial. De hecho, aún si el flujo del fluido es tremendamente
veloz, si no existe deslizamiento, su velocidad debe reducirse a cero en la superficie del sólido, por
lo que existe una pequeña región cerca de dicha superficie, conocida como capa lı́mite, donde las
fuerzas viscosas son tan importantes como las inerciales. Fuera de dicha capa, las fuerzas inerciales
dominan sobre las fuerzas viscosas y el flujo es potencial. Este resultado que vino a resolver la
paradoja de D’Alembert se le debe a Ludwig Prandtl. De hecho, el concepto de capa lı́mite es muy
útil para comprender la diferencia entre el flujo laminar y el flujo turbulento. En el flujo laminar
la energı́a cinética del fluido fuera de la capa lı́mite se balancea con las pérdidas por disipación
viscosa que ocurren en la capa lı́mite y por ello se pueden apreciar lı́neas rectas de fluido. Cuando
el flujo es turbulento este balance se pierde pues es mayor la energı́a cinética que tiene el fluido
que la tasa a la cual se disipa dicha energı́a por los esfuerzos viscosos en la capa lı́mite. En el flujo
turbulento se generan remolinos en cascada donde la energı́a cinética ingresa en los remolinos de
mayor tamaño y se transmite a los remolinos de escala menor donde se disipa de manera viscosa.
Es importante mencionar que, al aplicar estos principios de conservación, se obtienen modelos
abiertos pues tienen más incógnitas que ecuaciones. Para atender esta situación, se han desarrollado
(de manera experimental y teórica) ecuaciones de cerradura que completan los modelos basados en
los principios de conservación. En especı́fico, en el caso de la mecánica de fluidos, las ecuaciones de
cerradura que se necesitan son aquellas que proporcionan la relación entre la fuerza (o más bien el
esfuerzo, que es la fuerza por unidad de área) que se aplica y la deformación que sufre el fluido.
Como se verá más adelante, cuando esta relación es lineal, se sigue la ley de Newton de la viscosidad
y existen variantes interesantes de esta ley que se discutirán en su momento.
Por último, las ecuaciones de Navier-Stokes aplican al nivel de escala del continuo, del
cual hablaremos más adelante, y se conocen como ecuaciones locales instantáneas. En muchas
aplicaciones es necesario llevar a cabo un proceso de escalamiento para filtrar la información
redundante de las escalas más pequeñas y desarrollar modelos matemáticos válidos a un nivel
de escala macroscópico. De esta forma se podrán dimensionar equipos para el flujo de fluidos en
tuberı́as y arreglos de tuberı́as, medios porosos, entre muchas otras aplicaciones. Cabe agregar que
a nivel de escala del continuo es sólo necesario considerar las ecuaciones de conservación de masa
y de cantidad de movimiento. Sin embargo, a escala macroscópica debe considerarse además la
ecuación de la energı́a mecánica. Lo anterior se debe a que a la escala del continuo, la ecuación de
energı́a mecánica es simplemente el resultado del producto punto de la ecuación de Navier-Stokes
con la velocidad. Sin embargo, dicha relación tan simple no se conserva, en general, en aplicaciones
macroscópicas.
Los anteriores párrafos resumen el contenido de texto que involucra primeramente el desarrollo
de las ecuaciones que gobiernan el flujo de fluidos, para posteriormente presentar aplicaciones
que van desde las condiciones hidrostáticas, pasando por los régimenes de flujo reptante, laminar,
turbulento y potencial, para concluir con los balances macroscópicos de cantidad de movimiento.
11
Notas de clase: Francisco J. Valdés Parada
Curiosamente, las ecuaciones a la escala del continuo y en la macroescala comparten el mismo
punto de origen que son los principios fundamentales de conservación de masa y de cantidad de
movimiento. De hecho, dichos principios de conservación son macroscópicos y de ellos se extraen
las ecuaciones correspondientes en la escala del continuo. Antes de comenzar a abordar estos temas
es importante hablar un mismo lenguaje, por ello, a continuación se presenta un breve repaso de
vectores y tensores que será de utilidad para el resto del curso.
En la ec. (1.1) ai indica cualquiera de las tres componentes del vector a, y para abreviar la escritura
de esta expresión se introduce a continuación la convención de la sumatoria. Dicha convención
establece que la existencia de términos con ı́ndice repetido dos veces indica la suma de todos los
términos posibles. Por ejemplo, para el caso tridimensional, la suma incluye tres términos y, usando
la convención de la sumatoria, la ec. (1.1) se escribe como:
a = ai ei (1.2)
Bajo esta convención un ı́ndice que no se encuentra repetido en un término, tomará cualquier valor
permitido por la dimensión del espacio vectorial. Ası́, para un sistema tridimensional, se tiene que
Se debe hacer notar que cuando en un término aparezca un ı́ndice repetido más de dos veces, ésto
no implicará sumatoria. Por ejemplo,
Es importante mencionar que, a diferencia de los escalares, la división de vectores no está definida.
12
2 3 1 2
En la ec. (1.1)
y que representa ai indica cualquiera
la definición normalmente de las tres componentes del vector a , y para abreviar la
usada
escritura de esta expresión se introduce a continuación la convención de la sumatoria.
Dicha convención establece que a ·lab existencia
= a1 b1 + ade 2 b2términos
+ a3 b3 con índice repetido dos veces (1.7)
Note queindica
en ellacaso
sumaparticular
de todos los términos posibles. Por ejemplo: Para el caso tridimensional, la
suma incluye tres términos y usando a ·laa convención
= ai ai = |ade |2 la sumatoria la ec. (1.1) se escribe (1.8)
como:
por lo que la magnitud del vector a, está definida como
a =√
ai ei (1.2)
|a | = ai ai (1.9)
Una interpretación geométrica de este producto puede encontrarse aplicando la ley de los cosenos
a un triángulo como el mostrado en la Figura 1.2.3 El cuadrado de la longitud
Página J.A. del lado
Ochoa c = b−a
Tapia
puede escribirse como
|c |2 = |b − a |2 (1.10)
donde
|b − a |2 = (b − a) · (b − a) = |b |2 − 2(a · b) + |a |2 (1.11)
De la ley de los cosenos se tiene que
Por lo que, al igualar las ecs. (1.10) y (1.12), tomando en cuenta la ec. (1.11), resulta que
13
Notas de clase: Francisco J. Valdés Parada
En donde θab indica el ángulo definido por los vectores a y b. Debe notarse que el coseno es una
función simétrica o par, y por lo tanto
Como consecuencia
a·b=b·a (1.15)
Sin embargo, esta propiedad no es extensiva a tensores como se verá más adelante. El resultado
dado en la ec. (1.13) puede usarse para encontrar el valor de ei · ej ,
En donde θei ej es el ángulo formado por ei y ej . Como éstos son vectores unitarios su magnitud es
la unidad, y por lo tanto la ec. (1.16) se reduce a
Note que cos θei ej es igual a 1 sı́ i = j, e igual a 0 sı́ i 6= j. Esto es fácil de visualizar si se recuerda, de
cursos anteriores, que el producto escalar de un vector a con un vector unitario ei , es la proyección
de a en la dirección ei . Por lo tanto, sı́ a tiene la misma dirección que ei , su proyección es |a |, y
en contraste sı́ a es perpendicular a ei , su proyección es nula. Podemos ahora resumir los valores
de ei · ej en la siguiente forma
1, sı́ i = j
ei · ej = (1.18)
0, sı́ i 6= j
Es conveniente introducir ahora el sı́mbolo conocido como delta de Kronecker
1, sı́ i = j
δij = (1.19)
0, sı́ i 6= j
ei · ej = δij (1.20)
a · b = ai ei · bj ej = ai bj δij (1.21)
Resultado que, al utilizar la definición de la delta de Kronecker (ec. 1.19), se reduce a la ec. (1.6).
b
c
a
14
Notas de clase: Francisco J. Valdés Parada
1.7. Producto cruz o vectorial
El producto vectorial o cruz está definido por el determinante
e1 e2 e3
a × b = a1 a2 a3 (1.22)
b1 b2 b3
que se puede desarrollar para obtener
a × b = (a2 b3 − a3 b2 )e1 + (a3 b1 − a1 b3 )e2 + (a1 b2 − a2 b1 )e3 (1.23)
Para abreviar la escritura de la ec. (1.23) se introduce el sı́mbolo de permutación ijk definido de
la siguiente manera
0, si alguno de los ı́ndices está repetido
ijk = 1, si ijk es una permutación par de 123 (1.24)
−1, si ijk es una permutación non de 123
La permutación par de los ı́ndices ijk se logra cuando de los ı́ndices originales se encuentra la
secuencia 123 usando un número par de permutaciones. Es sencillo decidir si una permutación es
par o non si los ı́ndices ijk dan la secuencia 123 en una u otra de las maneras indicadas en la
Figura 1.3. Como puede notarse, ijk = +1 cuando los ı́ndices son 123 en la secuencia siguiendo
la dirección de las manecillas del reloj, y ijk = −1 cuando 123 están colocados en una secuencia
contraria a las manecillas del reloj. Por lo tanto
123 = 231 = 312 = +1
(1.25)
132 = 321 = 213 = −1
Usando el sı́mbolo de permutación y sus propiedades, no es difı́cil demostrar de la ec. (1.23),
que el producto vectorial entre a y b se escribe en notación indicial como
a × b = ijk ai bj ek (1.26)
La interpretación geométrica del producto vectorial a × b se puede encontrar a partir del uso
de |a × b |2 = (a × b) · (a × b) que lleva a
a × b = |a | |b |senθab e⊥ , para 0 ≤ θab ≤ π (1.27)
1 1
+
3 2 3 2
a) b)
15
Notas de clase: Francisco J. Valdés Parada
En donde θab es el ángulo definido por los vectores a y b. El vector unitario e⊥ es perpendicular al
plano que forman los vectores a y b, y con dirección de acuerdo a un sistema coordenado de mano
derecha. Note, de la ec. (1.27), que el módulo del vector a×b corresponde al área del paralelogramo
definido por los vectores a y b.
Por otro lado, una consecuencia directa de las propiedades de un determinante usadas en la
definición dada en la ec. (1.22) es que
a × b = −b × a (1.28)
También se puede llegar a esta conclusión a partir de la ec. (1.26) debido a las propiedades del
sı́mbolo de permutación. Ası́, consecuencia también de cualquiera de las ecs. (1.22), (1.26) ó (1.27),
es que los productos vectoriales entre los vectores unitarios ei , son
e1 × e2 = −e2 × e1 = e3
e2 × e3 = −e3 × e2 = e1
(1.29)
e3 × e1 = −e1 × e3 = e2
e1 × e1 = e2 × e2 = e3 × e3 = 0
ó
c1 c2 c3
c · a × b = a1 a2 a3 (1.31)
b1 b2 b3
Es posible, usando la ec. (1.27) y su interpretación geométrica, demostrar que |c·a×b | corresponde
al volumen del paralelepı́pedo definido por los vectores a, b y c. Ası́, si los vectores a, b y c son
coplanares, el producto c · a × b es cero.
Esta operación es completamente diferente a los productos escalar y vectorial puesto que cada uno
de los términos resultante tiene asociado dos vectores unitarios, este resultado es lo que se denomina
producto diádico o tensor de segundo orden. En forma expandida el producto ab, que se indicará
con A, es
A = ab = a1 b1 e1 e1 + a1 b2 e1 e2 + a1 b3 e1 e3
+ a2 b1 e2 e1 + a2 b2 e2 e2 + a2 b3 e2 e3 (1.33)
+ a3 b1 e3 e1 + a3 b2 e3 e2 + a3 b3 e3 e3
16
Notas de clase: Francisco J. Valdés Parada
Como puede notarse, un tensor de segundo orden tiene nueve componentes y dos direcciones para
cada una de sus componentes.
El producto ba, denominado B, es
B = ba = b1 a1 e1 e1 + b1 a2 e1 e2 + b1 a3 e1 e3
+ b2 a1 e2 e1 + b2 a2 e2 e2 + b2 a3 e2 e3 (1.34)
+ b3 a1 e3 e1 + b3 a2 e3 e2 + b3 a3 e3 e3
De la comparación de las ecuaciones (1.33) y (1.34) se puede ver que ab 6= ba, ésto porque, en
general, a2 b1 6= b2 a1 , a3 b1 6= b3 a1 , etc. La escritura del tensor A se puede simplificar utilizando la
nomenclatura
Aij = ai bj (1.35)
De tal forma que la ec. (1.33) es
A = Aij ei ej (1.36)
Donde Aij son los elementos del producto diádico o tensor de segundo orden.
En este texto se adopta la convención basada en los vectores unitarios más próximos.
Nótese que A · c = c · A solamente sı́ Aij = Aji , si el tensor tiene esta propiedad se
denomina tensor simétrico.
AT = ATij ei ej (1.40)
17
Notas de clase: Francisco J. Valdés Parada
en donde
ATij = Aji (1.41)
I = δij ei ej (1.42)
Note que es un tensor simétrico y que cumple con las siguientes fórmulas:
a·I=a (1.43)
A·I=A (1.44)
AT = −A (1.45)
f) El doble producto escalar de dos tensores de segundo orden también está sujeto a la
misma convención de anidamiento
A : B = (Aij ei ej ) : (Bkl ek el ) = Aij Bkl ei · (ej · ek )el = Aij Bkl (ei · el )(ej · ek )
= Aij Bkl δil δjk = Aij Bji (1.47)
18
Notas de clase: Francisco J. Valdés Parada
Ejemplo 1.9-1: Demuestre la siguiente identidad
c · A = AT · c
Solución: Usando notación indicial y tomando en cuenta el resultado del producto escalar
en la ec. (1.39) se tiene que
c · A = (Aij ci )ej
Más aún, el uso de la definición del transpuesto de un tensor dada en la ec. (1.41) permite
encontrar
c · A = (ATji ci )ej = (ATij cj )ei = AT · c
De este resultado no es difı́cil observar que si un tensor es simétrico A · c = c · A. Los tensores
simétricos aparecen frecuentemente en el estudio de mecánica de fluidos.
Ejemplo 1.9-2: Demuestre que el primer término en la Ec. (1.51) es un tensor simétrico y que
el segundo es un tensor antisimétrico.
Solución: Al combinar las ecs. (1.49)-(1.51), se obtienen las siguientes expresiones para los
tensores B y C
1
B= (Aij + Aji )ei ej (1.52a)
2
1
C = (Aij − Aji )ei ej (1.52b)
2
Dado que Aij + Aji = Aji + Aij y Aij − Aji = −(Aji − Aij ), entonces
1
B = (Aji + Aij )ei ej = BT (1.53a)
2
1
C = − (Aji − Aij )ei ej = −CT (1.53b)
2
Ejemplo 1.9-3: Demuestre que el doble producto escalar de un tensor simétrico con uno
antisimétrico es cero.
Solución: Para llevar a cabo esta demostración, suponga que el producto A : B es igual a
w, donde A es un tensor simétrico y B es antisimétrico. Usando las representaciones de estos
tensores en notación indicial se encuentra
Esta última expresión implica que w = −w, y por ello w = 0 es la única posibilidad de satisfacer
la igualdad. Por lo tanto, A : B = 0, sı́ A = AT y B = −BT .
19
Notas de clase: Francisco J. Valdés Parada
1.10. Operadores diferenciales
Ahora se mostrará la representación de operadores diferenciales en sistemas coordenados
Cartesianos usando notación indicial y la convención de la sumatoria. El operador nabla, ∇, es
un operador vectorial diferencial que se define como
∂ ∂ ∂ ∂
∇ = e1 + e2 + e3 = ei (1.54)
∂x1 ∂x2 ∂x3 ∂xi
Por lo que el gradiente de un escalar, ϕ, no requiere mayor manipulación y es
∂ϕ
∇ϕ = ei (1.55)
∂xi
Sin embargo, el gradiente del vector v = vi ei da un tensor de segundo orden como se muestra a
continuación
∂ ∂vj
∇v = ei (vj ej ) = ei ej (1.56)
∂xi ∂xi
∂v
Es claro que cada uno de los nueve componentes del tensor ∇v es ∂xji . Note además que la definición
del tensor identidad puede recuperarse a partir del gradiente del vector de posición r:
I = ∇r (1.57)
BT 1 = Bjikei ej ek (1.59)
BT 2 = Bikjei ej ek (1.60)
BT 3 = Bkjiei ej ek (1.61)
Para interpretar geométricamente al operador gradiente es conveniente recordar primero que
un vector se representa como una flecha cuyo largo corresponde a su magnitud y la flecha apunta
en la dirección del vector. De esta forma, el vector de la velocidad de un fluido representa con su
orientación la dirección del flujo en un punto dado y su largo representa la magnitud de la velocidad.
Note que la velocidad tiene tres componentes y la dirección nos indica cuál de ellas es predominante.
Aplicando estas ideas al operador gardiente de un escalar, ϕ, resulta entonces que la dirección de
este vector corresponde a la dirección en la cual ϕ experimenta el máximo cambio respecto a la
posición en un punto dado. De la misma forma, la magnitud del gradiente de una función escalar
representa la máxima razón de cambio de ϕ con la posición en un punto. Esta interprtación puede
extenderse a vectores y tensores.
El operador de divergencia, es el producto escalar del operador nabla aplicado a un vector
∂ ∂vi ∂v1 ∂v2 ∂v3
∇ · v = ei · (vj ej ) = = + + (escalar) (1.62)
∂xi ∂xi ∂x1 ∂x2 ∂x3
20
Notas de clase: Francisco J. Valdés Parada
Mientras que, para un tensor, la divergencia es
∂ ∂Aik
∇ · A = ei · (Ajk ej ek ) = ek (vector) (1.63)
∂xi ∂xi
∂A11 ∂A21 ∂A31 ∂A12 ∂A22 ∂A32 ∂A13 ∂A23 ∂A33
= + + e1 + + + e2 + + + e3
∂x1 ∂x2 ∂x3 ∂x1 ∂x2 ∂x3 ∂x1 ∂x2 ∂x3
El significado fı́sico de este operador es muy distinto que para un gradiente. Para entender este
operador considere un cubo cuyos lados miden ∆x, ∆y y ∆z, respectivamente. Consideremos por un
momento al vector de velocidad y supongamos entonces que dentro de ese cubo hay un movimiento
de fluido. Entonces, el flujo volumétrico que sale de las superficies del cubo se puede calcular como
m3
Z
Q̇ = n · v dA[=] (1.64)
s
∆A
donde n es el vector unitario perpendicular a cada cara del cubo y siempre está dirigido de la
superficie hacia afuera. Note además que
∆V = ∆x∆y∆z (1.66)
Para poder avanzar en nuestros desarrollos es conveniente desarrollar tanto el producto punto como
las integrales de superficie para obtener:
Z Z
n · v dA = (nx vx + ny vy + nz vz ) dA (1.67)
∆A ∆A
y+∆y
Z z+∆z
Z
= ( vx |x+∆x − vx |x )dzdy
y z
x+∆x
Z z+∆z
Z
+ ( vy |y+∆y − vy |y )dzdx
x z
Z y+∆y
x+∆x Z
+ ( vz |z+∆z − vz |z )dydx
x y
y+∆y
Z z+∆z x+∆x
Z z+∆z y+∆y
Z x+∆x
1 1 1
Z Z Z
hvx i = vx dzdy; hvy i = vy dzdx; hvz i = vz dxdy
∆y∆z ∆x∆z ∆y∆x
y z x z y x
(1.68)
21
Notas de clase: Francisco J. Valdés Parada
Entonces, se tiene que
Z
n · v dA = hvx i|x+∆x − hvx i|x ∆y∆z + hvy i|y+∆y − hvy i|y ∆x∆z + hvz i|z+∆z − hvz i|z ∆y∆x
∆A
(1.69)
Dividiendo esta ecuación entre el volumen del cubo y tomando el lı́mite cuando ∆V → 0 se obtiene
que
1 ∂vx ∂vy ∂vz
Z
lı́m n · v dA = + + (1.70)
∆V →0 ∆V ∂x ∂y ∂z
∆A
Note que la cantidad del lado derecho de esta ecuación puede representarse en forma vectorial como
∂vx ∂vy ∂vz
+ + =∇·v (1.71)
∂x ∂y ∂z
Lo cual se conoce como la divergencia del vector v. Este resultado entonces muestra que la
divergencia de un vector corresponde al flujo neto de ese vector multiplicado por el área por unidad
de volumen. De esta forma, ∇ · v representa el flujo volumétrico neto de salida por unidad de
volumen. Cuando ∇ · v > 0, el flujo se dice que es divergente y cuando ∇ · v < 0 se dice que es
convergente y cuando ∇ · v = 0 el flujo es incompresible. Más adelante explicaremos acerca del flujo
compresible e incompresible en función de la ecuación de continuidad.
En general cuando una función (vectorial o tensorial) tiene divergencia igual a cero se dice que
es solenoidal. Una función solenoidal es aquella que en la que no hay generación, esto es, todo lo
que entre por dicha función en cada punto saldrá.
ρv · ∇v = −∇p + ρg + |∇{z
· τ}
| {z } | {z } |{z}
esfuerzos esfuerzos esfuerzos esfuerzos
inerciales superficiales volumétricos viscosos
22
Notas de clase: Francisco J. Valdés Parada
es
∂p
ρv · ∇v1 = − + ρg1 + ∇ · τi1 ei
∂x1
o bien, en forma expandida,
∂v1 ∂v1 ∂v1 ∂p ∂τ11 ∂τ21 ∂τ31
ρ v1 + v2 + v3 =− + ρg1 + + +
∂x1 ∂x2 ∂x3 ∂x1 ∂x1 ∂x2 ∂x3
Ejercicio 1.1.
Demuestre las siguientes identidades
∂2
2 ∂ ∂ ∂ ∂
∇ = ∇ · ∇ = ei · ej = = (1.75)
∂xi ∂xj ∂xi ∂xi ∂xi ∂xi
∂ 2 vk ∂ 2 v1 ∂ 2 v2 ∂ 2 v3
∇2 v = ek = e1 + e2 + e3 (1.77)
∂xi ∂xi ∂xi ∂xi ∂xi ∂xi ∂xi ∂xi
2
∂ v1 ∂ 2 v1 ∂ 2 v1 ∂ 2 v2 ∂ 2 v2 ∂ 2 v2 ∂ 2 v3 ∂ 2 v3 ∂ 2 v3
= + + e1 + + + e2 + + + e3
∂x21 ∂x22 ∂x23 ∂x21 ∂x22 ∂x23 ∂x21 ∂x22 ∂x23
∇2 ϕ = 0 (1.78)
23
Notas de clase: Francisco J. Valdés Parada
magnotud de este vector es la magnitud del eje de rotación. De hecho, en mecánica de fluidos se
define al vector de vorticidad como:
ω =∇×v (1.80)
Este vector indica la tendencia de una partı́cula de fluido a rotar como se verı́a por un observador
que viaja con dicha partı́cula.
Ejercicio 1.2.
a × (b × c) = b(a · c) − c(a · b)
24
Notas de clase: Francisco J. Valdés Parada
le llamaremos del continuo. La diferencia entre estos dos niveles de escala consiste en que en el nivel
atómico, los átomos y las moléculas están separados por espacios vacı́os, por lo que para describir,
por ejemplo, la temperatura de una masa de moléculas tendrı́amos que cuantificar la temperatura
de cada molécula en dicha masa. Por otra parte, en el nivel del continuo, las propiedades (como
temperatura, velocidad, concentración, etc.) son funciones continuas del espacio y el tiempo, esto
es, ya no identificamos los espacios vacı́os entre átomos y moléculas. Lo anterior sugiere que las
propiedades que se manejan en el nivel de escala del continuo son un tipo de promediado de la
información a escala molecular.
Entonces, en el nivel del continuo ya no hablamos de moléculas y vacı́o sino de fases, las cuales
a su vez están separadas entre sı́ mediante interfases. En este contexto, entendemos como fase
una porción de un material cuya estructura intermolecular no cambia, en una cierta extensión del
espacio que ocupa el material. Ası́, en un sistema homogéneo sólo hay una fase y en un sistema
heterogéneo hay más de una fase, separada mediante una o más interfases. Hablaremos con más
precisión acerca de las interfases más adelante.
Evidentemente, llevar a cabo estudios al nivel del continuo es menos realista que hacerlo al
nivel atómico o molecular. Esta pérdida de precisión muchas veces es aceptable en aplicaciones de
la ingenierı́a y se ve compensada además con la sencillez que implican los modelos matemáticos
a este nivel de escala. De manera general, puede establecerse que la hipótesis del continuo es
aplicable cuando la longitud caracterı́stica del instrumento de medición, `, que estemos utilizando
para estudiar (experimentalmente o teóricamente) un sistema sea mucho mayor2 que la trayectoria
libre media entre moléculas3 , `mol . Como lo indica Whitaker (1992a), para que nuestro instrumento
de medición no pierda demasiada información, su longitud caracterı́stica debe ser también mucho
menor que la longitud caracterı́stica del sistema a estudiar, L, esto es:
`mol ` L (1.81)
Una amplia variedad de sistemas de interés en ciencias básicas e ingenierı́a son de naturaleza
jerárquica; esto es, los procesos que tienen lugar en un cierto nivel de escala están determinados por
los procesos que se llevan a cabo en niveles de escala inferiores. En el estudio de sistemas jerárquicos,
es posible identificar un número discreto de niveles de escala que permiten definir diversos continuos.
Dichos niveles suelen ir desde el molecular, pasando por el nanoscópico (Slattery y col., 2004),
hasta llegar, por ejemplo, a la escala de campo o industrial (Whitaker, 2009). De esta forma, puede
razonarse que lo que ocurre en una interfase a nivel microscópico pueda ser determinante en el
desempeño de una planta quı́mica (ver Fig. 11 en Whitaker, 2009). En ocasiones, el pasar de un
nivel de escala a otro involucra promediar ecuaciones diferenciales parciales4 , que gobiernan el
transporte en un nivel de escala inferior (por ejemplo, del orden de micrómetros), para obtener
ecuaciones diferenciales ordinarias5 que describen el sistema en un nivel de escala superior (por
ejemplo, de metros). Con este mismo proceso es posible incluso pasar de una ecuación diferencial
parcial a una ecuación algebraica para describir un determinado fenómeno. De esta forma, no es
sorprendente que las ecuaciones que gobiernen el transporte de masa en una operación unitaria o en
un reactor sean considerablemente más sencillas que las ecuaciones que gobiernan el transporte en
un nivel de escala microscópico. La reducción del nivel de complejidad de los modelos matemáticos
2
Esto es, al menos uno o dos órdenes de magnitud mayor
3
La trayectoria libre media es la distancia promedio que puede recorrer una molécula respecto a las demás sin
incurrir en una colisión.
4
Aquellas en las que la variable dependiente es función de más de una variable independiente.
5
Aquellas en las que la variable dependiente es función de sólo una variable independiente.
25
Notas de clase: Francisco J. Valdés Parada
con el incremento en la escala no es exclusiva de aplicaciones ingenieriles, por ejemplo, la ecuación
de Navier-Stokes para fluidos Newtonianos, que es válida a niveles de escala microscópicos, es
considerablemente más simple que la ecuación de Schrödinger a nivel molecular.
De acuerdo a diversos autores (Bear y Cheng, 2010; Pinder y Gray, 2008; Bird y col., 2007;
Whitaker, 1999), los modelos matemáticos en sistemas macroscópicos pueden deducirse a partir
de sus contrapartes microscópicas. Entre las herramientas necesarias para llevar a cabo estas
deducciones se encuentra el promediado y el escalamiento. Wood (2009) aclara la diferencia entre
la operación matemática de promediar (es decir, aplicar un operador de integración) y escalar. Esta
última se entiende como una reducción sistemática del número de grados de libertad involucrados en
el modelo de transporte de la escala más pequeña. Ası́, por ejemplo, si se estudia el flujo vehicular,
se puede eliminar información como el modelo o color de cada vehı́culo para describir de manera
macroscópica el aforo que tiene una cierta avenida. Esto quiere decir que, mediante un proceso de
promediado y escalamiento, es posible definir nuevos continuos a niveles de escala superiores que
el correspondiente al primer continuo.
Ejercicio 1.3.
En la desigualdad (1.81) se acota la longitud caracterı́stica del instrumento de medición entre dos
longitudes caracterı́sticas: `mf y L. Sin embargo, es natural pensar que existan distintas trayectorias
libres medias en las moléculas de una sustancia y que el sistema donde se encuentre tenga diversas
longitudes caracterı́sticas. Explique cómo eligirı́a de una distribución de valores de `mf y L, los que
utilizarı́a en la desigualdad (1.81). Justifique detalladamente su respuesta.
Ejercicio 1.4.
Proponga tres ejemplos donde aplique una separación de escalas caracterı́sticas que involucre tres
niveles de escala.
26
Notas de clase: Francisco J. Valdés Parada
2. Derivada temporal Euleriana: Esta descripción comprende el estudio de la dinámica de
una propiedad desde un punto fijo. Note que bajo este enfoque no se tiene interés en la masa
que entre o salga del punto fijo donde se observa. A este tipo de derivada también se le conoce
como derivada parcial.
3. Derivada temporal Lagrangiana: En este caso se estudia la evolución de una propiedad
asociada a una partı́cula de fluido. Una partı́cula se entiende como un dominio que siempre
tiene la misma cantidad de materia y que se mueve con la velocidad del fluido. Este tipo de
derivada se conoce como derivada material.
Notemos que los tres tipos de derivadas expuestos arriba están asociados a dominios que pueden
ser arbitrarios, fijos o materiales. A continuación mostraremos las expresiones matemáticas de estos
tres tipos de derivadas temporales.
Para entender a la derivada total consideremos una propiedad f la cual es función del tiempo t
y de la posición en la cual se encuentra el observador, la cual se determina con el vector de posición
r. Esto es, f (t, r). Supongamos, por el momento, que se utiliza un sistema coordenado cartesiano.
Entonces, el vector de posición es:
r = xex + yey + zez (1.82)
donde (x, y, z) ubican una posición, la cual puede cambiar con el tiempo. De esta forma, al tomar
la derivada de f con respecto al tiempo, se obtiene usando la regla de la cadena que
df df dx df dy df dz df
= + + + (1.83)
dt dt r dt dx t,y,z dt dy t,x,z dt dz t,x,y
En una derivada parcial (o Euleriana), se deriva una función respecto a sólo una de las variables
independientes y se mantienen constantes el resto. Bajo esta perspectiva, la ecuación anterior puede
escribirse como:
df ∂f ∂f ∂f ∂f
= + wx + wy + wz (1.84)
dt ∂t ∂x ∂y ∂z
donde se definieron las componentes del vector de velocidad del observador como:
dx dy dz
wx = ; wy = ; wz = (1.85)
dt dt dt
O bien en forma más compacta
df ∂f
= + w · ∇f (1.86)
dt ∂t
Note que el resultado dado en la ecuación (1.84) está compuesto por dos partes: la primera que
corresponde a la derivada parcial de f con respecto del tiempo, esto es a la cuantificación de los
cambios de la propiedad respecto a un sistema fijo, y a las contribuciones debidas al movimiento
del observador. Evidentemente, cuando el observador no se mueve, la derivada total se reduce a
la derivada parcial o Euleriana. De hecho, este tipo de descripción es la más adecuada cuando se
estudian sistemas que no modifican su posición ni tamaño con el tiempo o bien cuando se tienen
instrumentos de medición que están fijos en el espacio.
Por último, la definición de derivada material puede recuperarse también a partir de la ecuación
(1.84) para el caso particular en el que el observador (o el elemento de medición) se mueva a a la
misma velocidad que el fluido. Esto es,
Df df ∂f ∂f ∂f ∂f
= = + vx + vy + vz (1.87)
Dt dt w=v
∂t ∂x ∂y ∂z
27
Notas de clase: Francisco J. Valdés Parada
O bien en forma más compacta
Df ∂f
= + v · ∇f (1.88)
Dt ∂t
Fı́sicamente, el primer término de esta ecuación se suele asociar con la acumulación de la
propiedad f y los tres términos restantes se asocian con el transporte convectivo de f . Note que
aún en estado estacionario la derivada material (y también la derivada total) de una propiedad con
respecto al tiempo no es cero, sino que considera únicamente el transporte convectivo.
Es importante recalcar que las definiciones anteriores (incluso de la derivada parcial) están
referidas a dominios y por lo tanto a niveles de escala. Esto es, una derivada material no es igual
si está referida a una velocidad microscópica que cuando se utiliza una velocidad macroscópica.
Más aún, de las definiciones anteriores es claro que la interpretación de una derivada como la razón
de cambio en el tiempo de una propiedad no es suficiente para saber a qué nos referimos pues
es necesario saber la velocidad del observador y la escala a la cual se hacen las mediciones. Por
supuesto las derivadas de propiedades pueden también tomarse respecto al espacio y por ello es
necesario dedicar los siguientes párrafos a este tema.
O bien, Z
∇ · v dV = n · v dA (1.90)
∂A
Esta ecuación es válida para cualquier diferencial de volumen que compongan a un volumen mayor.
De hecho, se necesita un número infinito de diferenciales de volumen para reconstruir el volumen
mayor. Cuando se suman todos estos diferenciales de volumen, se obtiene entonces
∞
X X∞ Z
∇ · v dVk = n · v dA (1.91)
k=1 k=1∂A
k
Note que las integrales de las caras internas de los diferenciales de superficie se cancelan entre sı́ y
el lado derecho de la ecuación anterior toma la forma
X∞ Z X∞ Z Z
n · v dA = n · v dA ≡ n · v dA (1.92)
k=1∂A k=1∂A A
k ext
Por otro lado, el lado izquierdo de la ecuación (1.91) corresponde a la definición de una integral de
Rieman, esto es
X∞ Z
∇ · v dVk = ∇ · v dV (1.93)
k=1 V
28
Notas de clase: Francisco J. Valdés Parada
De esta forma, se obtiene la fórmula del teorema de la divergencia:
Z Z
∇ · v dV = n · v dA (1.94)
V A
Obviamente, esta fórmula no está limitada para el vector de velocidad, sino que funciona para
cualquier vector continuo f. De hecho, es de particular importancia considerar el caso en el que
f = cψ, donde c es un vector constante distinto de cero y ψ es un escalar que es función de la
posición. En este caso el teorema es:
Z Z
∇ · (cψ) dV = n · cψ dA (1.95)
V A
En la última igualdad se tomó en cuenta que c es independiente de la posición. Más aún, dado que
n · c = c · n, podemos escribir entonces
Z Z
c · ∇ψ dV = c · nψ dA (1.97)
V A
Dado que c no es cero, para que esta igualdad se satisfaga, se debe cumplir que
Z Z
∇ψ dV = nψ dA (1.99)
V A
La cual es la forma del teorema de la divergencia que aplica a un escalar que se conoce también
como la forma del gradiente del teorema de la divergencia.
Este teorema también puede aplicarse para el operador rotacional:
Z Z
∇ × a dV = n × a dA (1.100)
V A
29
Notas de clase: Francisco J. Valdés Parada
Ahora considere el caso en que dicho vector sea φ∇ϕ, el resultado es
Z Z
2
φ∇ ϕ + ∇φ · ∇ϕ dV = n · φ∇ϕ dA (1.102)
V A
Dado que ∇ϕ · ∇φ = ∇φ · ∇ϕ, es posible restar la ecuación anterior a la ecuación (1.102) para
obtener
Z Z
φ∇2 ϕ − ϕ∇2 φ = n · (φ∇ϕ − ϕ∇φ) dA
(1.104)
V A
Para poder avanzar es necesario conocer cómo cambia el diferencial de volumen en el tiempo. La
respuesta a esta pregunta la brinda la componente normal de la velocidad de desplazamiento en
cada punto de la superficie, esto es n · w. Si se integran en las superficies del elemento diferencial
arbitrario todas las contribuciones de n · w, se tiene entonces que
d∆V
Z
= n · w dA (1.109)
dt
∆A
30
Notas de clase: Francisco J. Valdés Parada
Obviamente ∆A es función del tiempo. Sustituyendo este resultado en la ecuación (1.108) y
tomando el lı́mite cuando ∆V → 0, resulta que
1 d dhf i 1
Z Z
lı́m f (r, t) dV = + hf i lı́m n · w dA (1.110)
∆V →0 ∆V dt dt ∆V →0 ∆V
∆V ∆A
1
Z
lı́m n · w dA = ∇ · w (1.111)
∆V →0 ∆V
∆A
1 d df
Z
f (r, t) dV = + f∇ · w (1.112)
dV dt dt
dV
Aquı́ se usó el lı́mite cuando ∆V → 0, note que bajo esta condición se cumple que hf i → f . El lado
izquierdo de esta ecuación puede reformularse al usar la definición de la derivada total como:
df ∂f ∂f
+ f∇ · w = + w · ∇f + f ∇ · w = + ∇ · (f w) (1.113)
dt ∂t ∂t
Entonces, la ecuación (1.112) toma la forma
d ∂f
Z
f (r, t) dV = + ∇ · (f w) dV (1.114)
dt ∂t
dV
Como en el teorema anterior, el paso siguiente es sumar un número infinito de elementos diferenciales
de volumen
∞ ∞
d ∂f
X Z X
f (r, t) dV = + ∇ · (f w) dVk (1.115)
dt ∂t
k=1 dVk k=1
Finalmente, usando el teorema de la divergencia en el último término del lado derecho, resulta el
teorema general del transporte:
d ∂f
Z Z Z
f (r, t) dV = dV + n · (f w) dA (1.117)
dt ∂t
V V A
Note que el resultado de tomar el lı́mite cuando ∆V → 0 en la ecuación (1.109) es una forma
particular de este teorema para el caso en que f = 1.
31
Notas de clase: Francisco J. Valdés Parada
Para el caso en el que el dominio de integración sea una región material, la derivada del lado
izquierdo es una derivada material y el teorema se conoce como el teorema del transporte de
Reynolds, el cual se expresa como
D ∂f
Z Z Z
f (r, t) dV = dV + n · (f v) dA (1.118)
Dt ∂t
V V A
Note que en este caso, el movimiento del volumen (y de sus superficies) está determinado por la
velocidad del fluido.
Por último, en el caso en el que el volumen es fijo, el teorema general del transporte se reduce a
∂ ∂f
Z Z
f (r, t) dV = dV (1.119)
∂t ∂t
V V
Ejercicio 1.5.
Utilice el teorema de la divergencia aplicado a un vector a en un dominio material:
Z Z
∇ · a dV = a · n dA (1.120)
Vm (t) Am (t)
junto con la definición de la derivada material dada en la Ec. (1.88) para mostrar que el teorema
del transporte de Reynolds puede expresarse como:
Z
D DA
Z
AdV = + A∇ · v dV (1.121)
Dt Dt
Vm (t) Vm (t)
Ejercicio 1.6.
Demuestre que si A = 1, el teorema del transporte de Reynolds indica que
dVm
Z
= ∇ · vdV
dt
Vm (t)
R
donde Vm = dV es el volumen del dominio material.
Vm (t)
32
Notas de clase: Francisco J. Valdés Parada
Para expresar la Ec. (1.122) en lenguaje matemático, se empieza por definir la masa, m ,
contenida en la región arbitraria Z
m= dm (1.123)
Va
dm = ρ dV (1.124)
Ahora bien, la tasa neta de entrada de masa por las fronteras de la región, se refiere a los kilogramos
por segundo de masa que pasan a través de las fronteras de la región. Para estimar esta cantidad
es necesario tomar en cuenta que la masa se mueve con una velocidad v, mientras que las fronteras
de la región se mueven con una velocidad w. Entonces la velocidad se debe cuantificar relativa a
la velocidad de las fronteras, esto es v − w. Si esta velocidad relativa se multiplica por la densidad
de la fase, el resultado es ρ(v − w) y representa la cantidad de masa por unidad de área y tiempo
que atraviesa las fronteras del dominio y es un vector que no necesariamente está alineado con las
fronteras del sistema. Para atender este problema, es necesario realizar el producto punto con el
vector normal a la superficie de la región arbitraria, esto es, n · ρ(v − w). Sin embargo, el resultado
ahora es la proyección hacia afuera de la región, como el interés está en la tasa neta de entrada a
la región arbitraria, es necesario hacer un cambio de signo e integrar en la superficie de la región.
El resultado es entonces
( )
Tasa neta de entrada de masa
Z
= − n · ρ(v − w) dA (1.127)
por las fronteras de la región
Aa
33
Notas de clase: Francisco J. Valdés Parada
En este momento, cabe preguntarse, ¿bajo qué condiciones es el resultado de una integral cero?,
algunas posibilidades son: 1) cuando el integrando es cero, 2) cuando los lı́mites de integración son
iguales, 3) cuando se integra en un punto. Dado que la región de integración no es un punto y sus
lı́mites no son iguales, se deduce que el integrando debe ser cero; por lo tanto, se puede extraer la
siguiente ecuación diferencial parcial
∂ρ
+ ∇ · (ρv) = 0 (1.130)
∂t
La cual se conoce como ecuación de continuidad y es una forma matemática de expresar que la
masa de un cuerpo se conserva. El vector de velocidad v en la Ec. (1.130) puede determinarse a
partir de la ecuación de conservación de cantidad de movimiento como se verá más adelante.
Vale la pena notar que, para casos en los que la densidad no cambie en la posición ni en el
tiempo6 , la Ec. (1.130) se reduce a
∇·v =0 (1.131)
Para concluir este apartado, es conveniente listar los pasos que se siguieron para obtener la ecuación
de continuidad:
1. Se escribió el principio de conservación en una ecuación con palabras.
3. Se buscó obtener una ecuación que estuviera expresada en la forma de una integral de volumen
igualada a cero.
4. Al concluir que el integrando de la ecuación del punto anterior debe ser cero, se obtuvo la
ecuación buscada.
Este procedimiento se utilizará más adelante para deducir la ecuación de conservación de cantidad
de movimiento.
Para concluir esta sección, es conveniente escribir la ecuación de continuidad en los tres sistemas
coordenados ortogonales más comunes:
Coordenadas cartesianas:
∂ρ ∂(ρv1 ) ∂(ρv2 ) ∂(ρv3 )
+ + + =0 (1.132)
∂t ∂x1 ∂x2 ∂x3
Coordenadas cilı́ndricas:
∂ρ 1 ∂(ρrvr ) 1 ∂(ρvθ ) ∂(ρvz )
+ + + =0 (1.133)
∂t r ∂r r ∂θ ∂z
Coordenadas esféricas:
∂ρ 1 ∂(ρr2 vr ) 1 ∂(ρvθ senθ) 1 ∂(ρvϕ )
+ 2 + + =0 (1.134)
∂t r ∂r rsenθ ∂θ rsenθ ∂ϕ
Ejercicio 1.7.
Deduzca la ecuación de continuidad partiendo de una región material y de una región fija. El
resultado final debe ser el mismo.
6
Como es el caso de algunos fluidos y sólidos a ciertas condiciones de temperatura y presión.
34
Notas de clase: Francisco J. Valdés Parada
Ejercicio 1.8.
Utilice la definición de la derivada material para demostrar que la ecuación de continuidad puede
escribirse como:
Dρ
+ ρ(∇ · v) = 0 (1.135)
Dt
Ejercicio 1.9.
Utilice este resultado junto con el dado en la Ec. (1.121) para obtener la forma especial del teorema
del transporte de Reynolds:
d DA
Z Z
ρAdV = ρ dV (1.136)
dt Dt
Vm (t) Vm (t)
fase-γ d i v isoria Aβ
γ
c i e
Superfi
βγ
n
I
@
@
@
R
@
γβ
n
fase-β
Figura 1.4: Esquema de una superficie divisoria Aβγ que separa a las fases homogéneas β y γ
incluyendo los vectores unitarios normales a Aβγ , nβγ = −nγβ .
35
Notas de clase: Francisco J. Valdés Parada
1.15.1. Condición de salto de masa en una interfase
Consideremos dos fases, β y γ separadas por una superficie divisoria Aβγ como se muestra en la
Figura 1.4. Para estudiar el sistema definimos un volumen arbitrario V = Vβ + Vγ . En este volumen
se cumple el principio de conservación de masa, el cual lo expresamos en forma integral como:
d
Z Z
ρ dV = − n · ρ (v − w) dA (1.137)
dt
V A
Donde ρ es la densidad en cualquier punto de V y no está restringida a una fase en particular. Para
entender la relación entre ρ, con las densidades en cada fase, ρβ y ργ , se define lo siguiente:
Z Z Z Z
ρ dV = ρβ dV + ργ dV + ρs dA (1.138)
V Vβ Vγ Aβγ
d
Z Z Z
ργ dV + ργ (vγ − w) · nγ dA + ργ (vγ − w) · nγβ dA = 0 (1.141a)
dt
Vγ Aγ Aβγ
d
Z Z Z
ρβ dV + ρβ (vβ − w) · nβ dA + ρβ (vβ − w) · nβγ dA = 0 (1.141b)
dt
Vβ Aβ Aβγ
La superficie divisoria puede desplazarse con una velocidad w en la dirección indicada por el vector
unitario nβγ = −nγβ como lo indica el producto w · nβγ . En la ecuación anterior nβ y nγ son
los vectores unitarios perpendiculares a las superficies Aβ y Aγ , respectivamente. Sumando las
ecuaciones anteriores resulta:
d d
Z Z Z Z
ργ dV + ρβ dV + ργ (vγ − w) · nγ dA + ρβ (vβ − w) · nβ dA
dt dt
Vγ Vβ Aγ Aβ
Z
=− [ρβ (vβ − w) − ργ (vγ − w)] · nβγ dA (1.142)
Aβγ
36
Notas de clase: Francisco J. Valdés Parada
Restando la ecuación anterior a la Ec. (1.140) resulta
d
Z Z
ρs dA = [ρβ (vβ − w) − ργ (vγ − w)] · nβγ dA (1.143)
dt
Aβγ Aβγ
Deseamos expresar esta ecuación como una sola integral de superficie, para lo cual requerimos
intercambiar diferenciación temporal e integración espacial. Dado que la integral es superficial en
lugar de volumétrica no podemos usar el teorema general del transporte. Para hacer frente a esta
situación, recurrimos al teorema del transporte superficial (ver Sección 3 en el trabajo de Whitaker,
1992b):
Z
d ∂ψs
Z
ψs dA = + ∇s · (ψs vs ) + ψs (∇s · nβγ )(w · nβγ ) dA (1.144)
dt ∂t
Aβγ Aβγ
donde:
∇s = ∇ · P; P = I − nβγ nβγ (1.145)
es el operador gradiente superficial y P es el tensor de proyección tangencial a la superficie divisoria.
Cabe mencionar que ∇s · nβγ = 2H, donde H es la curvatura media de la superficie divisoria
(Slattery y col., 2006).
Utilizando el teorema del transporte superficial, la Ec. (1.143) puede escribirse como:
Z
∂ρs
+ ∇s · (ρs vs ) + 2Hρs (w · nβγ ) − [ρβ (vβ − w) − ργ (vγ − w)] · nβγ dA = 0 (1.146)
∂t
Aβγ
Para situaciones en las cuales los términos del lado izquierdo de la ecuación anterior sean
despreciables, el resultado se reduce a la condición de continuidad del flujo de masa a través de una
superficie divisoria:
ργ (vγ − w) · nβγ = ρβ (vβ − w) · nβγ (1.148)
Ejercicio 1.10.
Demuestre que, si no hay transferencia de masa, la condición de salto dada por ργ (vγ − w) · nβγ =
ρβ (vβ − w) · nβγ se reduce a la continuidad de las velocidades de cada fase:
Ejercicio 1.11.
En muchos casos de estudio es aceptable imponer la condición de no deslizamiento, la cual implica
que:
vβ · tβγ = vγ · tβγ (1.150)
37
Notas de clase: Francisco J. Valdés Parada
donde tβγ es el vector unitario tangente a la superficie divisoria. Suponiendo que la fase β es un
sólido fijo e impermeable a la transferencia de masa, explique bajo qué condiciones puede deducirse
que
vγ = 0 en Aβγ
Note que la ecuación anterior aplica al vector velocidad y no sólo a su componente normal o
tangencial.
Ejercicio 1.12.
En flujos de gases la condición de no deslizamiento deja de ser válida y Navier propuso la siguiente
expresión para calcular la velocidad de deslizamiento:
λγ
tβγ · vγ = ξ nβγ · τ · tβγ (1.151)
µγ
Ejercicio 1.13.
Demuestre que la curvatura de un cı́rculo es 2H = 1/r. Para ello, parta de la definición de curvatura;
2H = [∇ · (I − nβσ nβσ )] · nβγ = y considere que nx = cos θ, ny = senθ, θ = atan(y/x) donde
x = r cos θ y y = rsenθ.
dr
=v (1.153)
dt
Esta ecuación diferencial puede resolverse sujeta a la siguiente condición inicial:
t = 0, r = r0 (1.154)
38
Notas de clase: Francisco J. Valdés Parada
La solución puede escribirse como
tZ0 =t
Ejercicio 1.14.
Considere una situación en la que las componentes del vector velocidad en coordenadas cartesianas
son (Slattery, 1999):
r1 r2
v1 = ; v2 = ; v3 = 0
1+t 1 + 2t
Demuestre que en el plano r3 , las lı́neas de camino son
1/2
r2 r1
= 2 −1
r20 zr0
39
Notas de clase: Francisco J. Valdés Parada
Ejercicio 1.15.
Demuestre que, para la distribución de velocidad dada en el ejercicio anterior, las lı́neas de corriente
son: (1+t)/(1+2t)
r2 r1
=
r20 r10
Ejercicio 1.16.
El concepto de la función de corriente ψ es especialmente útil para estudiar flujos bidimensionales.
En estos casos, el vector velocidad se expresa como
e3
v = ∇ψ × (1.161)
h3
donde h3 es un factor de escala. En coordenadas cartesianas h3 = 1, en coordenadas cilı́ndricas
axisimétricas, e3 = eθ y h3 = r mientras que en coordenadas esféricas e3 = eϕ y h3 = rsenθ.
Demuestre que la Ec. (1.161) satisface la ecuación de continuidad para flujo bidmiensional e
incompresible.
Ejercicio 1.17.
Utilice la Ec. (1.161) para demostrar que la función de corriente ψ es constante a lo largo de una
lı́nea de corriente. Esto es:
v · ∇ψ = 0 (1.162)
r = r(R(r0 , τ ), t) (1.164)
Ejercicio 1.18.
Muestre que, para el perfil de velocidad considerado en los ejercicios anteriores, la lı́nea de marcha
está especificada por
1+t
r1 = r10
1+τ
1 + 2t 1/2
r2 = r20
1 + 2τ
r3 = r30
40
Notas de clase: Francisco J. Valdés Parada
Todas estas representaciones son equivalentes para flujo en estado estacionario y solo pueden
distinguirse en estado transitorio. Más aún, todas estas lı́neas pueden observarse experimentalmente
y sólo las lı́neas de corriente tienen una aplicación mayoritariamente teórica. Por último, cabe
mencionar que las lı́neas de corriente pueden graficarse para cualquier campo vectorial, no solo
para la velocidad, sino que también están definidas para fluxes (de calor y masa) o incluso para el
vector de vorticidad.
Ejercicio 1.19.
Si para el flujo laminar las lı́neas de corriente para la velocidad se observan como lı́neas horizontales,
¿ cómo se visualizan las lı́neas de corriente para el vector de vorticidad para el mismo tipo de flujo
?
41
Notas de clase: Francisco J. Valdés Parada
42
Notas de clase: Francisco J. Valdés Parada
Contenido
2.1. Leyes de Newton y de Euler . . . . . . . . . . . . . . . . . . . . . . . 43
2.2. El lema de Cauchy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
2.3. Teorema fundamental de Cauchy . . . . . . . . . . . . . . . . . . . . 48
2.4. Primera ecuación de Cauchy . . . . . . . . . . . . . . . . . . . . . . . 51
2.5. Segunda ecuación de Cauchy: Simetrı́a del tensor de esfuerzos . . 52
2.6. Condición de salto de cantidad de movimiento . . . . . . . . . . . . 53
2.7. Esfuerzos viscosos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.7.1. Ley de Newton de la viscosidad . . . . . . . . . . . . . . . . . . . . . . 57
2.7.2. Fluidos Newtonianos generalizados . . . . . . . . . . . . . . . . . . . . 59
2.8. Tipos de flujo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
2.8.1. Flujo reptante . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.8.2. Flujo potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.9. La ecuación de energı́a mecánica . . . . . . . . . . . . . . . . . . . . 63
43
Notas de clase: Francisco J. Valdés Parada
donde las definiciones de la masa, velocidad y fuerzas no son siempre precisas en los libros de texto
de fı́sica básica. Lo que sı́ queda claro es que el vector f es la fuerza que se aplica sobre el cuerpo.
La tercera ley de Newton puede escribirse matemáticamente como:
donde f12 es la fuerza que ejerce el cuerpo 1 sobre el cuerpo 2. A pesar de la extensa popularidad
que gozan las leyes de Newton, no son las más indicadas para estudiar el movimiento y deformación
de los fluidos. Para ello, utilizaremos las dos leyes de la mecánica del continuo que se atribuyen a
Euler:
Dirijamos por el momento la atención a la primera ley de Euler, la cual la podemos expresar
en forma de ecuación con palabras como: la siguiente ecuación con palabras:
( ) ( )
La razón de cambio en el tiempo de la Las fuerzas que actúan
= (2.3)
cantidad de movimiento lineal de un cuerpo sobre el cuerpo
Claramente esta expresión aplica a una región material, por lo que es más conveniente enunciar el
principio de conservación para un dominio arbitrario, el cual se puede expresar en palabras como
La razón de cambio en el tiempo de
Flujo neta de entrada de
la cantidad de movimiento en = cantidad de movimiento por
un dominio arbitrario las fronteras del dominio
( ) ( )
Fuerzas que actúan Fuerzas que actúan en
+ + (2.4)
dentro del dominio las fronteras del dominio
Para traducir esta ecuación al lenguaje matemático, es conveniente recordar lo que se aprendió
en la deducción de la ecuación de continuidad. Recordando que la cantidad de movimiento es el
producto de la masa por la velocidad, tenemos entonces que
La razón de cambio en el tiempo de
d
Z
la cantidad de movimiento en = ρv dV (2.5)
dt
un dominio arbitrario
VA
Mientras que el flujo neto de entrada por las fronteras está dada por
Flujo neto de entrada de
Z
cantidad de movimiento por = − n · ρv(v − w) dA (2.6)
las fronteras del dominio
AA
44
Notas de clase: Francisco J. Valdés Parada
Ahora bien, las fuerzas que actúan dentro del dominio por unidad de masa se denotan por b y
corresponden a la gravedad o bien a campos electromagnéticos. Entonces,
( ) Z
Fuerzas que actúan
= ρb dV (2.7)
dentro del dominio
VA
Por último, se utiliza el vector de esfuerzos tn para representar las fuerzas superficiales por unidad
de área. Note que este vector requiere conocer sobre qué normal está siendo aplicado. Tenemos
entonces que ( ) Z
Fuerzas que actúan en
= tn dA (2.8)
las fronteras del dominio
AA
De esta forma, la ecuación de conservación de cantidad de movimiento puede expresarse como
d
Z Z Z Z
ρv dV + n · ρv(v − w) dA = ρb dV + tn dA (2.9)
dt
VA AA VA AA
Para una región material, la primera ley de Euler puede expresarse como:
D
Z Z Z
ρv dV = ρb dV + t(n) dA (2.10)
Dt
Vm Vm Am
Por último, para una región fija, este principio de conservación se expresa como
d
Z Z Z Z
ρv dV + n · ρvv dA = ρb dV + tn dA (2.11)
dt
Vf Af Vf Af
Antes de continuar es conveniente encontrar la forma en la que estas tres formulaciones dan
lugar a la misma expresión. Para ello, se utiliza el teorema general del transporte en la ecuación
(2.9) para obtener
∂(ρv)
Z Z Z Z
dV + n · ρvv dA = ρb dV + tn dA (2.12)
∂t
VA AA VA AA
Note que esta misma ecuación resulta de aplicar el teorema del transporte de Reynolds a la ecuación
(2.10). De la misma forma, dado que en la ecuación (2.11) la región de integración no cambia con
el tiempo, es posible introducir el operador de derivación parcial en el integrando y recuperar la
ecuación de arriba. Más aún, usando el teorema de la divergencia en el segundo término del lado
izquierdo de la ecuación anterior se tiene que
Z
∂(ρv)
Z Z
+ ∇ · (ρvv) dV = ρb dV + tn dA (2.13)
∂t
VA VA AA
Al desarrollar las derivadas temporales y espaciales del integrando del lado izquierdo de esta
ecuación se obtiene que
∂(ρv) ∂v ∂ρ Dv
+ ∇ · (ρvv) = ρ + v · ∇v + + ∇ · ρv v = ρ (2.14)
∂t ∂t ∂t Dt
| {z } | {z }
Dv/Dt =0
45
Notas de clase: Francisco J. Valdés Parada
De esta forma, la ecuación (2.13) puede expresarse como
Z
Dv
Z
ρ − ρb dV = tn dA (2.15)
Dt
VA AA
Ejercicio 2.1.
Comparando la Ec. (2.10) con la Ec. (2.1), podemos deducir que la velocidad a la que se refieren
las dos primeras leyes de Newton es la velocidad del centro de masa del cuerpo, es decir:
1
Z
vcm = ρv dV (2.17)
m
Vm
Esto quiere decir que la primera ley de Euler contiene a las primeras dos leyes de Newton. Para
mostrar la relación entre la tercera ley de Newton y la primera ley de Euler consideremos dos
cuerpos que se encuentran contenidos en un tercer cuerpo que está aislado en el sentido de que no
se ejercen fuerzas externas sobre éste. Bajo estas condiciones tenemos para los cuerpos 1 y 2:
d
Z
ρ1 v1 dV = f12 (2.18a)
dt
Vm1
d
Z
ρ2 v2 dV = f21 (2.18b)
dt
Vm2
Dado que: Z Z Z
ρvdV = ρ1 v1 dV + ρ2 v2 dV (2.19)
Vm Vm1 Vm2
puede concluirse que f12 = −f21 al sumar las ecs. (2.18) demostrando ası́ la tercera ley de Newton.
46
Notas de clase: Francisco J. Valdés Parada
n1
∆L
A
A
A
A
U
nlat
n2
Figura 2.1: Región arbitraria de espesor ∆L, las superficies superior e inferior son las mismas e
iguales a An y la superficie lateral es Alat .
superficies. Considere el cuerpo dibujado en la Figura 2.1, aplicando la primera ley de Euler en la
forma dada en la Ec. (2.15), tenemos
Z
Dv
Z Z Z
ρ − ρb dV = t(n1 ) dA + t(n2 ) dA + t(nlat ) dA (2.20)
Dt
VA An An Alat
Notemos que el volumen del cuerpo es VA = An ∆L. Para avanzar en nuestros desarrollos es
conveniente definir los siguientes operadores de promediado:
1
Z
hψi = ψdV (2.21a)
An ∆L
VA
1
Z
hψilat = ψdA (2.21b)
Alat
Alat
47
2.2 Cauchy’s fundamental theorem
We have learned something very important about the stress vector by applying Euler’s first law to the
slab illustrated in Figure 2-2, and we can learn more with an application to the tetrahedron shown in
Figure 2-3. Once again, you should think of the tetrahedron as a body cut out of a distinct body according
Notas de clase: Francisco J. Valdés Parada
Figure
Figura 2.2: Stresses
2-3. de
Vectores acting
esfuerzos on a sobre
actuando tetrahedron
un tetraedro.
Dado que los lı́mites de integración son arbitrarios, se extrae el siguiente resultado:
2
t(n1 ) = −t(n2 ) (2.24)
James, G. and James, R.C. 1959, Mathematics Dictionary (multilingual edition), “lemma, a theorem proved for the
use in theComo
proofnof1 =
another
−n2 =theorem”, D. van
n, podemos Nostrand
reescribir esteCo., Inc., New
resultado York.
como
El cual se conoce como el Lema de Cauchy y enuncia que los vectores de esfuerzos que actúan sobre
caras opuestas de la misma superficie en un punto dado son iguales en magnitud y opuestos en
dirección. Debe ser claro que este resultado no está restringido a geometrı́as como la usada en la
Figura 2.1.
48
Notas de clase: Francisco J. Valdés Parada
Aplicando la primera ley de Euler, en la forma dada por la Ec. (2.15) resulta:
Z
Dv
Z Z Z Z
ρ − ρb dV = t(n) dA + t−(i) dA + t−(j) dA + t−(k) dA (2.26)
Dt
∆V ∆A ∆Ax ∆Ay ∆Az
1
Z
hψi = ψ dV (2.27a)
∆V
∆V
1
Z
hψi∆A = ψ dA (2.27b)
∆A
∆A
1
Z
hψi∆Ax = ψ dA (2.27c)
∆Ax
∆Ax
1
Z
hψi∆Ay = ψ dA (2.27d)
∆Ay
∆Ay
1
Z
hψi∆Az = ψ dA (2.27e)
∆Az
∆Az
Más aún, del teorema del área proyectada (ver prueba en el Apéndice A de Whitaker, 1992a),
tenemos que
49
t (j) i Tyx j Tyy k Tyz , force per unit area acting on the y-surface (2-25b)
t (k ) i Tzx j Tzy k Tzz , force per unit area acting on the x-surface (2-25c)
Here the subscripts associated with Txx , Tzy , Tyx , etc., have a very definite meaning. In each case, the first
subscript represents the surface upon which the force is acting while the second subscript indicates the
direction in which the force Notas
is acting.
deThis ideaFrancisco
clase: is illustrated
J. inValdés
Figure Parada
2-4 where we have shown the
Figure
Figura 2.3: 2-4. Positive
Esquema de lascomponents of thedel
componentes stress vector
tensor deacting on theactuando
esfuerzos positive y-surface
sobre un cubo.
components of t (j) , all of which are positive.
El cual es el Teorema fundamental de Cauchy y es de especial importancia para nosotros porque nos
indica
2.3 la dependencia
Stress que tiene el vector de esfuerzos t(n) con el vector unitario normal n (además
in static fluids
del espacio
The y el tiempo).
commonly El resultado
accepted definition anterior
of a fluid puede
is that itrearreglarse como:
will deform continuously under the action of a
shear stress. This means that the stress on any surface must be normal to that surface when the fluid is in a
(n) = n · explicitly
static condition, and we identify thistcondition it(i) + jtas(j) + kt(k) (2.33)
Para hacer más compacto este resultado, se define el tensor de esfuerzos T como:
iiTxx + ijTxy + ikTxz
T = +jiTyx + jjTyy + jkTyz (2.36)
+kiTzx + kjTzy + kkTzz
t(n) = n · T (2.37)
Para tener una idea más clara del significado fı́sico de las componentes del tensor de esfuerzos,
consideremos el esquema mostrado en la Figura 2.3. Como puede notarse el primer subı́ndice de
50
Notas de clase: Francisco J. Valdés Parada
cada componente indica la superficie sobre la cual se aplica el esfuerzo y el segundo subı́ndice nos
dice la dirección a lo largo de la cual se aplica el esfuerzo. Ası́, Tyz es el esfuerzo que se aplica sobre la
superficie y en la dirección z. Notemos que las componentes diagonales del tensor de esfuerzos actúan
en la dirección perpendicular a cada superficie. Es decir, corresponden a los esfuerzos normales,
mientras que las componentes fuera de la diagonal actúan en la dirección tangencial a la superficie
y por ello se conocen como esfuerzos de corte o viscosos que ocasionan la deformación del fluido.
Para un fluido estático no hay esfuerzos viscosos y la única fuerza que se aplica sobre el cuerpo es
la presión en la dirección normal a la superficie del cuerpo. Esto es:
Por lo que, tras utilizar el teorema de la divergencia, la Ec. (2.39) puede escribirse como una sola
integral volumétrica: Z
Dv
ρ − ρb − ∇ · T dV = 0 (2.40)
Dt
VA
Ya que el integrando es una función continua y los lı́mites de integración son arbitrarios, extraemos
directamente la primera ecuación de Cauchy:
Dv
ρ = ρb + ∇ · T (2.41)
Dt
Es importante hacer notar que este resultado es válido para cualquier medio continuo, ya sea sólido,
lı́quido o gas. Para el caso particular de fluidos estáticos en los que la única fuerza volumétrica sea
la gravedad, tenemos que b = g, v = 0 y T = −pI y la ecuación anterior se reduce a la ecuación
fundamental de la hidrostática
Ejercicio 2.2.
Utilice la ecuación de continuidad y la definición de derivada material para demostrar que la primera
ecuación de Cauchy puede escribirse como sigue:
∂(ρv)
+ ∇ · (ρvv) = ρb + ∇ · T (2.43)
∂t
51
Notas de clase: Francisco J. Valdés Parada
2.5. Segunda ecuación de Cauchy: Simetrı́a del tensor de esfuerzos
Para probar la simetrı́a del tensor de esfuerzos, recurrimos a la segunda ley de Euler presentada
anteriormente como:
Dr × v
Z Z Z
ρ dV = r × ρb dV + r × t(n) dA (2.44)
Dt
VA VA AA
Dirijamos la atención ahora al término en el lado derecho, el cual puede escribirse, con ayuda del
teorema de la divergencia como:
Z Z Z
r × TT · n dA = ∇ · r × TT dV
r × n · T dA = (2.48)
AA AA VA
Sin embargo, debido a la primera ecuación de Cauchy (Ec. 2.41), el resultado anterior se reduce a:
Z
εijk Tij ek dV = 0 (2.51)
VA
Tomando en cuenta las propiedades del sı́mbolo de permutación, no es difı́cil concluir de esta
expresión que:
T12 = T21 ; T13 = T31 ; T23 = T32 (2.52)
Es decir, las componentes no diagonales del tensor total de esfuerzos son iguales y por lo tanto
T = TT (2.53)
En otras palabras, el tensor de esfuerzos es un tensor simétrico y es una restricción que impone la
segunda ley de Euler acerca de las fuerzas que actúan sobre un cuerpo. La Ec. (2.53) es la segunda
ecuación de Cauchy.
52
Notas de clase: Francisco J. Valdés Parada
Ejercicio 2.3.
Utilice notación indicial para demostrar la igualdad dada en la Ec. (2.50).
Al igual que en el caso de la condición de salto de masa, es conveniente relacionar las propiedades
que son válidas en todo el dominio con aquellas que están definidas en las fases y en la superficie
divisoria mediante:
Z Z Z Z
ρv dV = ρβ vβ dV + ργ vγ dV + ρs vs dA (2.55a)
V Vβ Vγ Aβγ
Z Z Z Z
ρb dV = ρβ bβ dV + ργ bγ dV + ρs bs dA (2.55b)
V Vβ Vγ Aβγ
Z Z Z I
n · T dA = nβ · Tβ dA + nγ · Tγ dA + ns · Ts dσ (2.55c)
A Aβ Aγ C
d d d
Z Z Z
ρβ vβ dV + ργ vγ dV + ρs vs dA
dt dt dt
Vβ Vγ Aβγ
Z Z
+ nβ · ρβ vβ (vβ − wβ ) dA + nγ · ργ vγ (vγ − wγ ) dA
Aβ Aγ
Z Z Z
= ρβ bβ dV + ργ bγ dV + ρs bs dA
Vβ Vγ Aβγ
Z Z Z
+ nβ · Tβ dA + nγ · Tγ dA + (∇s · Ts + 2Hnβγ · Ts ) dA (2.57)
Aβ Aγ Aβγ
53
Notas de clase: Francisco J. Valdés Parada
Note que se utilizó el teorema de la divergencia superficial en el último término (ver Ec. (2.21) en
Ochoa-Tapia y col., 1993):
I Z
ns · B dσ = (∇s · B + 2Hnβγ · B) dA (2.58)
C Aβγ
Dirigiendo ahora la atención a las fases, podemos escribir la forma macroscópica de la primera
ecuación de Cauchy para cada fase como:
d
Z Z Z
ρβ vβ dV + nβ · [ρβ vβ (vβ − wβ )] dA + nβγ · [ρβ vβ (vβ − wβγ )] dA
dt
Vβ Aβ Aβγ
Z Z Z
= ρβ bβ dV + nβ · Tβ dA + nβγ · Tβ dA (2.59a)
Vβ Aβ Aβγ
d
Z Z Z
ργ vγ dV + nγ · [ργ vγ (vγ − wγ )] dA + nγβ · [ργ vγ (vγ − wβγ )] dA
dt
Vγ Aγ Aβγ
Z Z Z
= ργ bγ dV + nγ · Tγ dA + nγβ · Tγ dA (2.59b)
Vγ Aγ Aβγ
tenemos entonces:
Z
∂ρs vs
+ ∇s · (ρs vs vs ) + ρs vs (∇s · nβγ )(w · nβγ ) dA
∂t
Aβγ
Z
− nβγ · [ρβ vβ (vβ − w) − ργ vγ (vγ − w)] dA
Aβγ
Z
= [ρs bs + (∇s · Ts + 2Hnβγ · Ts ) − nβγ · (Tβ − Tγ )] dA (2.62)
Aβγ
54
Notas de clase: Francisco J. Valdés Parada
De donde extraemos la condición de salto de cantidad de movimiento entre dos fases:
∂ρs vs
+ ∇s · (ρs vs vs ) + 2Hρs vs (w · nβγ ) = nβγ · [ρβ vβ (vβ − w) − ργ vγ (vγ − w)]
| ∂t
{z } | {z } |
transporte
{z } |
efecto de cambios de la
{z }
flujo de momento desde las fases
acumulación superficial superficie divisoria hacia la superficie divisoria
superficial
Para casos en los que la acumulación, transporte y esfuerzos volumétricos superficiales sean
despreciables respecto a los otros términos en la condición de salto, la ecuación anterior se reduce
a:
nβγ · [ρβ vβ (vβ − w) − ργ vγ (vγ − w)] = − (∇s · Ts + 2Hnβγ · Ts ) + nβγ · (Tβ − Tγ ) (2.64)
Más aún, el tensor de esfuerzos superficiales Ts suele expresarse en términos de la tensión superficial,
γ, como Ts = γI y entonces:
nβγ · [ρβ vβ (vβ − w) − ργ vγ (vγ − w)] = − (∇s γ + 2Hnβγ γ) + nβγ · (Tβ − Tγ ) (2.65)
55
Notas de clase: Francisco J. Valdés Parada
Ejercicio 2.4.
Demuestre que, si no hay transferencia de masa a través de la superficie divisoria, la condición de
salto de cantidad de movimiento se reduce a:
Ejercicio 2.5.
Explique bajo qué condiciones el término ∇s γ = 0 en la condición de salto de cantidad de
movimiento.
Ejercicio 2.6.
Explique bajo qué condiciones la condición de salto de cantidad de movimiento se reduce a la
continuidad de esfuerzos:
nβγ · Tβ = nβγ · Tγ (2.67)
Ejercicio 2.7.
La presión capilar suele definirse como pc = 2Hγ. Explique bajo qué condiciones la condición de
salto de cantidad de movimiento se reduce a
pc = pγ − pβ (2.68)
56
Notas de clase: Francisco J. Valdés Parada
2.7. Esfuerzos viscosos
Hasta este punto hemos llevado nuestros desarrollos sin especificar una forma particular del
tensor de esfuerzos T. Para un fluido en movimiento, dicho tensor se descompone en su contribución
isotropa (−pI) y su contribución viscosa (τ ):
T = −pI + τ (2.69)
Note que una consecuencia directa de la simetrı́a del tensor de esfuerzos es que el tensor de esfuerzos
viscosos sea también simétrico, esto es:
τ = ττ (2.70)
Sustituyendo la Ec. (2.69) en la primera ecuación de Cauchy da lugar a la siguiente expresión:
Dv
ρ = −∇p + ρb + ∇ · τ (2.71)
Dt
que se conoce como la Ecuación de movimiento de Cauchy. A menudo es conveniente escribir esta
ecuación en términos de la presión dinámica, P, la cual se define como:
∇P = ∇p − ρb (2.72)
57
Notas de clase: Francisco J. Valdés Parada
Por último, para el caso en que el flujo sea incompresible, la ecuación de continuidad indica que
∇·v =0 (2.77)
y por tanto la Ec. (2.76) se reduce a:
h i
τ = 2µD = µ ∇v + (∇v)T (2.78)
Ejercicio 2.9.
Para obtener la ecuación de Navier Stokes en la forma dada en la Ec. (2.80) se supuso que el
coeficiente de viscosidad es constante. Sin involucrar esta suposición, explique bajo qué condiciones
se puede aún obtener esta ecuación para casos en los que la viscosidad dependa de la posición.
Ejercicio 2.10.
De acuerdo a la ecuación de continuidad, el orden de magnitud del término ∇ · v es:
1 ∂ρ
∇·v =O + ∇ ln ρ · v (2.81)
ρ ∂t
Tomando en cuenta este estimado, deduzca bajo qué suposiciones es posible recuperar la Ec. (2.80)
sin suponer que el flujo es incompresible.
Ejercicio 2.11.
El vector de vorticidad se define como:
ω =∇×v (2.82)
En este ejercicio debes demostrar que tanto el tensor de vorticidad, Ω, como el vector de vorticidad,
ω, contienen la misma cantidad de información. Para lograrlo, puedes demostrar que en coordenadas
cartesianas el tensor de vorticidad puede expresarse como:
0 ω3 −ω2
1
Ω= −ω3 0 ω1 (2.83)
2
ω2 −ω1 0
58
Notas de clase: Francisco J. Valdés Parada
2.7.2. Fluidos Newtonianos generalizados
De acuerdo a la ley de Newton de la viscosidad, la deformación de un fluido es linealmente
proporcional al esfuerzo aplicado. En la práctica existen varios fluidos que no obedecen esta relación
lineal. Sin embargo, la relación de esfuerzo y deformación aún puede expresarse de acuerdo a la
primera igualdad de la Ec. (2.78) con la diferencia de que ahora el coeficiente de viscosidad será
dependiente del esfuerzo. Para tener más claro lo anterior, es necesario contar con una medida
escalar del tensor de deformación, para ello definimos:
p
γ= 2 (D : D) (2.84)
τ = 2µ(γ)D (2.85)
Por lo que µ ya no puede suponerse como una constante y recibe el nombre de viscosidad aparente.
Bajo estas condiciones, la ecuación de Navier-Stokes es ahora:
Dv
= −∇p + ρg + ∇ · µ(γ) ∇v + (∇v)T
ρ (2.86)
Dt
Existen diversas expresiones de la viscosidad aparente en la literatura, una de las más populares es
el modelo de ley de potencia también conocido como modelo de Ostwald-de Waele (Ostwald, 1929):
Por último, un modelo que engloba a los plásticos de Bingham y a la ley de potencias es el modelo
de Herschel y Bulkley (1926):
59
Notas de clase: Francisco J. Valdés Parada
Introduction 23
Shear thinning fluids have the interesting property that the apparent viscosity decreases with increasing
rate of deformation, i.e., as the rate of deformation is increased the fluid becomes thinner. The apparent
viscosity, µ app , is determined by the relation
Figura 2.4: Bocina con un fluido no newtonianos en cámara lenta. Ver video en:
Tyx = µ app d v x dy
https://www.youtube.com/watch?v=RkLn2gR7SyE ( ) (1-51)
which can be determined by curve shown in Figure 1-10. One can always find a region of linear behavior
at low rates of deformation (think about a Taylor series14 expansion in d v x / dy ), thus many shear thinning
plastic
Bingham
g
i nnin
h
ar t
she
Tyx
ng
ian en
i
w to n ick
Ne r th
hea
s
dvx
dy
Figura 2.5:Figure
Relación1-10. Behavior
de esfuerzo of non-Newtonian
y deformación fluids
en los diversos tipos de fluidos
fluids exhibit what are referred to as lower and upper Newtonian regions. High molecular weight liquids,
such as polymers, often exhibit shear thinning behavior and this represents an important characteristic in
60
the design of coating devices where shear thinning behavior exacerbates coating film non-uniformities.
In addition to shear thinning fluids, there are fluids which behave in precisely the opposite manner, i.e.,
the apparent viscosity increases with increasing rate of deformation. One convenient example is a
Notas de clase: Francisco J. Valdés Parada
Ejercicio 2.12.
Demuestre, usando notación indicial, que el escalar D : D es igual a la traza del tensor D · D, esto
es, D : D = tr(D · D).
Ejercicio 2.13.
Definiendo r
1
τ= tr(τ · τ )
2
Demuestre que
τ = µ(γ)γ
Ejercicio 2.14.
El inverso de la viscosidad aparente se conoce como fluidez, ϕ. Demuestre que la ley de potencias
puede expresarse como:
ϕ = m−1/n τ (1−n)/n
Ejercicio 2.15.
El modelo de Ellis está dado por la siguiente expresión:
" α−1 #
1 τ
ϕ= 1+
µ0 τ1/2
donde los parámetros µ0 , τ1/2 y α se fijan empı́ricamente. Explique bajo qué condiciones este modelo
se reduce al modelo de Ostwald de Waele.
61
Notas de clase: Francisco J. Valdés Parada
Para entender el significado fı́sico de estos números, es conveniente reconocer que la razón
L∗ /v ∗ puede entenderse como un tiempo caracterı́stico asociado a la convección de cantidad de
movimiento, tc ; mientras que la razón L∗ 2 /ν puede asociarse a un tiempo caracterı́stico de difusión
de momento, tν . Es decir, tν indica el tiempo que tardan los esfuerzos viscosos en transmitirse a lo
largo del sistema. De esta forma, el número de Reynolds es la razón entre el tiempo caracterı́stico
de difusión de momento respecto al tiempo caracterı́stico de convección de momento:
L∗ 2 /ν tν
Re = ∗ ∗
= (2.92)
L /v tc
Claramente, si Re > 1, toma más tiempo a la difusión de momento transmitirse que a la convección
y el proceso estará gobernado por esfuerzos inerciales, mientras que en el caso contrario prevalecerán
los esfuerzos viscosos.
De acuerdo a lo anterior, el número de Strouhal puede expresarse como:
t∗
St = (2.93)
tc
Es decir, es la razón entre el tiempo caracterı́stico de transporte en el sistema entre el tiempo
caracterı́stico del transporte de momento por convección. Entonces, un proceso será transitorio si
St 1 y estará en estado estacionario (o en ocasiones cuasi-estacionario) cuando St 1.
El número de Ruark (cuyo inverso se conoce como número de Euler, Eu = 1/Ru) es una relación
de la energı́a cinética por unidad de volumen respecto a la presión aplicada sobre el fluido. En un
proceso sin fricción, Ru = 1.
Por último, el número de Froude es una relación entre la velocidad debida a la inercia y la
magnitud de la fuerza volumétrica debida a la gravedad. Si dicha fuerza es despreciable, se tiene
F r 1.
Consideremos ahora las diferentes versiones de la ecuación de Navier Stokes para dos casos
lı́mite:
62
Notas de clase: Francisco J. Valdés Parada
flujo de fluidos no viscosos (µ = 0). Bajo esta condición, la ecuación de Navier Stokes se reduce a
la ecuación de cantidad de movimiento de Euler:
∂v
ρ + ρv · ∇v = −∇p + ρg (2.96)
∂t
Adicionalmente, para casos en los que la energı́a cinética por unidad de volumen del fluido supere
a la presión, esto es, para Ru 1, la ecuación anterior se reduce a:
∂v
ρ + ρv · ∇v = ρg (2.97)
∂t
Otro caso lı́mite de la ecuación de Euler (Ec. 2.96) corresponde a situaciones en las que St 1 y
corresponde a transporte en estado estacionario (o incluso en estado cuasi-estacionario):
ρv · ∇v = −∇p + ρg (2.98)
La cual es una ecuación diferencial no lineal de primer orden para la velocidad y que dificulta su
solución analı́tica.
Para situaciones de flujo que no se encuentren entre los extremos del flujo reptante y el
potencial se debe utilizar la forma completa de la ecuación de Navier-Stokes. Esta ecuación se
puede simplificar en base a los valores de los números de Strouhal, Ruar o Froude.
Ejercicio 2.16.
Explique bajo qué condiciones la ecuación de Stokes conduce a la ecuación de Laplace para la
presión:
∇2 p = 0
Ejercicio 2.17.
Explique bajo qué condiciones la ecuación de Stokes conduce a la ecuación de Laplace para el vector
de vorticidad (ω = ∇ × v)
∇2 ω = 0
∂(ρv)
+ ∇ · (ρvv) = −∇p + ρg + ∇ · τ (2.100)
∂t
Haciendo el producto punto con el vector de velocidad v,
∂(ρv)
· v + ∇ · (ρvv) · v = −(∇p) · v + ρg · v + (∇ · τ ) · v (2.101)
∂t
63
Notas de clase: Francisco J. Valdés Parada
donde:
Además,
(∇ · τ ) · v = ∇ · (τ · v) − τ T : ∇v (2.106)
Por lo que
∂ 12 ρv 2
1 2
+∇· ρv v = −∇ · (pv) + p(∇ · v) + ρg · v + ∇ · (τ · v) − τ T : ∇v (2.107)
∂t 2
3. −∇ · (pv): Tasa de trabajo de presión hecho por los alrededores sobre el fluido.
Ejercicio 2.18.
ρ ∂v 2
µ
+ ∇ · (vv ) = −∇ · (pv) + ρg · v + ∇2 v 2 − µ∇vT : ∇v
2
(2.108)
2 ∂t 2
64
Notas de clase: Francisco J. Valdés Parada
Una forma alternativa de la ecuación de energı́a mecánica resulta de sustituir al vector de
gravedad por el gradiente de la energı́a potencial por unidad de masa (φ̂ = gh) como:
g = −∇φ̂ (2.109)
Entonces:
ρg · v = −∇φ̂ · ρv = −∇ · (ρvφ̂) + φ̂∇ · (ρv) (2.110)
Con base en la ecuación de continuidad, el resultado anterior se puede expresar como:
∂ φ̂ρ
ρg · v = −∇ · (ρvφ̂) − (2.111)
∂t
donde se supuso que la energı́a potencial no cambia con el tiempo. Bajo estas condiciones, la
ecuación de energı́a mecánica toma la siguiente forma:
∂ 1 2 1 2
ρv + ρφ̂ + ∇ · ρv + ρφ̂ v = −∇ · (pv) + p(∇ · v) + ∇ · (τ · v) − τ T : ∇v (2.112)
∂t 2 2
65
Notas de clase: Francisco J. Valdés Parada
66
Notas de clase: Francisco J. Valdés Parada
3 | Hidrostática
Contenido
3.1. Ecuación fundamental de la hidrostática . . . . . . . . . . . . . . . 67
3.2. Presión hidrostática en lı́quidos y gases . . . . . . . . . . . . . . . . 68
3.3. Manometrı́a . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71
3.4. Fuerzas sobre superficies . . . . . . . . . . . . . . . . . . . . . . . . . 75
Tβ = −Ipβ + τβ (3.3)
| {z } |{z}
esfuerzos normales esfuerzos tangenciales
67
Notas de clase: Francisco J. Valdés Parada
donde se consideró a la gravedad como la única fuerza volumétrica.
Antes de continuar es interesante notar que, para un sólido (fase σ), las ecuaciones que rigen la
transferencia de masa total y de cantidad de movimiento son:
∂ρσ
+ ∇ · (ρσ vσ ) = 0 (3.6a)
∂t
∂vσ
ρσ + ρσ vσ · ∇vσ = −∇pσ + ρσ g + ∇ · τσ (3.6b)
∂t
donde τσ es el tensor de deformación, el cual se expresa, mediante la ley de Hooke, como
∇pβ = ρβ g (3.9)
Note que resulta la misma ecuación para casos en los que la fase sólida sea inmóvil e indeformable.
g = −ge3 (3.10)
68
Notas de clase: Francisco J. Valdés Parada
pβ = −ρβ gz + c (3.13)
Note que para llegar a este resultado se supuso que la densidad del fluido es constante. Más adelante
presentaremos un ejercicio donde este no es el caso. El valor de la constante c depende de la situación
particular a estudiar. Por ejemplo, considere el esquema mostrado en la Figura 3.1, si se toma como
condición de frontera:
en z = 0, pβ = pa (3.14)
la ecuación anterior se expresa como:
p β = p a − γβ z (3.15)
Entonces, a una distancia z = b, dado que el fluido es aire, se tiene
pβ = pa − γaire b (3.16)
pβ = pa + γagua h (3.17)
La gravedad especı́fica del aire es 11.8 N/m3 , mientras que la del agua es de 9790 N/m3 , por lo
que un cambio de altura en una columna de agua se traduce en un mayor cambio de presión que
en una columna de aire.
Considere ahora el sistema mostrado en la Figura 3.2, dado que el sistema coordenado se
encuentra en el fondo del lago, la condición de frontera es:
en z = H, p = pa (3.18)
69
Notas de clase: Francisco J. Valdés Parada
z=H
z=h
z=0
Por lo que:
c = pa + ρβ gH (3.19)
y la solución particular de la forma integrada de la ecuación fundamental de la hidrostática es
Ejercicio 3.1.
El lago de Chapala tiene una profundidad de 6m, calcule la presión manométrica y la presión
absoluta en el fondo del lago. Considere la presión atmosférica como 101325 Pa y la densidad del
agua como 1000 kg/m3 :
Ejercicio 3.2.
Deduzca una expresión que prediga el cambio de presión de un gas ideal con la altura bajo
condiciones estáticas, en el cual la temperatura cambia con la posición de acuerdo a la siguiente
relación lineal:
Tβ = T0 − Bz (3.22)
donde T0 es la temperatura en la superficie (z = 0) y B es una constante.
70
Notas de clase: Francisco J. Valdés Parada
Solución: De acuerdo a la ecuación de los gases ideales, la densidad está relacionada con la
presión como sigue:
M pβ
ρβ = (3.23)
R(T0 − Bz)
Sustituyendo este resultado en la ecuación fundamental de la hidrostática se obtiene
dpβ M gdz
=− (3.24)
pβ R(T0 − Bz)
que al integrarla da lugar a:
Mg
Bz RB
pβ = pa 1 − (3.25)
T0
Mg
Para el aire M = 0.028 kg/mol y B = 0.0065 K/m, por lo que RB = 5.1.
Ejercicio 3.3.
Utilice el resultado del problema anterior para predecir la presión del aire a una altura de 10 km,
para un lugar donde la temperatura en la superficie sea de 30 o C y compara su resultado con el
correspondiente a condiciones isotérmicas.
Solución. Utilizando el resultado del problema anterior tenemos:
Este resultado tiene un error de 14.7 %. Es de esperarse que este error aumente conforme se
incrementa la altura desde la superficie a la cual se desea calcular la presión.
3.3. Manometrı́a
Considere el sistema formado por aceite, agua, glicerina y mercurio como el mostrado en la
Figura 3.3. Para el aceite, la ecuación fundamental de la hidrostática da como resultado:
O bien, evaluado en z = z2 :
p2 − p1 = −ρo g(z2 − z1 ) (3.27)
Para el agua, el resultado es:
p3 − p2 = −ρw g(z3 − z2 ) (3.28)
71
simply having the pressure increase downward and decrease upward. Thus they
following mnemonic, or memory, device, first suggested to the writer by Professo
Oil, ρo
z2 p2 – p1 = – ρog(z 2 – z1)
Water, ρw
z z3 p3 – p2 = – ρw g(z 3 – z 2)
Glycerin, ρG
z4 p4 – p3 = – ρG g(z 4 – z 3)
Mercury, ρM
g. 2.8 Evaluating pressure z5 p5 – p4 = – ρM g(z 5 – z 4)
anges through a column of multi- Sum = p5 – p1
e fluids.
O bien,
p3 − p1 = −ρo g(z2 − z1 ) − ρw g(z3 − z2 ) (3.29)
Siguiendo este procedimiento, se obtiene:
O bien
p1 + ρo g(z1 − z2 ) + ρw g(z2 − z3 ) + ρG g(z3 − z4 ) + ρw g(z4 − z5 ) = p5 (3.31)
Note que una regla práctica para recordar este resultado es que se empieza en un punto donde se
conoce la presión y se agrega el valor de ρg y la altura correspondiente a cada sección de fluido
hasta llegar al punto deseado.
72 ChapterEjercicio
2 Pressure3.4.
Distribution in a Fluid
Flow device
(a) (b)
1 h
2
E2.3
Solution
Solution
Using our “up-down” concept as in Eq. (2.32), start at (a), evaluate pressure changes around the
U-tube, and end up at (b):
Notas
pa de
1gLclase:
Francisco J. Valdés Parada
1gh 2gh 1gL pb
Considere el sistema mostrado en la figura, calcule la diferencia de presión entre los puntos a y b.
or Solución: pa pb (2 1)gh Ans.
ρ3
Jump across
z 2, p2 z 2, p2
ρ1
zA, pA A
B zB, pB
Jump across
z1, p1 z1, p1
ed multiple- Jump across
ate pA to pB. z 3, p3 z 3, p3
ecially prac- ρ2
homework
ρ4
m.
Considere el sistema mostrado en la figura, calcule la diferencia de presión entre los puntos a y b.
Solución:
Ejercicio 3.6.
73
EXAMPLE 2.4
Pressure gage B is to measure the pressure at point A in a water flow. If the pres
kPa, estimate the pressure at A, in kPa. Assume all fluids are at 20°C. See Fig. E
Notas de clase: Francisco J. Valdés Parada
Mercury 6 cm
A
5 cm
Water
flow 11 cm
4 cm
Fluid Statics 43
E2.4
dA
z w dx (2-81)
Solution
Considere el sistema mostrado en la figura, calcule la presión en el punto A sabiendo que p =
B
3 = 133, 100 N/m3 , γaceite = 8720 N/m3
and this allows87kPa,
us toγexpress Eq.N/m
agua = 9790 2-80, γinmercurio
the form
First list the specific weights from Table 2.1 or Table A.3:
Solución: La fórmula del cambio de presión entre los puntos A y B es:
xL
water 9790 N/m3 mercury 133,100 N/m3 oil 8720 N/
k fpfluid w
g h(0.07)
a + γagua (0.05) − γmercurio 1 h−
( xγ)aceite
dx(0.06) = pB (2-82)
Now proceed fromx A0 to B, calculating the pressure change in each fluid and add
de donde se obtiene que ppA W(z)W M(z)M O(z)O pB
= 96.35kPa.
A
Here we have made use of the representation, z h( x) , to locate the fluid-solid interface as indicated in Figure
A x(9790 (133,100 that (8720 N/m3)(
3 3
2-10 where we have also illustrated
or pthat rangesN/m
from 0 to
)( 0.05 A little
L. m) thought will
N/mindicate
)(0.07 m)
where we replace N/m2 by its short name, Pa. The value zM 0.07 m is the
change in the mercury (11 cm 4 cm). Solving for the pressure at point A, we
pA 96,351 Pa 96.4 kPa
The intermediate six-figure result of 96,351 Pa is utterly fatuous, since the meas
cannot be made that accurately.
2-10. Volume
Figure Figura ofsobre
3.4: Fuerza wateruna
located above
presa de alturathe
h1 barrier
k f fluid
g V (2-83)
Notas de clase: Francisco J. Valdés Parada
3.4. Fuerzas sobre superficies
Considere el sistema mostrado en la Figura 3.4; de acuerdo a la ecuación fundamental de la
hidrostática, la presión dentro del agua en la presa es:
pβ = pa − ρβ g(z − h1 ) (3.32)
Se desean conocer las componentes del vector de fuerza que ejerce el fluido sobre el sólido. Para
ello se recurre a la siguiente definición:
Z
f = n · TdA (3.33)
A
T = −pI (3.34)
por lo tanto, Z Z Z
f=− npdA = − npa dA + nρβ g(z − h1 )dA (3.35)
A A A
El primer término de la ecuación anterior es cero. Este resultado puede verificarse con ayuda del
teorema de la divergencia: Z Z
npa dA = ∇pa dV = 0 (3.36)
A V
Para calcular las componentes del vector de fuerza, se hace el producto punto con las componentes
x y z de la Ec. (3.35), el resultado para la componente x es:
Z
fx = ρβ g(z − h1 )i · ndA (3.37)
A
Aquı́ de tomó en cuenta el hecho de que la normal va dirigida de la superficie de la presa hacia
el agua, por lo que su componente x es negativa. Pensando que la presa tiene un espesor w, el
resultado es ahora
z=h
Z 1
h1
fx = −w ρβ g(z − h1 )dz = ρβ g wh1 (3.39)
2
z=0
Es decir, es el esfuerzo evaluado en el centroide de la presa multiplicada por el área perpendicular
a la dirección x.
Dirigiendo ahora la atención a la componente z del vector de fuerza, tenemos
Z
fz = ρβ g(z − h1 )k · ndA (3.40)
A
75
Notas de clase: Francisco J. Valdés Parada
dada la geometrı́a de la presa, tenemos que
Entonces:
x=L
Z
fz = w ρβ g(z − h1 )dx (3.42)
Chapte
x=0
Para evaluar esta ecuación se necesita una forma explı́cita de cómo cambia la coordenada z con x;
por ejemplo si z = h1 x/L, el resultado es
x=L
Z x
fz = w ρβ gh1 − 1 dx = −(wh1 L/2)ρβ g (3.43)
L
x=0
Ejercicio 3.7.
Repita el cálculo de fz para el caso el que se use una aproximación polinomial de segundo orden
para describir la forma de la presa. Para ello, proponga un polinonio del tipo:
z = a + bx + cx2
x = 0, z=0
x = L, z = h1
Figure 2.4a.dzLog dam
x = L, =0
dx
The total force exerted
Muestre que laat the surface
expresión illustrated
resultante de z es in Figure 2.4b can be expressed as
h1 x
f surface z = 2 −
L t (n ) dA
L
x (1)
Recalcule fz y compare el resultado con el que
A se obtuvo en la ecuación (3.43).
and if we take n to be
Ejercicio the outwardly directed normal vector, f surface represents the force exerted by
3.8.
76
Figure 2.4b. Eulerian cut
the surroundings on the log. Since we are interested only in forces over and above those created by the
constant ambient stress, we make use of Eq. 2-68 to express Eq. 1 as
Notas de clase: Francisco J. Valdés Parada
Calcule la fuerza horizontal por unidad de la longitud del tronco mostrado en la figura tomando
H = 1.6m, α = 0.375, ρβ = 997 kg/m3 .
Solución.
z=H z=αH
fx 1
Z Z
= −γ (H − z)dz − (αH − z)dz = − γH 2 (1 − α2 ) = −10, 758N/m
w 2
z=0 z=0
77
Notas de clase: Francisco J. Valdés Parada
78
Notas de clase: Francisco J. Valdés Parada
Contenido
4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79
4.2. Flujo entre dos placas paralelas (Flujo de Poiseuille) . . . . . . . . 80
4.3. El problema de Beavers y Joseph . . . . . . . . . . . . . . . . . . . . 84
4.4. Flujo de dos fluidos inmiscibles entre dos placas paralelas . . . . . 86
4.5. Flujo entre dos placas que liberan gas . . . . . . . . . . . . . . . . . 90
4.6. Flujo en un tubo cilı́ndrico (Flujo de Hagen-Poiseuille) . . . . . . 93
4.7. Flujo en la sección anular de un tubo . . . . . . . . . . . . . . . . . 95
4.8. Viscosı́metro de Couette . . . . . . . . . . . . . . . . . . . . . . . . . 101
4.1. Introducción
Es conveniente comenzar este capı́tulo explicando cada elemento de su tı́tulo. Por flujo laminar
entendemos aquel flujo donde los esfuerzos inerciales y los viscosos son importantes pero la tasa
de disipación de energı́a compensa a la tasa de energı́a cinética del fluido. Por flujo uniforme nos
referimos a que el vector de velocidad sólo tiene una componente que no es cero y es aquella que va
en la misma dirección del flujo. Por flujo unidimensional nos referimos a que dicha componente del
vector de velocidad que no es nula depende de una sola variable independiente de la posición. Por
último, por flujo en estado estacionario, entendemos que la velocidad es independiente del tiempo.
Además de las suposiciones anteriores, en la mayor parte de este capı́tulo abordaremos problemas
que involucran fluidos Newtonianos.
Las ecuaciones gobernantes del transporte de masa y cantidad de movimiento bajo estas
condiciones son:
∇·v= 0 (4.1a)
2
ρv · ∇v = −∇p+ ρg + µ∇ v (4.1b)
79
d These three axioms, along with the laws of thermodynamics (see Chapter 11), provide the framework for
our analysis of fluid flow phenomena.
i t (k ) 0 , y0 (3-2
3.1 Uniform laminar flow in a channel.
In this section we consider laminar flow in an infinitely wide channel caused by the application of a
om this we conclude that the shear stress at the centerline must be zero,
pressure gradient. This flow is illustrated in Figure 3-1 where the length of the entrance region is
Notas de clase: Francisco J. Valdés Parada
i t (k ) 0 , yLe 0 , or zx 0 , no-slip
y0 (3-
z
d Eq. 3-20 simplifies to
p x
zx gx z (3-
x
r negative values of z this result indicates that zx is positive, while for positive values of z we see tha
no-slip
91
∂ 2 vx v
x ∂ 2 vx v
x
= O ; = O (4.6)
∂y 2 W2 ∂z 2 h2
Por lo que, en base a la desigualdad de longitudes:
h2 W 2 (4.7)
∂p d2 vx
0=− +µ 2 (4.9)
∂x dz
Aplicando la condición de no deslizamiento en las superficies de las placas, tenemos las siguientes
dos condiciones de frontera:
En z = ±h, vx = 0 (4.10)
Antes de comenzar con la solución del problema, es necesario contar con información acerca del
gradiente de presión. Suponiendo que se conoce la presión a la entrada y salida del sistema, se tiene
entonces que:
∂p pL − p0
= (4.11)
∂x L
donde p0 y pL representan las presiones a la entrada y salida del sistema, respectivamente y L es
la longitud de las placas. Como puede notarse, pL < p0 y el gradiente de presión es negativo. De
esta forma, la ecuación de movimiento es ahora:
d2 vx pL − p0
2
= (4.12)
dz µL
81
Notas de clase: Francisco J. Valdés Parada
Dada la simetrı́a del problema, puede demostrarse que el problema puede resolverse solo en la mitad
del sistema, por lo que las condiciones de frontera son:
dvx
En z = 0, = 0; (4.13)
dz
En z = h/2, vx = 0 (4.14)
La solución de este problema de valor a la frontera es
p0 − pL 2 4z 2
vx = h 1− 2 (4.15)
8µL h
El valor de la velocidad máxima se encuentra en z = 0 y es:
p0 − pL 2
vmáx = vx (0) = h (4.16)
8µL
El valor de la velocidad promedio es:
z=h/2 z=h/2 z=h/2
4z 2
Z
1 2 2vmáx 2vmáx h h 2vmáx
Z Z
hvx iβ = vx dz = vx dz = 1− 2 dz = − =
h h h h h 2 6 3
z=−h/2 z=0 z=0
(4.17)
La definición del flujo volumétrico es:
Q̇ = hvx iβ Ax (4.18)
donde Ax = hW es el área de la sección transversal del flujo. Por lo tanto, el flujo volumétrico es:
2vmáx
Q̇ = Wh (4.19)
3
Este resultado puede usarse para conocer el flujo másico:
2ρvmáx
ṁ = ρQ̇ = Wh (4.20)
3
Por último, la fuerza que ejerce el fluido sobre cualquiera de las placas es:
Z Z
f = −npdA + n · µ(∇v + ∇vT )dA (4.21)
Az Az
82
Notas de clase: Francisco J. Valdés Parada
Ejercicio 4.1.
Demuestre que la fuerza ejercida sobre la superficie en z = −h es la misma que en z = h.
Ejercicio 4.2.
dvx
z = 0, =0
dz
h
z= , vx = 0
2
Dado que el problema es lineal, su solución formal se puede expresar en función de la fuente como:
b dp
vx = −
µ dz
donde b = b(z) es una función cuyo problema de valor a la frontera se obtiene de sustituir la
expresión anterior en la ecuación diferencial y condiciones de frontera para vx , el resultado es
d2 b
1+ =0
dz 2
db
z = 0, =0
dz
h
z= ,b = 0
2
Muestre que la solución de este problema es
h2 z 2
b= −
8 2
Por lo tanto, la solución formal es ahora
h2 z 2
1 dp
vx = − −
µ 8 2 dz
O bien, al aplicarle el operador de promediado espacial
h2 dp
hvx iβ = −
12µ dz
Este resultado representa el promedio intrı́nseco de la velocidad. La ley de Darcy, está expresada
en términos del promedio superficial de la velocidad:
K dp
hvx i = −
µ dz
83
Notas de clase: Francisco J. Valdés Parada
donde los operadores de promedio superficial e intrı́nseco están relacionados por
hvx i = εhvx iβ
donde ε es la fracción volumétrica del fluido contenido en una celda unitaria. En el caso de estudio,
dicha celda puede suponerse que tiene una altura H, por lo que ε = h/H. Entonces, puede concluirse
que, en este caso, la permeabilidad es
K ε3
=
H2 12
Evidentemente, en cuanto más separadas estén las placas mas permeable será el sistema.
dp d2 vx
0=− +µ 2 (4.24)
dx dz
Mientras que en el medio poroso aplica la ley de Darcy:
K dp
vD = − (4.25)
µ dx
donde K representa la permeabilidad del medio poroso. Dado que el gradiente de presión se supone
constante, la velocidad de Darcy es constante. Las condiciones de frontera para este problema son:
En z = h, vx = 0 (4.26)
dvx α α K dp
En z = 0, = √ (vx − vD ) = √ vx + (4.27)
dz K K µ dx
Donde α es una constante conocida. Note que cuando el medio poroso es impermeable (K = 0), la
condición de frontera anterior se reduce a la condición de no deslizamiento.
La solución general de este problema es:
1 dp 2
vx = z + cz + c2 (4.28)
2µ dx
z2
1 dp 2 z
vx = − h 1 − 2 − ch 1 − (4.29)
2µ dx h h
84
etain the parallel flow assumption, which leads, through the continuity equatio
o the independence of u upon x. It follows through (2) that /3 also does n
epend on x,and so depends only on the properties of the fluid and the permeab
material. Notas de clase: Francisco J. Valdés Parada
Permeable
material
IGURE
1. Velocity profile for t h e rectilinear flow in a horizontal channel formed by
Figura 4.2: Sistema de Beavers y Joseph que onsiste en un canal parcialmente lleno con un medio
permeable
poroso. lower wall (y = 0) and an impermeable upper wall (y = h).
ndependent quantity containing the dimensions 85 of mass and time. This sugge
hat p is independent of the viscosity of the fluid, a result which is strong
mplied by our experimental observations. Furthermore, a length scale chara
Notas de clase: Francisco J. Valdés Parada
el resultado " #
dvx 1 dp 2z α 2 − σ2
= h + =0 (4.34)
dz 2µ dx h (1 + ασ) σ
De esta ecuación se obtiene el valor crı́tico como
α σ2 − 2
z
=λ= (4.35)
h 2 (1 + ασ) σ
Por lo tanto, la velocidad máxima es
1 dp 2
vmáx = − h (1 − λ)2 (4.36)
2µ dx
Ejercicio 4.3.
Demuestre que, para el caso en el que el medio poroso sea impermeable, es decir, cuando K →
0(σ → ∞), la velocidad máxima se localiza a la mitad del canal.
Ejercicio 4.4.
Muestre que la velocidad promedio en el canal es
" #
1 h2 4 α 2 − σ 2 dp
hvx i = − +
µ 4 3 (ασ + 1) σ dx
Ejercicio 4.5.
Muestre que las fuerzas en las superficies del canal son:
" #
α 2 − σ2
dvx hW L dp
fx |z=h = −µW L =− 2+
dz z=h 2 (ασ + 1) σ dx
" #
hW L α 2 − σ 2
dvx dp
fx |z=0 = µW L =
dz z=0 2 (ασ + 1) σ dx
Como puede notarse, la superficie en contacto con el medio poroso tiene una menor fuerza de
fricción. Por ello, es que en algunos diseños de aviones, álabes y autos se considera recubrir sus
superficies de medio poroso para reducir la fricción. Note que en el caso en que σ → ∞ ambas
fuerzas se igualan.
86
This equation is obtained for both phase I and phase 11. Integration of Eq. 2.5-1 for the
two regions gives
Velocity
distribution,
Shear stress
or momentum-
flux distribution
Fig. 2.5-1 Flow of two immiscible fluids between a pair of horizontal plates under
Figura 4.3:ofFlujo
the influence de dos
a pressure fluidos inmiscibles entre dos placas paralelas.
gradient.
The adjacent
Las cuales están sujetas a lasflow of gases and condiciones
siguientes liquids in conduits been reviewed by A. E. Dukler and M.
dehasfrontera:
Wicks, 111, in Chapter 8 of Modern Chemical Engineering, Vol. 1, "Physical Operations," A. Acrivos (ed.),
Reinhold, New York (1963).
En z = −h/2; vI = 0 (4.37c)
En z = h/2; vII = 0 (4.37d)
En z = 0;
vI = vII (4.37e)
dvI dvII
En z = 0; µI = µII (4.37f)
dz dz
Las tres primeras condiciones de frontera se deducen de la condición de no deslizamiento, la cual
implica que la velocidad del fluido es la velocidad de la cpa adyacente. La última condición de
frontera se deduce de la condición de salto de transporte de cantidad de movimiento, la cual bajo
condiciones de flujo en estado estacionario puede escribirse como sigue:
∇ · (ρ v v ) + 2Hρ v (w · nI,II ) = nI,II · [ρI vI (vI − wI,II ) − ρII vII (vII − wI,II )]
| s {zs s s} | s s {z } | {z }
transporte efecto de cambios de la flujo de momento desde las fases
superficial superficie divisoria hacia la superficie divisoria
∇s · (ρs vs vs ) · ex = nI,II · (ρI vI vI − ρII vII vII ) · ex + ∇s · Ts · ex − nI,II · (TI − TII ) · ex (4.39)
Note que la contribución volumétrica se eliminó porque la gravedad actúa en una dirección
perpendicular al eje x. Ahora bien, dado que nI,II = ez y se ha hecho la suposición de flujo
desarrollado, la ecuación anterior se reduce a
∇s · (ρs vs vs ) · ex = ∇s · Ts · ex − nI,II · (TI − TII ) · ex (4.40)
Dirigiendo la atención al término del lado izquierdo, recordando la condición de incompresibilidad,
∂ ∂
que ∇s = ∇ · (I − nI,II nI,II ) = ∂x ex + ∂y ey y que vs = vx ex + vy ey , tenemos entonces lo siguiente
∂ ∂
∇s · (ρs vs vs ) · ex = ∇s · (ρs vs vx ) = ρs (vx vx ) + (vy vx ) = 0 (4.41)
∂x ∂y
87
Notas de clase: Francisco J. Valdés Parada
Si además se desprecia la contribución de los esfueroz superficiales respecto a los esfuerzos normales,
la condición de salto se reduce a:
0 = ez · (TI − TII ) · ex (4.42)
Por último, al sustituir la definición del tensor total de esfuerzos, resulta la condición de frontera
dada en la Ec. (4.37f).
Al integrar una vez las ecuaciones diferenciales obtenemos
dvI (p0 − pL )z
=− + c1 (4.43a)
dz µI L
dvII (p0 − pL )z
=− + c3 (4.43b)
dz µII L
Por lo que las soluciones generales de este problema son:
(p0 − pL )z 2
vI = − + c1 z + c2 (4.44a)
2µI L
(p0 − pL )z 2
vII =− + c3 z + c4 (4.44b)
2µII L
Usando las condiciones de frontera en z = ±h/2 obtenemos:
(p0 − pL )h2 h
c2 = + c1 (4.45a)
8µI L 2
(p0 − pL )h 2 h
c4 = − c3 (4.45b)
8µII L 2
Por lo que las soluciones toman la forma siguiente:
(p0 − pL )h2 4z 2
h 2z
vI = 1 − 2 + c1 +1 (4.46a)
8µI L h 2 h
(p0 − pL )h2 4z 2
h 2z
vII = 1 − 2 + c3 −1 (4.46b)
8µII L h 2 h
Más aún, en base a la condición de frontera de la continuidad de los esfuerzos, tenemos que:
µI c1 = µII c3 (4.47)
(p0 − pL )h2 4z 2
h 2z
vI = 1 − 2 + c1 +1 (4.48a)
8µI L h 2 h
(p0 − pL )h2 4z 2
µI h 2z
vII = 1 − 2 + c1 −1 (4.48b)
8µII L h µII 2 h
Para calcular el valor de la última constante a determinar, se utiliza la condición de frontera de
continuidad de la velocidad:
(p0 − pL )h2 h (p0 − pL )h2 µI h
+ c1 = − c1 (4.49)
8µI L 2 8µII L µII 2
88
Notas de clase: Francisco J. Valdés Parada
O bien
(p0 − pL )h2 µI − µII
h
c1 = (4.50)
2 8µI L µI + µII
Por lo tanto,
En este caso, no hay simetrı́a en el sistema, por lo que es necesario calcular las velocidades promedio
en cada sección del sistema y los valores de la fuerza sobre cada superficie. Las velocidades promedio
en cada fase son:
z=0
2
Z
I
hvI i = vI dz (4.52a)
h
z=−h/2
z=h/2
2
Z
hvII iII = vII dz (4.52b)
h
z=0
Ejercicio 4.6.
Demuestre que la posición z = δ donde se ubica la velocidad máxima es:
h µI − µII
δ=
4 µI + µII
Ejercicio 4.7.
Calcule los valores de la componente x del vector de fuerza en las superficies z = h/2 y z = −h/2.
Respuesta:
(p0 − pL )W h µI + 3µII
∂vII
fx |z=h/2 = −µII W L =
∂z z=h/2 4 µI + µII
(p0 − pL )W h 3µI + µII
∂vI
fx |z=−h/2 = µI W L =−
∂z z=−h/2 4 µI + µII
89
Notas de clase: Francisco J. Valdés Parada
Figura 4.4: Flujo entre dos placas paralelas que liberan un gas de manera asimétrica.
y = h, v = −v0+ ey (4.54)
Se desean obtener las expresiones para las componentes x y y del vector de velocidad ası́ como de la
presión. Puede suponer que en x = 0 la velocidad es cero. Es decir que hay una pared impermeable.
En este caso las ecuaciones de continuidad y de cantidad de movimiento son
∂vx ∂vy
+ =0 (4.56a)
∂x ∂y
∂p
0=− + µ∇2 vx (4.56b)
∂x
∂p
0=− + µ∇2 vy (4.56c)
∂y
donde el Laplaciano aplica en 2D. Sin embargo, suponiendo que el largo, L, del canal es mucho
mayor que su ancho, 2h, se pueden despreciar los cambios de la velocidad en la dirección x respecto
a los cambios en la dirección y y entonces
∂p ∂ 2 vx
0=− +µ 2 (4.57)
∂x ∂y
∂p ∂ 2 vy
0=− +µ 2 (4.58)
∂y ∂y
90
Notas de clase: Francisco J. Valdés Parada
El gradiente de presión en x es independiente de y por lo que
y 2 ∂p
vx = + c1 y + c2 (4.59)
2µ ∂x
Tomando en cuenta las condiciones de frontera sabemos que en y = ±h, vx = 0 y por lo tanto
h2 ∂p y2
vx = − 1− 2 (4.60)
2µ ∂x h
Ahora bien, con la información con la que se cuenta no se sabe cómo cambia la presión en la
dirección y, ni se sabe el valor de la presión, por lo que no se gana mucho al integrar la ecuación
de movimiento en la dirección y. Por ello, es más conveniente dirigir la atención a la ecuación de
continuidad, la cual al integrarla respecto a x da como resultado
∂vy
vx = − x (4.61)
∂y
Aquı́ se tomó en cuenta la suposición que se pueden despreciar los cambios de vy en la dirección x
y la condición de frontera que en x = 0, la velocidad es cero. Bajo estas condiciones, debe ser claro
que vx /x es una función que solo depende de y. Esto es,
h2 ∂p y2
∂vy
= 1− 2 (4.62)
∂y 2µx ∂x h
2
h ∂p +
Note que 2µx ∂x debe ser una constante. Integrando esta ecuación, sujeta a que en y = h, vy = −v0
se tiene que
h3 ∂p h2 ∂p y2
+
vy = −v0 − + 1− 2 y (4.63)
3µx ∂x 2µx ∂x 3h
Usando ahora la condición de frontera que en y = −h, vy = v0− se obtiene que
∂p 3µ
= − 3 v0+ + v0− x
(4.64)
∂x 2h
Con este resultado se pueden expresar ahora los perfiles de velocidad en ambas direcciones como
3 v0+ + v0−
y2
vx = 1− 2 x (4.65)
4h h
∂ 2 vy ∂p 3µ
= 3 v0+ + v0− y
µ 2
= (4.67)
∂y ∂y 2h
Por lo tanto, dada la siguiente expresión que es el resultado de usar la regla de la cadena:
∂p ∂p
dp = dx + dy (4.68)
∂x ∂y
91
Notas de clase: Francisco J. Valdés Parada
es posible calcular a la presión al integrar esta ecuación respecto a un valor de referencia:
y2 x2
3µ +
v0 + v0−
p = pref + − (4.69)
4h h2 h2
Ejercicio 4.8.
Muestre que los valores de la fuerza en la dirección x evaluados en y = ±h son iguales y valen
x=L
Z z=W x=L
Z z=W
3 v0+ + v0− L2
∂vx
Z Z
fx |y=±h = n · T · ex dzdx = ∓µ dzdx = µW
∂y y=±h 4h2
x=0 z=0 x=0 z=0
Ejercicio 4.9.
Determine ahora los valores de la fuerza en la dirección y evaluados en y = ±h y muestre que no
son iguales y se rigen por la siguiente expresión
x=L
Z z=W
x=L ! x=L
∂v
Z Z Z
y
fy |y=±h = n · T · ey dzdx = ∓W − p|y=±h + 2µ dx = ±W
p|y=±h dx
∂y y=±h
x=0 z=0 x=0 x=0
2
3µ + L
v + v0− 1 − 2
= ± W L pref +
4h 0 3h
Ejercicio 4.10.
Considere ahora el caso en que hay acumulación en la fase gaseosa, de manera que la ecuación de
continuidad es ahora
∇ · v = q(y)
donde q es una función continua de y. Suponiendo que en este caso si es razonable imponer la
condición de no deslizamiento (v = 0) en las paredes del sistema, demuestre que la solución del
problema es
y2
3
vx = 1 − 2 hqix
2 h
η=y η=h
2
3 y 1
Z Z
vy = − 1 − 2 hqiy + qdy − qdy
2 3h 2
η=−h η=y
y2 x2
3µ
p = pref + hqi − + µq
2 h2 h2
Donde pref es la presión de referencia en el punto x = 0 y y = 0, donde además puede suponerse
que q(0) = 0.
92
Flow Through a Circular Tube
Notas de clase: Francisco J. Valdés Parada
4rzI r + A r = flux
of z-momentum
out at r + Ar
+ Tube wall
of z-momentum
outatz=L
1 ∂rvr Eq.1(2.3-8)
When we divide ∂vθ by ∂vz2.irLAr and take the limit as Ar + 0, we get
+ + =0 (4.70)
r ∂r r ∂θ ∂z
Pero, dado que en condiciones de flujo desarrollado se tiene que vr = vθ = 0, en este caso, la
ecuación anterior se reduce a:
The expression on the∂vzleft side is the definition of the first derivative of r4,, with res
=0 (4.71)
to r. Hence Eq. 2.3-7∂z
may be written as
Por su parte, la componente z de la ecuación de movimiento es:
1 ∂ 2 vz ∂ 2 vz
∂vz vθ ∂vz ∂vz ∂p 1 ∂ ∂vz
ρ vr + + vz =− + ρgz + µ r + 2 + (4.72)
∂r r ∂θ Now∂z we have ∂z
to evaluate ther components
∂r ∂r ∂θ2
4,r and +,,
∂z 2 Eq. 1.7-1and Appendix B
from
93
Notas de clase: Francisco J. Valdés Parada
note el signo positivo en el término gravitatorio debido a la elección del sistema coordenado.
Tomando en cuenta las consideraciones anteriores, este resultado se simplifica a:
1 ∂ 2 vz
∂p 1 ∂ ∂vz
0=− + ρgz + µ r + 2 (4.73)
∂z r ∂r ∂r r ∂θ2
El segundo término de esfuerzos viscosos puede eliminarse de esta última ecuación al suponer flujo
axisimétrico. Esto último puede probarse al integrar la ecuación anterior desde θ = 0 hasta θ = 2π.
Por lo tanto, la ecuación diferencial simplificada es:
∂p µ ∂ ∂vz
0=− + ρg + r (4.74)
∂z r ∂r ∂r
Antes de continuar, es conveniente especificar que el gradiente de presión está dado por
∂p p0 − pL
=− (4.75)
∂z L
Por lo que introduciendo el cambio de variable
P = p − ρgz (4.76)
tenemos que:
∂P ∂p P0 − PL
− =− + ρg = (4.77)
∂z ∂z L
y entonces, dado que vz sólo depende de r, la ecuación de movimiento se reduce a:
P0 − PL
1 d dvz
r =− (4.78)
r dr dr µL
Las condiciones de frontera a las que está sujeta esta ecuación diferencial son:
dvz
en r = 0, =0 (4.79a)
dr
en r = R, vz = 0 (4.79b)
La primera condición de frontera es una condición de simetrı́a que permite resolver el problema sólo
en la mitad del cilindro y la segunda condición de frontera viene de la condición de no deslizamiento.
La solución de este problema de valor a la frontera es:
(P0 − PL )R2 r2
vz = 1− 2 (4.80)
4µL R
De esta expresión podemos deducir que la velocidad máxima es (vz (0)):
(P0 − PL )R2
vmáx = (4.81)
4µL
Además, la velocidad promedio, es
θ=2π
R r=R
R
vz rdrdθ r=R
2 vmáx
Z
β θ=0 r=0
hvz i = = 2 vz rdr = (4.82)
θ=2π
R r=R
R R 2
rdrdθ r=0
θ=0 r=0
94
Notas de clase: Francisco J. Valdés Parada
Con este resultado obtenemos de inmediato el flujo másico como:
cuya componente z es Z Z
fz = − n · ez p + n · τ · ez dA (4.85)
Ar Ar
donde
dvz (P0 − PL )R
=− (4.87)
dr r=R
2µL
Por lo tanto, la fuerza en z es
Como puede notarse, las fuerzas viscosas están balanceadas con las fuerzas superficiales de presión
y las fuerzas volumétricas gravitacionales.
Ejercicio 4.11.
Muestre que en este caso la permeabilidad está dada por
K ε2
=
H2 8π
Compare este resultado con el obtenido para el flujo de Poiseuille y con el resultante de la ecuación
de Kozeny-Carman
d2p ε3
KKC =
180(1 − ε)2
donde dp = 6V /Aβσ
95
54 Chapter 2 Shell Momentum Balances and Velocity Distributions in Laminar Flow
Notas de clase: Francisco J. Valdés Parada
Fig. 2.4-1 The momentum
and velocity distribution
flow in a cylindrical annu
momentum flux changes
value of r for which the v
Velocity
distribution maximum.
Shear stress
or momentum-
flux distribution
the
Y, bajo condiciones de t uuniforme
flujo b e t h aunidireccional,
t is, in the direction opposed
la ecuación de to segravity.
movimiento reduce a:We make the sa
52.3: v, = v,(r),∂pv, = 0,µv,∂ = 0, and
∂vz
p = p(z). Then when we make a m
0=− − ρg + r (4.90)
over a thin cylindrical
∂z shell
r ∂r of∂rliquid, we arrive at the following differe
Note el signo negativo en el término gravitacional debido a la elección del sistema coordenado. El
gradiente de presión es:
∂p p0 − pL
=− (4.91)
∂z L
This
En este caso, definimos differs
el cambio de from
variableEq.
para2.3-10 only
la presión como: +
in that 9 = p pgz here, since the
co
direction opposed to gravity (i.e., z is the same as the h of footnote 1 i
P = p + ρgz (4.92)
of Eq. 2.4-1 gives
tenemos que:
∂P ∂p P0 − PL
= + ρg = − (4.93)
∂z ∂z L
Entonces, la ecuación de movimiento toma la forma:
just as in Eq.
1 d2.3-1
dv1.
z
P0 − PL
r =− (4.94)
The constant
r dr drC, cannot µLbe determined
immediately, since we ha
about the momentum 96
flux at the fixed surfaces r = KR and r = R. A
there will be a maximum in the velocity curve at some (as yet unknow
which the momentum flux will be zero. That is,
Notas de clase: Francisco J. Valdés Parada
La cual está sujeta a las siguientes condiciones de frontera de no deslizamiento:
en r = κR, vz = 0 (4.95a)
en r = R, vz = 0 (4.95b)
(P0 − PL )R2 r2
2 ln(r/R)
vz = 1− 2 − 1−κ (4.97)
4µL R ln(k)
Para encontrar la velocidad máxima, es conveniente tratar el problema como un problema de
optimización. Para ello, se deriva la ecuación anterior con respecto a r y se iguala a cero el resultado:
" #
dvz (P0 − PL )R2 2r 1 − κ2
= − 2− =0 (4.98)
dr 4µL R r ln(k)
Utilizando el criterio de la segunda derivada puede demostrarse que este valor crı́tico corresponde
a un máximo. La velocidad máxima se obtiene al sustituir r/R = λ en la ecuación (4.97):
(P0 − PL )R2
1 − λ2 1 − ln λ2
vmáx = (4.100)
4µL
(P0 − PL )R2 1 − κ2
2
hvz i = 1+κ − (4.102)
8µL ln (1/κ)
ρπ(P0 − PL )R4 (1 − κ2 )2
4
ṁ = ρAz hvz i = 1−κ − (4.103)
8µL ln (1/κ)
97
Notas de clase: Francisco J. Valdés Parada
Por último, la fuerza que se ejerce sobre cada pared del tubo es:
dvz
fz |r=κR = 2πκRLµ (4.104a)
dr r=κR
dvz
fz |r=R = −2πRLµ (4.104b)
dr r=R
Usando la primera igualdad de la Ec. (4.98), las ecuaciones anteriores toman la siguiente forma:
" #
2 1 − κ2
2
fz |r=κR = π(P0 − PL )R −κ − (4.105a)
2 ln(k)
" #
2 1 − κ2
fz |r=R = π(P0 − PL )R 1 + (4.105b)
2 ln(k)
Ejercicio 4.12.
Deduzca la ecuación (4.97) a partir de la ecuación (4.96).
Ejercicio 4.13.
Considere el sistema mostrado en la figura, en el cual un fluido viscoso que fluye de abajo hacia
arriba se derrama por la parte superior de un tubo cilı́ndrico. Suponiendo condiciones de flujo en
estado estacionario, newtoniano e incompresible, determine el perfil de velocidad por fuera del tubo
mientras va cayendo. Puede despreciar los efectos de entrada y salida.
98
2 3 4
and then performing a long division. The first term in the resulting series will be Eq. 2B.
tion: In the derivation it is necessary to use the first four terms of the Taylor series in Eq
2B.7 Flow of a film on the outside of a circular tube. In a gas absorption experiment, a visco
Velocity
distribution
inside tube Velocity
distribution
outside
in film
∆r
z-Momentum z-Momentum
into shell out of shell
L of thickness
∆r ∆r
of thickness
Gravity force
acting on
2πr∆rL
the volume
CL
Fig. 2B.7 Velocity distribution vz (r) and z
R
momentum balance for the flow of a fallin
aR on the outside of a circular tube.
Respuesta:
ρgR2 r2
2 r
vz = 1 − 2 + 2a ln
4µ R R
Ejercicio 4.14.
Considere el sistema mostrado en la figura, en el cual un tubo cilı́ndrico de radio κR se mueve para
arrastrar el fluido contenido en otro compartimento cilı́ndrico de radio R sin que esto genere un
gradiente de presión en la dirección horizontal. Muestre que el perfil de velocidad es
vz ln r ln r
= 1 + κR = R
v0 ln κ ln κ
Además, muestre que la velocidad promedio es
r=R
1 − κ2
2v0 r v0
Z
hvz i = 2 ln rdr = − 2κ2
R (1 − κ2 ) ln κ R 2 (1 − κ2 ) ln 1/κ
r=κR
99
(c) Show that the result in (b) simplifies to Eq. 2.2-26 if the film thickness is very small.
2B.8 Annular flow with inner cylinder moving axially. A cylindrical rod of radius nR moves axi-
ally with velocity v0 along the axis of a cylindrical cavity of radius R as seen in Fig. 2B.8. The
pressure at both ends of the cavity is the same, so that the fluid moves through the annular
region solely because of the rod motion.
Notas de clase: Francisco J. Valdés Parada
Cylinder of inside
Fluid at modified radius R Fluid at modified
pressure 0 pressure 0
Rod of radius κR
L
moving with velocity v0
Fig. 2B.8 Annular flow with the inner cylinder moving axially.
Ejercicio 4.15.
Un viscosı́metro de cilindro descendente consiste en un tubo cerrado (de radio R) por sus dos
(a) Find the velocity distribution in the narrow annular region.
extremos que contiene un cilindro interno (de radio κR) que se deja caer (ver figura). Suponiendo
(b) de
que la velocidad Find the mass
la caı́da del rate of flow
cilindro through
interno the annular
es constante region.
y vale v0 . Dirigiendo la atención al
espacio comprendido entre
(c) Obtain theel viscous
cilindroforce
interno y elon
acting externo,
the roddemuestre lo siguiente:
over the length L.
(d) Show
1. Que el perfil that the está
de velocidad resultdado
in (c)por:
can be written as a “plane slit” formula multiplied by a “curvature
correction.” Problems of this kind arise in studying the performance of plastics extrusion dies
1
for wire coating. R2 dP r2 R2 dP
ln r/R
vz = − 1 − 2 − v0 − 1 − κ2
v 4µ
(r) dz ln(r∕R)R 4µ dz ln κ
ln n
Answers: (a) z =
v0
0R2 iv0
2. Que la velocidad promedio es [ ]
1 −n2
− 2n2
2 ln (1∕n) R dP
(b)R2wdP= 2
κ2 1
hvz i = − 2
1 − κ + v0 − 1−κ 2
+
(c) Fz = −204Lv0 ∕ ln(1∕n)
8µ dz 4µ dz (1 − κ2 ) ln κ2
204Lv0
( ) del cilindro interno (πκ2 R2 v0 )
s s
3. Suponiendo que el fluido es incompresible, 1el flujo
1 volumétrico
s
(d) F = −
puede igualarse con el z flujo volumétrico 1 − − 2
+ · · ·
2en el12espacio entre los dos cilindros. Utilice esta
where s = 1 − n (see Problem 2B.6)
observación para deducir que el gradiente de presión está dado por
R2 dP v0
=
2B.9 Analysis of a capillary 4µflow dz meter. 2
(1 −Determine
κ ) + (1 +the ln κ of flow (in lbm ∕hr) through the capillary
κ2 )rate
flow meter shown in Fig. 2B.9. The fluid flowing in the inclined tube is water at 20∘ C, and
4. Utilice este
the resultado
manometer en fluid
el resultado
is carbon del tetrachloride
punto 1 para (CCl re-escribir el perfil de velocidad 3como
4 ) with density 1.594 g∕cm . The capillary
diameter is 0.010 in. Note: Measurements of Hand L are sufficient to calculate the flow rate; p
r2
need not be measured. zWhy?v 1 − R 2 + 1 + κ2 ln r/R
=−
v0 (1 − κ2 ) + (1 + κ2 ) ln κ
1 J.
B. Paton, P. H. Squires, W. H. Darnell, F. M. Cash, and J. F. Carley, Processing of Thermoplastic
Materials, E. C. Bernhardt (ed.), Reinhold, New York (1959), Chapter 4.
100
in which y = r∕R and y is a dummy variable of integration. The authors give no derivation, nor
do they give a literature citation. As the referee you feel obliged to derive the formula and list
any restrictions implied.
R vz(r)
κR
Container
Falling wall
cylinder
–v0
Cylindrical H
slug descends
with speed v0
Cylindrical container κR R
filled with fluid
z (b)
y
x
(a)
90 Chapter 3 TheFig.
Equations of Change for Isothermal Systems
2C.4 (a) A falling-cylinder viscometer with a tightly fitting solid cylinder moving
vertically. The cylinder is usually equipped with fins to maintain centering within the tube.
The fluid completely fills the tube and the top and bottom are closed. (b) Magnification of the
-1
gap between the slug and the container wall, illustrating the velocity profile vz (r).
Torsion wire with torsion constant k,
As the cylindrical slug descends in the −z direction, the (incompressible) fluid is pushed
ve is a function
upward through r
the gapofbetween the falling slug and the container walls. A sketch of the
resulting fluid velocity profile is illustrated in Fig. 2C.4(b). Note that a portion of the fluid will
be moving upward, while another portion of the fluid is moving downward. Note also that
the pressure at the bottom of the falling slug (p0 ) will be greater than the pressure at the top
6 J. Lohrenz, G. W. Swift, and F. Kurata, AIChE Journal, 6, 547–550 (1960) and 7, 6S (1961); E. Ashare,
R. B. Bird, and J. A. Lescarboura, AIChE Journal, 11, 910–916 (1965); F. J. Eichstadt and G. W. Swift, AIChE
Journal, 12, 1179–1183 (1966); M. C. S. Chen and J. A. Lescarboura, AIChE Journal, 14, 123–127 (1968).
Ir.">-
(a)
Fluid inside
Fig. 3.6-1. (a)Tangential laminar flow of an incompressible fluid in the space between two cylinders; the outer
Figura
one is moving with an angular velocity In,4.7:
.(b) AViscosı́metro de Couette
diagram of a Couette viscometer. One measures the angular
velocity Ino of the cup and the deflection oBof the bob at steady-state operation. Equation 3.6-31 gives the vis-
cosity p in terms of a, and the torque T, = kt&.
SOLUTION In the portion of the annulus under101consideration the fluid moves in a circular pattern. Rea-
sonable postulates for the velocity and pressure are: v, = v,(r), v, = 0, v, = 0, and p = p(r, 2).
We expect p to depend on z because of gravity and on r because of the centrifugal force.
For these postulates all the terms in the equation of continuity are zero, and the compo-
nents of the equation of motion simplify to
Notas de clase: Francisco J. Valdés Parada
la siguiente ecuación:
1 ∂ 2 vθ ∂ 2 vθ
∂vθ vθ ∂vθ ∂vθ vr vθ 1 ∂p ∂ 1 ∂(rvθ ) 2 ∂vr
ρ vr + + vz + =− + ρgθ + µ + 2 + +
∂r r ∂θ ∂z r r ∂θ ∂r r ∂r r ∂θ2 ∂z 2 r2 ∂θ
(4.107)
Suponiendo flujo desarrollado, tenemos que vr = vz = 0, además sabemos que gθ = 0 y que no los
cambios de presión se dan principalmente en las direcciones r y z, la ecuación anterior se reduce a:
1 ∂ 2 vθ ∂ 2 vθ
∂ 1 ∂(rvθ )
+ 2 + =0 (4.108)
∂r r ∂r r ∂θ2 ∂z 2
Para simplificar esta ecuación usamos los siguientes estimados de orden de magnitud:
∂ 2 vθ v
θ
2
= O 2
(4.109a)
∂z L
∂ 1 ∂(rvθ ) vθ
=O (4.109b)
∂r r ∂r R2 (1 − κ2 )
Suponiendo que:
R2 (1 − κ2 ) L2 (4.110)
entonces,
1 ∂ 2 vθ
∂ 1 ∂(rvθ )
+ =0 (4.111)
∂r r ∂r r2 ∂θ2
Como última simplificación, notemos que, bajo la suposición de flujo axisimétrico, la ecuación
anterior se reduce a:
d 1 d(rvθ )
=0 (4.112a)
dr r dr
Para resolver esta ecuación imponemos las siguientes condiciones de frontera:
En r = κR, vθ = 0 (4.112b)
En r = R, vθ = Ω0 R (4.112c)
102
Notas de clase: Francisco J. Valdés Parada
La fuerza que ejerce el fluido sobre la pared interna es:
Z
fθ = n · T · eθ dA (4.115)
A
Ejercicio 4.17.
Resuelva de nuevo el problema pero ahora considere el caso en el que el cilindro interno gira con
una velocidad angular Ωi y el cilindro exterior se manteiene fijo. Esto es:
En r = κR, vθ = Ωi κR
En r = R, vθ = 0
Muestre que, en este caso, el perfil de velocidad angular está dado por
R r
−
r R
vθ = Ωi κR
1
−κ
κ
103
Notas de clase: Francisco J. Valdés Parada
104
Notas de clase: Francisco J. Valdés Parada
Contenido
5.1. Herramental matemático . . . . . . . . . . . . . . . . . . . . . . . . . 106
5.1.1. Transformada de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . 106
5.1.2. Acerca de la solución de ecuaciones diferenciales de segundo orden . . 107
5.1.3. Identidades de funciones trigonométricas hiperbólicas . . . . . . . . . 108
5.1.4. Solución de ecuaciones no homogéneas . . . . . . . . . . . . . . . . . . 109
5.1.5. Ecuación diferencial de Bessel . . . . . . . . . . . . . . . . . . . . . . . 109
5.2. Aceleración repentina del fluido comprendido entre dos placas . 111
5.2.1. Análisis a tiempos cortos . . . . . . . . . . . . . . . . . . . . . . . . . 111
5.3. Análisis en estado transitorio . . . . . . . . . . . . . . . . . . . . . . 114
5.4. Flujo transitorio en un tubo circular . . . . . . . . . . . . . . . . . . 117
5.4.1. El sistema y sus fronteras . . . . . . . . . . . . . . . . . . . . . . . . . 117
5.4.2. Solución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119
5.5. Uso del método de expansión en funciones propias para resolver
problemas transitorios . . . . . . . . . . . . . . . . . . . . . . . . . . 120
105
Notas de clase: Francisco J. Valdés Parada
5.1. Herramental matemático
5.1.1. Transformada de Laplace
En este capı́tulo abordaremos la solución de problemas en estado transitorio. Para resolver
estos problemas, una alternativa (de muchas) es utilizar el método de la transformada de Laplace.
En esta sección, revisaremos rápidamente los elementos esenciales de esta alternativa. El lector
experimentado en esta técnica puede saltarse este material.
El método de la transformada de Laplace consiste en la aplicación de un operador integral a la
ecuación diferencial que permite tratar la parte temporal de manera algebraica y en consecuencia,
la ecuación diferencial se vuelve ahora más sencilla para su solución. Una vez resuelta la ecuación
diferencial, se utilizan las condiciones de frontera (escritas también en el dominio de Laplace)
para encontrar la solución particular. Por último, se aplica el operador de transformada inversa de
Laplace para expresar a la solución del problema en el dominio del tiempo.
La transformada de Laplace de una función del tiempo f (t) se define como
Z∞
L {f (t)} = f (t)e−st dt = f (s) (5.1)
0
γ+iT
1
Z
−1
f (t) = L {f (s)} = lı́m esT f (s) ds (5.2)
2πi T →∞
γ−iT
donde γ es un número real que se elige para garantizar la convergencia de la integral en el dominio
del tiempo. Existen valores tabulados de las transformadas de Laplace de muchas funciones y sus
derivadas, por lo que raras veces es necesario llevar a cabo las operaciones de integración indicadas
en las definiciones del operador de transformada de Laplace y de su inversa. Para el problema que
deseamos resolver, podemos considerar la siguiente fórmula
∂f
L = sf (s) − f (0) (5.3)
∂t
donde f (0) se refiere a la condición inicial. Más aún, como el operador de transformada de Laplace
afecta a la dependencia temporal, es posible intercambiar diferenciación e integración directamente
cuando se aplica a una derivada espacial. Esto es,
∂2f ∂2f
L = (5.4)
∂x2 ∂x2
Otras fórmulas de transformada de Laplace útiles para los desarrollos que se presentan en este
capı́tulo son
106
Notas de clase: Francisco J. Valdés Parada
f (s) f (t)
a
a
s
a
at
s2
√
e−a s
a
erfc √
s 2 t
√ n
2 2
sinh( sx) x 2 P∞ (−1) nπx n π t
√ + n=1 sin exp − 2
s sinh( sa) a a n a a
√ ∞
J0 (ix s) 1 2 P J0 (λn x/a)
2
√ (x − a2 ) + t + 2a2 3
exp(−λ2n t/a2 )
s J0 (ia s) 4 n=1 λn J1 (λn )
donde a es una constante y en la última fórmula, λn son cada una de las raı́ces positivas de la
ecuación
J0 (λn ) = 0 (5.5)
Además erfc(x) es la función error complementaria, la cual se define como
ξ=x
2
Z
2
erfc(x) = 1 − erf(x) = 1 − √ e−ξ dξ (5.6)
π
ξ=0
ξ=∞ √
π
Z
2
e−ξ dξ = (5.7)
2
ξ=0
se tiene que erf(∞) = 1, por lo que es sencillo deducir que la función error es una función creciente
acotada entre cero y uno, mientras que la función error complementaria es una función decreciente
acotada también entre cero y uno.
d2 y
− αy = 0 (5.8)
dx2
donde α es una constante positiva conocida. Para resolver esta ecuación se propone la siguiente
solución (también conocida como función de prueba)
y = emx (5.9)
donde m es una constante que se debe determinar. Sustituyendo esta propuesta en la ecuación (5.8)
resulta que
m2 − α emx = 0
(5.10)
107
Notas de clase: Francisco J. Valdés Parada
Para satisfacer esta ecuación debe ser cierto que el término entre paréntesis es cero, ya que de otra
forma la función de prueba es cero y la solución del problema serı́a cero. Entonces, tenemos la
siguiente ecuación algebraica
m2 − α = 0 (5.11)
√
de donde se deduce que m = ± α. Dado que hay dos valores de m que satisfacen la ecuación, la
solución se expresa como una superposición lineal
√ √
αx
y = c1 e + c2 e− αx
(5.12)
donde c1 y c2 se determinan a partir de las condiciones de frontera. Esta forma de expresar la
solución es muy útil cuando se manejan problemas en dominios semi-infinitos. Cuando el dominio
de solución es finito, es conveniente recordar las definiciones del seno y el coseno hiperbólicos
ex − e−x ex + e−x
sinh x = ; cosh x = (5.13)
2 2
Note que al combinar estas definiciones se obtiene que
ex = sinh x + cosh x
(5.14)
e−x = cosh x − sinh x
Usando estos resultados en la ecuación (5.42) se obtiene que
√ √
y = A sinh αx + B cosh αx (5.15)
donde A = c1 − c2 y B = c1 + c2 son dos constantes a determinar mediante las condiciones de
frontera.
108
Notas de clase: Francisco J. Valdés Parada
5.1.4. Solución de ecuaciones no homogéneas
En algunos problemas, es normal encontrarnos con casos como el siguiente:
d2 y
− αy = f (x) (5.21)
dx2
Para resolver este problema, la solución se parte en dos: la solución del problema homogéneo (es
decir, tomando f (x) = 0) y la solución del problema particular. Es decir,
y = yh + yp (5.22)
√ √
Donde en este caso yh = A sinh( αx) + B cosh( αx), mientras que yp depende de la forma
matemática que tenga f (x). Por ejemplo, digamos que f (x) = ax, en este caso corresponde a
un polinomio de primer orden y se propone
yp = K0 + K1 x (5.23)
Donde Jn (x) es la función de Bessel de primera clase de orden n. Para el caso en el que n sea entero,
la solución general es √ √
y = AJn ( αx) + BYn ( αx), ∀n = 0, 1, 2, . . . (5.30)
109
Notas de clase: Francisco J. Valdés Parada
a) b)
c) d)
Figura 5.1: Evaluación de las funciones Bessel a) de primera y b) segunda clase y de las funciones
modificadas de Bessel c) de primera y d) segunda clase.
donde Yn (x) es la función de Bessel de segunda clase de orden n. La evaluación de estas funciones
se encuentra en la figura 5.1, donde puede apreciarse que son funciones oscilatorias decrecientes.
Ahora bien, en ocasiones se presenta la siguiente forma de la ecuación diferencial,
d2 y dy
x2 2
+x − (αx2 − n2 )y = 0, n≥0 (5.31)
dx dx
la cual se conoce como la ecuación diferencial de Bessel modificada. En este caso las soluciones
generales son √ √
y = AIn ( αx) + BI−n ( αx), ∀n 6= 0, 1, 2, . . . (5.32a)
√ √
y = AIn ( αx) + BKn ( αx), ∀n = 0, 1, 2, . . . (5.32b)
donde In (x) y Kn (x) son las funciones de Bessel modificadas de primera y segunda clase,
respectivamente. La relación entre la función de Bessel modificada de primer orden y Jn (x) es
110
Notas de clase: Francisco J. Valdés Parada
6
y
y b 6
6
6
6 x
-
6
x vmax
a) b)
Figura 5.2: a) El sistema y sus fronteras en la configuración de referencia (t ≤ 0). El fluido está
en reposo y sus fronteras sólidas también. b) En un tiempo t > 0, la lámina sólida se mueve a
velocidad constante vmax . El fluido en contacto con ella se mueve también con la misma velocidad
vx (y = 0, t = 0) = vmax (condición de no deslizamiento).
∇·v=0 (5.35)
Suponiendo que, a pesar de lo transitorio del proceso, el flujo se puede suponer desarrollado y
por lo tanto vy = vz = 0. Bajo estas condiciones, la ecuación anterior se reduce a
∂vx
=0 (5.37)
∂x
111
Notas de clase: Francisco J. Valdés Parada
Dirigiendo la atención a la componente x de la ecuación de Navier-Stokes, obtenemos
∂vx ∂p ∂ 2 vx ∂ 2 vx
ρ =− +µ 2 +µ 2 (5.39)
∂t ∂x ∂y ∂z
Note que aquı́ además se ignoró la contribución del esfuerzo viscoso en la dirección x ya que la
velocidad es independiente de esta dirección. Más aún, en este problema el movimiento se debe
al desplazamiento de la placa y no por un gradiente de presión inducido. Por tanto, la ecuación
anterior se puede reducir a
∂vx ∂ 2 vx ∂ 2 vx
ρ =µ 2 +µ 2 (5.40)
∂t ∂y ∂z
Con el ánimo de simplificar más esta ecuación, consideremos los siguientes estimados de órdenes
de magnitud
∂ 2 vx v
x
µ 2 =O µ 2 (5.41a)
∂y b
∂ 2 vx v
x
µ = O µ 2 (5.41b)
∂z 2 W
donde W es la profundidad de las placas. Suponiendo que hay una separación de longitudes
caracterı́sticas del tipo b2 W 2 , la ecuación (5.40) se reduce a
∂vx ∂ 2 vx
ρ =µ 2 (5.42)
∂t ∂y
A pesar del progreso obtenido, el problema ya no puede simplificarse más y hemos terminado con
una ecuación diferencial parcial que requiere dos condiciones de frontera y una condición inicial.
Para tiempos posteriores a la condición inicial, se puede imponer la condición de no deslizamiento
en la placa inferior:
y = 0, vx = vmáx (5.43)
donde, por ahora, se supondrá que vmáx es una constante conocida. Dado que estamos interesados
en lo que ocurre a tiempos cortos, el dominio de solución del problema puede tratarse como un
dominio semi-infinito, en donde se puede suponer que el fluido que se encuentra suficientemente
lejos de la placa inferior aún no ha percibido este movimiento. Bajo esta suposición, es razonable
imponer la siguiente condición de frontera
y → ∞, vx = 0 (5.44)
Cuando t = 0, vx = 0 (5.45)
112
Notas de clase: Francisco J. Valdés Parada
Con esta información disponible, podemos proceder a aplicar el operador de transformada de
Laplace a la ecuación diferencial (5.42); el resultado es
∂ 2vx ρ
2
− sv x = 0 (5.46)
∂y µ
donde ya se tomó en cuenta a la condición inicial. Las condiciones de frontera pueden escribirse en
el dominio de Laplace como
vmáx
y = 0, v x = (5.47)
s
y → ∞, v x = 0 (5.48)
Dado que el dominio de solución es semi-infinito, es conveniente expresar la solución de esta
ecuación en términos de funciones exponenciales como sigue
r r
ρ ρ
v x = c1 exp sy + c2 exp − sy (5.49)
µ µ
Ejercicio 5.1.
Deduzca el perfil de velocidad dado en la ecuación (5.52).
Ejercicio 5.2.
La definición de una función error es:
u=x
2
Z
2
erf(x) = √ e−u du
π
u=0
113
Notas de clase: Francisco J. Valdés Parada
1
0.8
0.6
vx
vmáx νt
0.4
b2 = 10−4
10−3
0.2 10−2
10−1
∞
0
10−4 10−3 10−2 10−1 100
y/b
Figura 5.3: Dinámica de los perfiles de velocidad obtenidos a partir de la evaluación de la ecuación
vx
(5.51). El perfil en estado estacionario (lı́nea roja) resulta de usar la expresión vmáx = 1 − yb .
Derive ambos lados de la ecuación anterior respecto a x y utilice la regla de Leibniz para probar
que
derf(x) 2 2
= √ e−x
dx π
Utilice este resultado para probar que
2
e−x
Z
erf(x)dx = xerf(x) + √
π
Tomando en cuenta lo anterior, calcule el valor promedio de la velocidad
y=b
1
Z
hvx i = vx dx
b
y=0
114
Notas de clase: Francisco J. Valdés Parada
Como se mostró, esta suposición puede fallar conforme aumenta el tiempo, lo que motiva a revisar
el planteamiento del problema para que la solución no esté limitada por la suposición de medio
semi-infinito.
La versión simplificada de la ecuación diferencial de conservación de cantidad de movimiento
sigue siendo válida, aquı́ la re-escribimos en términos de la viscosidad cinemática como
∂vx ∂ 2 vx
=ν (5.54)
∂t ∂y 2
La cual ahora está sujeta a las siguientes condiciones de frontera e inicial en un dominio finito
y = 0, vx = vmáx (5.55)
y = b, vx = 0 (5.56)
Para resolver este problema recurrimos de nuevo al método de la transformada de Laplace, por lo
que aplicando este operador al problema en el dominio del tiempo obtenemos el siguiente problema
en el dominio de Laplace
∂ 2 v̄x s
2
− v̄x = 0 (5.57a)
∂y ν
vmáx
y = 0, v̄x = (5.57b)
s
y = b, v̄x = 0 (5.57c)
Como el dominio es finito, es conveniente expresar solución general en términos de funciones
trigonométricas hiperbólicas como sigue
r r
s s
v̄x = A sinh y + B cosh y (5.58)
ν ν
115
Notas de clase: Francisco J. Valdés Parada
1
0.8
0.6
vx
vmáx νt
0.4 b2 = 10−4
10−3
10−2
0.2
10−1
100
0
10−4 10−3 10−2 10−1 100
y/b
Figura 5.4: Comparación de la dinámica de los perfiles de velocidad predichos con la solución
transitoria dada en la ecuación (5.61) (lı́neas) y aquellos resultantes de usar la suposición de medio
semiı́nfinito (cı́rculos, excepto los cı́rculos rojos). Los cı́rculos rojos corresponden a la solución en
estado estacionario.
116
Notas de clase: Francisco J. Valdés Parada
0.5
0.4 exacta
aproximada
hvx i 0.3 estacionaria
vmax
0.2
0.1
0
10−5 10−4 10−3 10−2 10−1 100 101
τ
Figura 5.5: Comparación de la dinámica la velocidad promedio adimensional hvx i/vmax obtenida
con las soluciones exacta (Ec. 5.63), aproximada (Ec. 5.53) y en estado estacionario. Se utilizaron
1000 términos en la sumatoria. En el eje horizontal τ = tν/b2 .
117
Notas de clase: Francisco J. Valdés Parada
En la Figura 5.6 b) se muestra la condición inicial de flujo, que consiste en el cambio de la presión,
donde, para el caso mostrado en la figura, p0 > pa y el flujo procede en la dirección descendente.
Debido a la condición de no-deslizamiento, el fluido que cae escurre con mayor facilidad por el
centro del tubo que por sus paredes, a las cuales tiende a quedar adherido.
Las suposiciones para plantear el modelo son las mismas que en el caso anterior, excepto que no
se pueden despreciar los efectos de la presión y la gravedad sobre el flujo. Sin embargo, si el flujo
es suficientemente lento, puede seguirse suponiendo que el tiempo caracterı́stico de la convección
de momento es mucho mayor que el tiempo caracterı́stico del proceso y el tiempo caracterı́stico de
la difusión de momento. De esta forma, la ecuación de Navier-Stokes se reduce a:
∂v
ρ = −∇p + ρg + µ∇2 v (5.66)
∂t
Como las fuerzas de presión y de gravedad actúan en la dirección z, es conveniente extraer la
componente en dicha dirección de la ecuación anterior, es decir:
∂vz ∂p
ρ =− + ρg + µ∇2 vz (5.67)
∂t ∂z
donde el operador Laplaciano en coordenadas cilı́ndricas es:
1 ∂ 2 vz ∂ 2 vz
2 1 ∂ ∂vz
∇ vz = r + 2 + (5.68)
r ∂r ∂r r ∂θ2 ∂z 2
Suponiendo flujo axisimétrico, podemos ignorar la dependencia angular de la velocidad, por lo que
usando la definición de la presión modificada tenemos
∂ 2 vz
∂vz ∂P µ ∂ ∂vz
ρ =− + r +µ 2 (5.69)
∂t ∂z r ∂r ∂r ∂z
P = p0 P = p0
HH HH
jr
H jr
H
z ? z ?
HH
P = p0
z=L P = pa
r=a
a) b)
Figura 5.6: Esquemas del sistema en: a) su configuración de referencia y b) la condición inicial
donde la presión a la salida cambió de p0 a la presión atmosférica pa .
118
Notas de clase: Francisco J. Valdés Parada
Para simplificar esta ecuación usamos los siguientes estimados de órdenes de magnitud
µ ∂ ∂vz µv
z
r =O (5.70a)
r ∂r ∂r R2
∂ 2 vz µv
z
µ 2
= O 2
(5.70b)
∂z L
De esta forma, suponiendo que R2 L2 , podemos despreciar la contribución axial de los esfuerzos
cortantes y la ecuación de movimiento se reduce a
∂vz ∂P µ ∂ ∂vz
ρ =− + r (5.71)
∂t ∂z r ∂r ∂r
donde ∂P
∂z es una constante. Para completar el planteamiento del problema tenemos las siguientes
condiciones de frontera
vz está definida ∀r, ∀t (5.72)
en r = R, vz = 0, ∀t (5.73)
Y la condición inicial corresponde a condiciones hidrostáticas
cuando t = 0, vz = 0 (5.74)
5.4.2. Solución
Para resolver el problema recurrimos de nuevo al método de la transforma de Laplace. Aplicando
este operador a la ecuación diferencial se obtiene
1 ∂ ∂v̄z ρ 1 ∂P 1
r − sv̄z = (5.75)
r ∂r ∂r µ µ ∂z s
Dado que el dominio de solución incluye al punto r = 0 y en ese punto K0 (x) tiende a menos
infinito, entonces debe cumplirse que B = 0. Usando la condición de frontera en r = R, v̄z = 0,
resulta que √
r s
1 ∂P I0 √ν 1 ∂P 1
v̄z = √ − (5.77)
ρ ∂z s2 I R√ s ρ ∂z s2
0 ν
119
Notas de clase: Francisco J. Valdés Parada
Aplicando el operador de transformada inversa, se obtiene la siguiente expresión
∞
" #
R2 ∂P r2
X J0 (λn r/R) 2 νt
vz = − 1− 2 −8 exp −λn 2 (5.79)
4µ ∂z R λ3n J1 (λn ) R
n=1
R2
Note que − 4µ ∂P
∂z corresponde al valor máximo de la velocidad en estado estacionario, denotando
dicho valor como vmáx , obtenemos
∞
r2
vz X J0 (λn r/R) 2 νt
=1− 2 −8 exp −λn 2 (5.80)
vmáx R λ3n J1 (λn ) R
n=1
En la figura 5.7 se muestra la evaluación de esta expresión. Como puede notarse, los efectos
transitorios toman lugar más claramente en el intervalo de tiempo adimensional (τ = νt/R2 )
comprendido entre 10−2 y 100 . Más aún, de la evaluación de la dinámica de los perfiles de velocidad
transitorios se puede notar la naturaleza parabólica a distintos tiempos, cuyo lado recto aumenta
con el tiempo.
Ejercicio 5.3.
Demuestre que el esfuerzo cortante en las paredes del tubo está dado por la siguiente expresión
∞
" #
∂vz 2µvmáx X 1 νt
τrz |r=R = µ = −1 + 4 exp −λ2n 2
∂r r=R R λ2n R
n=1
Evalúe esta expresión para determinar si el esfuerzo aumenta o disminuye con el tiempo.
120
Notas de clase: Francisco J. Valdés Parada
0.5 νt
= 10−2
b2
0.05
r 0.10
0
R 0.20
∞
−0.5
−1
vx /vmáx
a)
0.5
0.4
transitoria
estacionaria
hvz i 0.3
vref
0.2
0.1
0
10−5 10−4 10−3 10−2 10−1 100 101
τ
b)
Figura 5.7: Dinámica de la velocidad a) puntual y b) promedio adimensional para el flujo en una
tuberı́a. Los perfiles transitorios se obtuvieron utilizando más de 30000 términos en la sumatoria.
En b) la velocidad de referencia se tomó como la velocidad máxima en estado estacionario; además,
τ = νt/R2 .
121
Notas de clase: Francisco J. Valdés Parada
Para ello, retomemos el último problema para la situación más general en que el gradiente de
presión puede ser función conocida del tiempo. De esta forma el problema es
∂vz ∂P µ ∂ ∂vz
ρ =− + r (5.84a)
∂t ∂z r ∂r ∂r
vz está definida ∀r, ∀t (5.84b)
en r = R, vz = 0, ∀t (5.84c)
cuando t = 0, vz = 0, ∀r (5.84d)
Claramente el gradiente de presión es la fuente del problema. Para avanzar hacia la solución, es
conveniente considerar la versión homogénea del problema, esto es
ρ ∂vz,H 1 ∂ ∂vz,H
= r (5.85a)
µ ∂t r ∂r ∂r
vz,H está definida ∀r, ∀t (5.85b)
en r = R, vz,H = 0, ∀t (5.85c)
cuando t = 0, vz,H = 0, ∀r (5.85d)
Dirigiendo la atención a la ecuación diferencial, se puede proponer una solución del tipo:
De acuerdo a lo anterior, se deduce que, en este caso, el signo de la constante debe ser negativo
y por tanto
ρ 1 dF 1 d dG
= r = −λ2 (5.90)
µ F dt Gr dr dr
122
Notas de clase: Francisco J. Valdés Parada
De donde se obtienen las siguientes ecuaciones diferenciales ordinarias
dF 2 d dG
= −νλ F ; r + λ2 Gr = 0 (5.91)
dt dr dr
Dirigiendo la atención a la ecuación diferencial de G, notamos que su estructura es de tipo Sturm-
Liouville y corresponde a la ecuación diferencial de Bessel. La solución general de este problema
es
G = AJ0 (λn r) + BY0 (λn r) (5.92)
Dado que el dominio incluye a r = 0 y Y0 (0) → −∞, se deduce que B = 0. Por otro lado, la
condición de frontera en r = R lleva a que
Si la constante A fuera cero, entonces G se reduce a cero y la solución del problema serı́a la trivial.
Como este no es el caso, debe cumplirse que
J0 (λn R) = 0 (5.94)
Dado que la función Bessel de primera clase es una función oscilatoria, existe una infinidad de
soluciones a esta ecuación. Por lo tanto, la solución del problema puede expresarse como sigue
∞
X
G= An J0 (λn r) (5.95)
n=1
donde J0 (αn r/R) es la función propia y αn son los valores propios que resultan de resolver la
ecuación J0 (αn ) = 0. Ahora bien, tomando en cuenta la solución propuesta en la ecuación (5.86),
entonces
∞
X
vz (r, t) = An (t)φn (r) (5.97)
n=1
donde φn (r) = J0 (αn r/R). Note que, debido a la condición inicial, An (0) = 0.
Ahora bien, para desarrollar la solución del problema nos ayudaremos de la fórmula de Green,
Z Z
2 2
φ∇ ϕ − ϕ∇ φ dV = n · (φ∇ϕ − ϕ∇φ) dA (5.98)
V A
123
Notas de clase: Francisco J. Valdés Parada
Note que, debido a las condiciones de frontera para vz y φm , el lado derecho de la ecuación anterior es
cero. Más aún, usando las ecuaciones diferenciales correspondientes para las variables dependientes,
la ecuación anterior se reduce a
r=R
Z
1 ∂vz 1 ∂P
φm + + λ2m φm vz rdr = 0 (5.100)
ν ∂t µ ∂z
r=0
Ejercicio 5.4.
Muestre que, en el caso en el que el gradiente de presión sea constante, el resultado dado en la
ecuación (5.105) se reduce a la ecuación (5.80).
Ejercicio 5.5.
Reconsidere el problema visto en esta sección para el caso en que la condición inicial sea ahora
r2
cuando t = 0, vz = −vmáx 1 − 2
R
Fı́sicamente, esta condición corresponde al caso en el que el flujo en el tubo halla legado a un estado
estacionario y que, en un tiempo t > 0, se invierta la dirección del gradiente de presión.
124
Notas de clase: Francisco J. Valdés Parada
Contenido
6.1. Solución aproximada de la placa que se acelera súbitamente . . . 126
6.1.1. Revisión del problema y la solución en un medio semi-infinito . . . . . 126
6.1.2. Método aproximado de von Kármán-Pohlhausen . . . . . . . . . . . . 126
6.2. Análisis cerca de la entrada: Planteamiento del problema . . . . . 130
6.2.1. Análisis en la zona de velocidad máxima . . . . . . . . . . . . . . . . . 130
6.2.2. Análisis en la zona cercana a las paredes fijas . . . . . . . . . . . . . . 130
125
Notas de clase: Francisco J. Valdés Parada
6.1. Solución aproximada de la placa que se acelera súbitamente
6.1.1. Revisión del problema y la solución en un medio semi-infinito
En el capı́tulo anterior resolvimos el problema de un fluido comprendido entre dos placas que
estaba inicialmente en reposo y posteriormente la placa inferior se mueve a una velocidad constante
vmáx . Después de hacer las simplificaciones correspondientes el problema se redujo al siguiente
∂vx ∂ 2 vx
=ν (6.1a)
∂t ∂y 2
En y = 0, vx = vmáx (6.1b)
En y → ∞, vx = 0 (6.1c)
Cuando t = 0, vx = 0 (6.1d)
Usando el método de la transformada de Laplace se encontró el siguiente resultado
vx y
= erfc √ (6.2)
vmáx 2 νt
Como puede notarse, el lado izquierdo de la ecuación es adimensional y sus valores están acotados
entre 0 y 1. Además, el argumento de la función error complementaria es también adimensional. En
adelante nos referiremos a la expresión dada en la ecuación (6.2) como la solución analı́tica exacta
a pesar de que está limitada por la suposición
√de medio semi-infinito. La evaluación de esta solución
en términos de la variable adimensional y/2 νt se presenta en la figura 6.1. En esta figura además
se ubica la zona en donde la velocidad experimentó sus cambios más importantes. En este contexto,
entenderemos por cambios importantes a la zona en la que ocurren el 99 % de los cambios. Es a
esta región de cambios a la que denominaremos capa lı́mite. Para este caso en particular notamos
que la capa lı́mite corresponde
√ a la zona en la que la velocidad adimensional disminuyó de 1 hasta
0.01. El valor de y/2 νt donde termina la capa lı́mite es aproximadamente 1.8, esto es,
√
δ = 3.6 νt (6.3)
donde δ representa el espesor de la capa lı́mite y es, como se esperaba, una función del tiempo.
A continuación explicaremos un método aproximado para llegar a este resultado que nos será de
mucha utilidad en problemas posteriores.
126
Notas de clase: Francisco J. Valdés Parada
100
10−1 6
vx Capa lı́mite
vmáx
10−2
10−3
En y = δ, vx ≈ 0 (6.5b)
Para encontrar las otras dos condiciones de frontera podemos tomar en cuenta que después de la
capa lı́mite ya no se transporta el esfuerzo cortante, por lo tanto es razonable imponer que
∂vx
En y = δ, =0 (6.5c)
∂y
Por último, en y = 0 la velocidad es constante, esto es, no cambia con el tiempo, por lo que la
ecuación diferencial (6.1a) se reduce en esta frontera a
∂ 2 vx
En y = 0, =0 (6.5d)
∂y 2
Usando las condiciones de frontera se deduce que a = 1, b = −3/2, c = 0 y d = 1/2 y por tanto la
propuesta de solución es
y 3
vx 3 y 1
=1− + (6.6)
vmáx 2 δ(t) 2 δ(t)
A pesar del progreso obtenido, la solución aún no es utilizable pues está expresada en términos
de δ, la cual es una incógnita. Para atender esta situación, es necesario desarrollar y resolver una
ecuación para la dinámica del espesor de la capa lı́mite. Con este objetivo en mente, se integra
definidamente la ecuación de movimiento desde y = 0 hasta y = δ obtenemos
y=δ
∂vx y=δ
∂vx ∂vx
Z
dy = ν = −ν (6.7)
∂t ∂y y=0 ∂y y=0
y=0
127
Notas de clase: Francisco J. Valdés Parada
A partir de la ecuación (6.6), se deduce que
∂vx 3vmáx
=− (6.8a)
∂y y=0
2δ(t)
" #
y 3 dδ
∂vx 3vmáx y
= − (6.8b)
∂t 2δ (t) δ(t) δ(t) dt
Sustituyendo estos resultados en la ecuación (6.7) obtenemos
y=δ" #
y 3
y dδ
Z
− dy =ν (6.9)
δ(t) δ(t) dt
y=0
128
Notas de clase: Francisco J. Valdés Parada
Solución exacta
Solución aproximada
100
10−1
vx 10−2
vmáx
10−3
10−4
10−5
0 0.5 1 1.5 2 2.5
y
√
2 νt
Figura 6.2: Comparación de los perfiles de velocidad obtenidos con las soluciones analı́ticas exacta
y aproximada.
129
r v dV r b dV
d
Euler II: r t (n ) dA (3-3)
dt
Vm (t ) Vm (t ) Am ( t )
These three axioms, along with the laws of thermodynamics (see Chapter 11), provide the framework for
our analysis of fluid flow phenomena.
Le no-slip
no-slip
2. Viscosidad constante.
ρv · ∇v = −∇p (6.16)
cuya componente x es
du ∂p
ρU 2 u =− (6.17)
dx ∂x
donde u = u∗ /U .
130
Notas de clase: Francisco J. Valdés Parada
∂ 2 ux ∂ 2 ux
∂ux ∂vx 1 ∂p
ρU ux + uy =− +µ + µ (6.18b)
∂x ∂y U ∂x ∂x2 ∂y 2
∂ 2 uy ∂ 2 uy
∂uy ∂uy 1 ∂p
ρU ux + uy =− +µ 2 +µ 2 (6.18c)
∂x ∂y U ∂y ∂x ∂y
∂ux ∂uy
+ =0 (6.20a)
∂x ∂y
∂ 2 ux
∂ux ∂ux 1 ∂p
ρU ux + uy =− +µ 2 (6.20b)
∂x ∂y U ∂x ∂y
∂ 2 uy
∂uy ∂uy 1 ∂p
ρU ux + uy =− +µ 2 (6.20c)
∂x ∂y U ∂y ∂y
De las últimas dos ecuaciones podemos deducir los siguientes estimados de órdenes de magnitud:
1 L ρU 2
∂p
=O 1+ (6.21a)
∂x Reδ δ L
∂p ∂p δ
=O (6.21b)
∂y ∂x L
donde Reδ = ρU δ/µ. De esta forma, dado que, en general la presión es función de x y de y, tenemos
δ2
∂p ∂p ∂p
dp = dx + dy = O L 1+ 2 (6.22)
∂x ∂y ∂x L
∂p
dp ≈ dx (6.23)
∂x
Esto es, la presión no depende de y y su derivada respecto a x puede calcularse mediante la Ec.
(6.17). Lo anterior hace innecesaria la Ec. (6.20c) y nuestro problema es ahora:
ξ=y
∂ux
Z
uy = − dξ (6.24a)
∂x
ξ=0
131
Notasde clase: Francisco J. Valdés Parada
∂ 2 ux
∂ux ∂ux du
ρU ux + uy −u =µ 2 (6.24b)
∂x ∂y dx ∂y
O bien,
ξ=Y
∂ux ∂ux ∂ux du ∂ 2 ux
Z
ReH ux − dξ − u = (6.25)
∂X ∂Y ∂X dX ∂Y 2
ξ=0
donde ReH = HρU/µ, Y = y/H y X = x/H. Integrando definidamente esta ecuación desde Y = 0
hasta Y = ∆ = δ/H, resulta:
YZ=∆ YZ=∆ ξ=Y YZ=∆
∂ux ∂ux ∂ux du ∂ux
Z
ReH ux dY − dξdY − u dY = − (6.26)
∂X ∂Y ∂X dX ∂Y Y =0
Y =0 Y =0 ξ=0 Y =0
Aquı́ se tomó en cuenta el hecho que en Y = ∆, ∂ux /∂Y = 0. Al usar la fórmula de integración
por partes en el segundo término del lado izquierdo de la ecuación anterior se tiene:
Y =∆
YZ=∆ ξ=Y ξ=Y YZ=∆ YZ=∆
∂ux ∂ux ∂ux ∂ux ∂ux
Z Z
dξdY = ux (Y ) dξ − ux dY = (u − ux ) dY (6.27)
∂Y ∂X ∂X ∂X ∂X
Y =0 ξ=0 ξ=0 Y =0 Y =0
Y =0
132
Notas de clase: Francisco J. Valdés Parada
La solución de esta ecuación no es posible por métodos analı́ticos exactos, por ello recurrimos a un
método aproximado que se basa en notar que el perfil de velocidad en la zona de cambios cerca de
la pared es parabólico. Por ello, se propone la siguiente solución:
2
Y Y
ux (X, Y ) = a0 (X) + a1 (X) + a2 (X) (6.33)
∆ ∆
Esta ecuación debe satisfacer al menos las siguientes condiciones de frontera:
Y = 0, ux = 0 (6.34a)
Y = ∆, ux = u(x) (6.34b)
∂ux
Y = ∆, =0 (6.34c)
∂Y
Al utilizar estas condiciones resulta que a0 = 0, a1 = 2u y a2 = −u; de esta forma la Ec. (6.33)
queda como: " 2 #
Y Y
ux (X, Y ) = u(X) 2 − (6.35)
∆ ∆
Utilizando esta expresión se obtienen los siguientes resultados:
∂ux 2u
= (6.36a)
∂Y Y =0
∆
YZ=∆
2 2
[ux (u − ux )]dY = u ∆ (6.36b)
15
Y =0
YZ=∆
2u∆
ux dY = (6.36c)
3
Y =0
YZ=∆
u∆
(u − ux ) dY = (6.36d)
3
Y =0
−1 3 2 du u∆ d∆
ReH = ∆ + (6.39)
10 dX 15 dX
La cual tiene dos incógnitas u y ∆, para eliminar la primera recurrimos al principio de conservación
de masa, el cual requiere que:
A1 hvx i1 = A2 hvx i2 (6.40)
133
Notas de clase: Francisco J. Valdés Parada
donde
1
Z
hvx ii = vx dA (6.41)
Ai
Ai
Tomando como punto 1 la entrada del sistema (x = 0, donde vx = U ) y como punto 2 cualquier
otro punto en la dirección x, tenemos:
Z y=2H Z z=W Z y=2H Z z=W
U dzdy = vx dzdy (6.42)
y=0 z=0 y=0 z=0
en la última igualdad se tomó en cuenta la simetrı́a del problema. La ecuación anterior puede
expresarse como:
YZ=1 YZ=∆ YZ=1
∆
1= ux dY = ux dY + ux dY = u(x) 1 − (6.44)
|{z} 3
Y =0 Y =0 Y =∆ u(x)
En X = 0, ∆=0 (6.47)
En X = Le /H, ∆=1 (6.48)
1
Se recomienda hacer el cambio de variable ξ = 1 − ∆/3
134
Notas de clase: Francisco J. Valdés Parada
7 | Flujo reptante
Contenido
7.1. Flujo reptante alrededor de una esfera . . . . . . . . . . . . . . . . . 135
∇·v=0 (7.1a)
0 = −∇p + ρg + µ∇2 v (7.1b)
O bien, en términos de ∇P = ∇p − ρg:
∇2 ω = 0 (7.3)
donde ω es el vector de vorticidad, el cual se define como ω = ∇×v. Además, el vector de velocidad
puede expresarse en términos de la función de corriente, ψ de acuerdo a la siguiente expresión:
−1 ∂ψ
e3 1 ∂ψ
v = ∇ψ × = e1 + e2 (7.4)
h3 h2 h3 ∂u2 h1 h3 ∂u1
En coordenadas esféricas se tiene que h1 = 1, h2 = r y h3 = rsen(θ), por lo que la ecuación anterior
conduce a:
1 ∂ψ
vr = (7.5a)
r2 sen(θ)
∂θ
−1 ∂ψ
vθ = (7.5b)
rsen(θ) ∂r
Sustituyendo la definición de la velocidad en términos de la función de corriente en el vector de
vorticidad y sustituyendo la expresión resultante en la Ec. (7.3), se obtiene, en general,
E4ψ = 0 (7.6)
135
Notas de clase: Francisco J. Valdés Parada
136
Notas de clase: Francisco J. Valdés Parada
Tomando en cuenta las ecs. (7.5), se deduce que
∂ψ ∂ψ
r = R; = 0; =0 (7.12a)
∂θ ∂r
Lo cual indica que la función de corriente en la superficie de la esfera es una constante, la cual se
escoge como cero, esto es:
r = R; ψ = 0 (7.12b)
Dirigiendo la atención al lı́mite cuando r → ∞, de acuerdo a las ecs. (7.5), se tiene:
∂ψ
r → ∞; = r2 v∞ sen(θ) cos(θ) (7.13a)
∂θ
∂ψ
r → ∞; = rv∞ sen2 (θ) (7.13b)
∂r
Al integrar cualquiera de las ecuaciones anteriores se llega al mismo resultado:
r2 v∞ sen2 (θ)
r → ∞; ψ= (7.14)
2
Entonces, el problema está definido por la ecuación diferencial de cuarto orden (7.6), la cual está
sujeta a las condiciones de frontera dadas por las ecs. (7.12) y (7.14). Esta última es la única fuente
del problema y exhibe a única dependencia con el ángulo, la cual es sen2 (θ), lo que sugiere proponer
la siguiente solución:
ψ = f (r)sen2 (θ) (7.15)
Note que
d2 f
2 2 f
E ψ = sen (θ) 2
−2 2 (7.16)
dr r
Entonces:
d2 f d2 d2 f
2
1 ∂sen2 (θ)
4 2 2 f 2 2f d f 2f sen(θ) ∂
E ψ=E sen (θ) 2
−2 2 = sen (θ) 2 − 2 + − 2
dr r dr dr2 r dr2 r r2 ∂θ sen(θ) ∂θ
| {z }
−2sen2 (θ)/r2
(7.17)
Tras hacer el álgebra correspondiente, resulta la siguiente ecuación diferencial ordinaria de cuarto
orden:
d4 f 4 d2 f 8 df 8
4
− 2 2 + 3 − f =0 (7.18)
dr r dr r dr r4
O bien,
d4 f d2 f df
r4 4 − 4r2 2 + 8r − 8f = 0 (7.19)
dr dr dr
La cual es una ecuación del tipo Euler-Cauchy y su solución puede encontrarse al usar la función
de prueba f = rn , la cual al sustituirla en la ecuación anterior lleva a:
137
Notas de clase: Francisco J. Valdés Parada
Al utilizar las condiciones de frontera no redundantes (r = R, f = 0; r = R, df /dr = 0 y r → ∞,
f = r2 v∞ /2) se deduce que A = v∞ R3 /4, B = −3v∞ R/4, C = v∞ /2 y D = 0. Note que la
condición de frontera en r → ∞ permite deducir los valores de C y D. De esta forma, la función
de corriente es:
v∞ r2 sen2 (θ) 3 R 1 R3
ψ= 1− + (7.22)
2 2r 2 r3
Podemos ahora sustituir este resultado en las ecs. (7.5) para obtener:
3 R 1 R3
1 ∂ψ
vr = 2 = v∞ cos(θ) 1 − + (7.23a)
r sen(θ) ∂θ 2r 2 r3
−1 ∂ψ 3 R 1 R3
vθ = = −v∞ sen(θ) 1 − − (7.23b)
rsen(θ) ∂r 4r 4 r3
El siguiente paso en el análisis consiste en calcular la presión, para ello nos valemos de las
componentes r y θ de la ecuación de Stokes:
µ ∂2 2
∂P µ ∂ ∂vr
= 2 2 (r vr ) + 2 sen(θ) (7.24a)
∂r r ∂r r sen(θ) ∂θ ∂θ
1 ∂P µ ∂ 2 ∂vθ µ ∂ 1 ∂ 2µ ∂vr
= 2 r + 2 (vθ sen(θ)) + 2 (7.24b)
r ∂θ r ∂r ∂r r ∂θ sen(θ) ∂θ r ∂θ
Al llevar a cabo las operaciones correspondientes se obtiene:
∂P R
= 3µv∞ cos(θ) 3 (7.25a)
∂r r
∂P 3 R
= µv∞ sen(θ) 2 (7.25b)
∂θ 2 r
Dado que la presión depende de r y θ, se tiene con ayuda de la regla de la cadena:
∂P ∂P
dP = dr + dθ (7.26)
∂r ∂θ
Esta ecuación debe integrarse semidefinidamente partiendo de una condición conocida, la cual es
la siguiente:
r → ∞, θ = π/2; P = p0 (7.27)
Integrando la Ec. (7.26) se obtiene:
UZ=P ξ=r ζ=θ
∂P ∂P
Z Z
dU = dξ + dζ (7.28)
∂ξ θ=π/2 ∂ζ r
U =p0 ξ=∞ ζ=π/2
| {z }
=0
Cuyo resultado es
3 R
p = p0 − ρgz − µv∞ cos(θ) 2 (7.29)
2 r
donde z = r cos(θ); por lo que
3 R
p = p0 − ρgr cos(θ) − µv∞ cos(θ) 2 (7.30)
2 r
138
Notas de clase: Francisco J. Valdés Parada
El vector de fuerza en la superficie de la esfera es:
φ=2π
Z θ=π
Z
f|r=R = n · T|r=R R2 sen(θ)dθdφ (7.32)
φ=0 θ=0
tomando el producto punto en ambos lados de la ecuación anterior con el vector unitario ez , da
lugar a:
φ=2π
Z θ=π
Z
fz |r=R = n · T · [cos(θ)er − sen(θ)eθ ]|r=R R2 sen(θ)dθdφ (7.33)
φ=0 θ=0
En la ecuación anterior se tomó en cuenta que ez = cos(θ)er − sen(θ)eθ : Dado que n = er , tenemos:
φ=2π
Z θ=π
Z
fz |r=R = [Trr cos(θ) − Trθ sen(θ)]|r=R R2 sen(θ)dθdφ (7.34)
φ=0 θ=0
donde
∂vr
Trr = − p + 2µ (7.35a)
∂r
∂ vθ µ ∂vr
Trθ = µr + (7.35b)
∂r r r ∂θ
3
− p|r=R = − p0 + ρgR + µv∞ cos(θ) (7.36a)
2R
R3
∂vr 3 R
= v∞ cos(θ) − 4 =0 (7.36b)
∂r r=R 2 r2 r r=R
∂ vθ
3
µr =− µv∞ sen(θ) (7.36c)
∂r r r=R 2R
3 R 1 R3
µ ∂vr µv∞ sen(θ)
=− 1− + =0 (7.36d)
r ∂θ r=R R 2r 2 r3 r=R
139
Notas de clase: Francisco J. Valdés Parada
De esta forma, la fuerza es:
θ=π
1
Z
fz |r=R = − p0 sen(θ) cos(θ)dθ
2πR2
θ=0
| {z }
=0
θ=π
3
Z
+ ρgR + µv∞ sen(θ) cos2 (θ)dθ (7.37)
2R
θ=0
| {z }
=2/3
θ=π
3
Z
+ µv∞ sen3 (θ)dθ
2R
θ=0
| {z }
4/3
Por lo tanto:
4πR3
fz |r=R = ρg + 6πRµv∞ (7.38)
3
El primer término representa la fuerza boyante que ejerce el fluido sobre la esfera y el segundo es
consecuencia del movimiento del fluido y su arrastre viscoso y se conoce como la ley de Stokes.
Para el caso en que la esfera caiga en un fluido estático con velocidad terminal vt , la fuerza total
es igual al peso de la esfera, es decir:
4πR3
fz |r=R = ρs g (7.39)
3
Por lo tanto
4πR3
(ρs − ρ)g = 6πRµvt (7.40)
3
O bien,
2R2
µ= (ρs − ρ)g (7.41)
9vt
La expresión anterior puede utilizarse para calcular la viscosidad de un fluido.
140
Notas de clase: Francisco J. Valdés Parada
8 | Flujo potencial
Contenido
8.1. Flujo no viscoso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
8.2. Flujo potencial alrededor de un cilindro . . . . . . . . . . . . . . . . 143
8.3. Flujo potencial alrededor de una esfera . . . . . . . . . . . . . . . . 145
∇2 φ = 0 (8.4)
∇ × ∇φ = 0 (8.5)
w=∇×v (8.6)
141
Notas de clase: Francisco J. Valdés Parada
La ecuación (8.5) indica que, en problemas de flujo potencial, se satisface que el flujo es además
irrotacional. Por lo tanto,
w=∇×v=0 (8.7)
Dado que la definición de la función de corriente es:
e3
v = ∇ψ × (8.8)
h3
se tiene:
e3
∇ × ∇ψ × =0 (8.9)
h3
la cual da lugar a:
E2ψ = 0 (8.10)
donde
∂2ψ ∂2ψ
+ Cartesiano (x, y)
∂x2 ∂y 2
1 ∂2ψ
1 ∂ ∂ψ
r + Cilı́ndrico (r, θ)
r ∂r ∂r r2 ∂θ2
∂2ψ
1 ∂ ∂ψ
r + Cilı́ndrico (r, z)
r ∂r ∂r ∂z 2
∂ 2 ψ senθ ∂
1 ∂ψ
+ 2 Esférico (r, θ)
∂r2 r ∂θ senθ ∂θ
A partir de la Ec. (8.8), se obtiene la relación entre la función de corriente y las componentes
del vector de velocidad, las cuales son:
∂ψ ∂ψ
vx = ; vy = − Cartesiano (x, y)
∂y ∂x
1 ∂ψ ∂ψ
vr = ; vθ = − Cilı́ndrico (r, θ)
r ∂θ ∂r
1 ∂ψ 1 ∂ψ
vz = ; vr = − Cilı́ndrico (r, z)
r ∂r r ∂z
1 ∂ψ 1 ∂ψ
vr = ; vθ = − Esférico (r, θ)
r2 senθ ∂θ rsenθ ∂r
Utilizaremos estas expresiones para estudiar el flujo potencial alrededor de un cilindro y una esfera
en lo que resta de este capı́tulo.
142
Notas de clase: Francisco J. Valdés Parada
8.2. Flujo potencial alrededor de un cilindro
Considere el flujo potencial, en estado estacionario de un fluido Newtoniano e incompresible
alrededor de un cilindro cuyo largo L excede en dos órdenes de magnitud a su radio R. De acuerdo
a lo visto anteriormente, para el flujo irrotacional, se tiene la siguiente versión de la ecuación de
Laplace:
1 ∂2ψ
1 ∂ ∂ψ
r + 2 2 =0 (8.11)
r ∂r ∂r r ∂θ
donde las relaciones entre las componentes del vector de velocidad y la función de corriente son:
1 ∂ψ ∂ψ
vr = ; vθ = − (8.12)
r ∂θ ∂r
En base a la condición de no penetración tenemos:
En r = R; n · v = vr = 0 (8.13)
O bien,
∂ψ
En r = R; =0 (8.14)
∂θ
Mientras que, suficientemente lejos de la superficies del cilindro (r R), se tiene un perfil plano y
paralelo al eje x:
r → ∞; v = v∞ ex (8.15)
dado que:
ex = cos θer − senθeθ (8.16)
Deducimos de inmediato que:
La cual es la fuente del problema de valor a la frontera y ayude a proponer una solución de la
forma:
ψ = f (r)senθ (8.20)
La cual al ser sustituida en la Ec. (8.11) da lugar a:
1 d df f
r − 2 =0 (8.21)
r ∂r ∂r r
O bien:
d2 f df
r2 +r −f =0 (8.22)
dr2 dr
143
Notas de clase: Francisco J. Valdés Parada
La cual es una ecuación de Euler-Cauchy, cuya solución puede encontrarse al utilizar la función de
prueba f = rn , lo cual de lugar a:
f = c1 r−1 + c2 r (8.24)
En r = R; f = 0; (8.25)
En r → ∞; f = v∞ r (8.26)
Con las cuales deducimos que c2 = v∞ y c1 = −v∞ R2 . Entonces, la función de corriente es:
r R
ψ = v∞ senθR − (8.27)
R r
R2
vr = v∞ cosθ 1 − 2 (8.28a)
r
R2
vθ = −v∞ senθ 1 + 2 (8.28b)
r
Por último, debemos determinar la presión en el fluido como función del radio y del ángulo. Para
ello, partimos de:
∂P ∂P
dP = dr + dθ (8.30)
∂r ∂θ
La cual debe integrarse bajo la suposición de que P = 0 en r → ∞ y θ = π:
ξ=r ζ=θ
∂P ∂P
Z Z
P= dξ + dζ (8.31)
∂ξ θ=π ∂ζ
ξ=∞ ζ=π
donde los gradientes de presión se obtienen de las componentes radial y angular de la ecuación de
Euler:
vθ2
∂vr vθ ∂vr ∂P
ρ vr + − =− (8.32a)
∂r r ∂θ r ∂r
∂vθ vθ ∂vθ vθ vr 1 ∂P
ρ vr + + =− (8.32b)
∂r r ∂θ r r ∂θ
144
Notas de clase: Francisco J. Valdés Parada
Tras sustituir los perfiles de velocidad y sus derivadas obtenemos:
2 R2 2
∂P 2ρv∞ 2 R
− = 1 − cos θ 2 (8.33a)
∂r r3 r
1 ∂P 2
4v senθ cos θ R 2
− = ∞ (8.33b)
r ∂θ r r2
Usando estos resultados en la Ec. (8.31) y llevando a cabo las integrales correspondientes resulta
finalmente:
2 R2 R2
ρv∞ 2
P= 2 − 4sen θ − 2 (8.34)
2 r2 r
145
Notas de clase: Francisco J. Valdés Parada
146
Notas de clase: Francisco J. Valdés Parada
9 | Capa lı́mite
Contenido
228 9.1. Planteamiento del problema . . . . . . . . . . . . . . . . . .Chapter
. . . . 7. 148
9.2. Solución exacta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 151
9.3. Solución aproximada . . . . . . . . . . . . . . . . . . . . . . . . . . . 155
( x) 3.6 x / u (7-123)
and in the following paragraphs we will see that this is an excellent estimate.
u∞
y
d (x)
e
v 2
v v p g v (7-125)
t
v
Assumption: 0 (7-126)
t
v
Restriction: 2 v 147 (7-127)
t
We could simplify our problem further by assuming that the flat plate is infinitely wide, or we could restrict
the width, w, according to
Notas de clase: Francisco J. Valdés Parada
9.1. Planteamiento del problema
Considere el flujo bidimensional en estado estacionario de un fluido newtoniano e incompresible
alrededor de una placa de largo L como se muestra en la Figura 9.1. Antes de chocar con la placa,
el fluido llevaba una velocidad uniforme v = u∞ ex . Al igual que en el problema del flujo cerca de
la entrada, se pueden distinguir dos regiones: una cercana al obstáculo (de espesor δ) y otra en
la zona lejos del obstáculo donde la velocidad es uniforme. Las ecuaciones gobernantes de masa y
cantidad de movimiento en la región cercana al obstáculo son:
∂vx ∂vy
+ =0 (9.1a)
∂x ∂y
∂ 2 vx ∂ 2 vx
∂vx ∂vx ∂p
ρ vx + vy =− +µ 2 +µ 2 (9.1b)
∂x ∂y ∂x ∂x ∂y
∂vy ∂vy
2
∂ vy ∂ 2 vy
∂p
ρ vx + vy =− − ρg + µ 2 + µ 2 (9.1c)
∂x ∂y ∂y ∂x ∂y
Un ejemplo de la solución numérica de estas ecuaciones se muestra en la figura 9.2 para un caso
en el que la velocidad de entrada sea de 1 m/s y pasa alrededor de una placa de 1 cm de largo.
Estas condiciones corresponden a un número de Reynolds (Re = ρu∞ L/µ) de aproximadamente
105 . Note que la capa lı́mite no es simétrica, que existe una zona de recirculación del fluido tanto
cerca de la entrada de la placa como al final de ella. También note que antes de chocar con el
obstáculo se produce otra capa lı́mite de considerable tamaño y que el espesor de la capa lı́mite
es oscilante, lo cual no concuerda con el esquema de la figura 9.1. Si se reduce en dos órdenes de
magnitud la velocidad de entrada (es decir, se fija en 0.01 m/s), el valor del número de Reynolds
baja a 103 y los resultados corresponden a los mostrados en la figura 9.3. Como puede observarse
ahora la forma de la capa lı́mite sı́ corresponde a la mostrada en la figura 9.1. Además, ya no hay
recirculación del fluido y el espesor de la capa lı́mite antes del obstáculo disminuyó; pero el de la
capa lı́mite alrededor del obstáculo aumentó. Lo anterior tiene sentido fı́sico pues cuando el flujo
es laminar, los esfuerzos viscosos son tan importantes como los inerciales y por ello la capa lı́mite
a lo largo del obstáculo aumenta.
En esta sección buscamos llevar a cabo una solución aproximada de este problema. Para
simplificar este modelo, es conveniente recurrir a los siguientes estimados de orden de magnitud:
vx =O(v∞ ) (9.2a)
x =O(L) (9.2b)
y =O(δ) (9.2c)
148
Notas de clase: Francisco J. Valdés Parada
a)
b)
Figura 9.2: a) Lı́neas de corriente de velocidad resultante de la solución numérica de las ecuaciones
de continuidad y Navier-Stokes en donde el fluido es agua. b) Filtrado de lı́neas de corriente en la
zona de capa lı́mite. La velocidad no perturbada es de 1m/s.
149
Notas de clase: Francisco J. Valdés Parada
a)
b)
Figura 9.3: a) Lı́neas de corriente de velocidad resultante de la solución numérica de las ecuaciones
de continuidad y Navier-Stokes en donde el fluido es agua. b) Filtrado de lı́neas de corriente en la
zona de capa lı́mite. La velocidad no perturbada es de 0.01 m/s.
150
Notas de clase: Francisco J. Valdés Parada
Pero como la presión no cambia con y, puede extenderse el resultado de la ecuación anterior para
toda y. Bajo estas condiciones, el modelo se reduce a:
∂vx ∂vy
+ =0 (9.11a)
∂x ∂y
∂ 2 vx
∂vx ∂vx
ρ vx + vy =µ 2 (9.11b)
∂x ∂y ∂y
La cual debe resolverse sujeta a las condiciones de frontera:
y = 0, vx = vy = 0 (9.12a)
y → ∞, vx = u∞ (9.12b)
O bien, en términos de las variables adimensionales:
vj x y
uj = (j = x, y); X= ; Y = (9.13)
u∞ L L
obtenemos:
∂ux ∂uy
+ =0 (9.14a)
∂X ∂Y
∂ux ∂ux 1 ∂ 2 ux
ux + uy = (9.14b)
∂X ∂Y ReL ∂Y 2
Y = 0; ux = uy = 0 (9.14c)
Y → ∞; ux = 1 (9.14d)
donde ReL = ρv∞ L/µ. Note que a pesar de las simplificaciones realizadas el problema es aún
no lineal debido al término inercial. A continuación presentaremos dos alternativas de solución.
La primera utiliza una descomposición de similaridad para reducir el problema a una ecuación
diferencial ordinaria no lineal de tercer orden que debe resolverse numéricamente y que fue propuesta
por Blasius. El segundo método es la solución analı́tica aproximada de von Kármán-Pohlhausen.
∂ψ ∂ 2 ψ ∂ψ ∂ 2 ψ 1 ∂3ψ
− = (9.16)
dY ∂X∂Y ∂X ∂Y 2 ReL ∂Y 3
Para resolver este problema, se propone:
151
Notas de clase: Francisco J. Valdés Parada
donde η = aX b Y c . De esta forma,
∂ψ df
= g acX b Y c−1
∂Y dη
2
∂ ψ dg df df d2 f
= acX b Y c−1 + g abcX b−1 Y c−1 + g 2 a2 bcX 2b−1 Y 2c−1
∂X∂Y dX dη dη dη
∂ψ dg df
= f + g abX b−1 Y c
∂X dX dη
2 (9.18a)
∂ ψ df d2 f
2
= g ac (c − 1) X b Y c−2 + g 2 a2 c2 X 2b Y 2(c−1)
∂Y dη dη
3
∂ ψ df df
3
= g ac (c − 1) (c − 2) X b Y c−3 + g a2 bc (c − 1) X 2b−1 Y 2c−2
∂Y dη dη
d f2 3
d f
+g 2 a2 c2 2 (c − 1) X 2b Y 2c−3 + g 3 a3 c3 X 3b Y 3(c−1)
dη dη
Como puede notarse, al escoger c = 1, las ecuaciones anteriores se reducen a:
∂ψ df
= g aX b
∂Y dη
∂2ψ dg df df d2 f
= aX b + g abX b−1 + g 2 a2 bX 2b−1 Y
∂X∂Y dX dη dη dη
∂ψ dg df
= f + g abX b−1 Y (9.18b)
∂X dX dη
2
∂ ψ 2
d f
2
= g 2 a2 c2 X 2b
∂Y dη
3
∂ ψ d3 f
3
= g 3 a3 X 3b
∂Y dη
De esta forma, la ecuación de cantidad de movimiento queda como:
2
d2 f d3 f
df 1 dg gb 1 dg
ReL + − f ReL = (9.19)
dη X b a dX X b+1 a dη 2 aX b dX dη 3
Para encontrar una solución, los términos entre paréntesis de la ecuación anterior deben ser
constantes,
ReL dg
= c1 (9.20a)
aX b dX
1 dg gb c2
b
+ b+1 = (9.20b)
X a dX X a ReL
Al resolver la primera ecuación (fijando la constante de integración como cero);
ac1 X b+1
g= (9.21)
ReL (b + 1)
Mientras que, de la segunda ecuación obtenemos:
(c2 − c1 ) X b+1 a
g= (9.22)
Re b
152
Notas de clase: Francisco J. Valdés Parada
Al igualar estos últimos resultados da lugar a:
2b + 1
c2 = c1 (9.23)
b+1
Notemos que al fijar b = −1/2 resulta que c2 = 0 y las expresiones anteriores para g se reducen a:
2ac1 X 1/2
g= (9.24)
ReL
√
Para simplificar más esta expresión se eligen: c1 = 1/2 y a = ReL , de esta forma:
r r
X ReL
g= ; η=Y (9.25)
ReL X
d3 f 1 d2 f
+ f =0 (9.26)
dη 3 2 dη 2
La cual se conoce como Ecuación de Blasius y está sujeta a las siguientes condiciones de frontera:
df
η = 0; = 0; f = 0 (9.27a)
dη
df
η → ∞; =1 (9.27b)
dη
Dada la no linealidad de este problema, no tiene solución analı́tica y debemos recurrir a técnicas
de integración numérica como el método de Euler o el método de Runge-Kutta. En ambos casos se
requiere de un esquema iterativo en donde se busca el valor en η = 0 para d2 f /dη 2 que garantice
que se satisfaga la segunda condición de frontera. Al llevar a cabo las iteraciones encontramos:
d2 f
η = 0; = 0.332056 (9.28)
dη 2
Los perfiles de velocidad resultantes de la solución numérica se muestran en la Figura 9.4. Para
obtenerlos, utilizamos:
df
η −f
df dη
ux = ; uy = √ (9.29)
dη 2 X ReL
Por último calculamos el valor del esfuerzo en la superficie de la placa:
ReL d2 f
∂vx µv∞ ∂ux µv∞ ∂ ∂ψ µv∞
τ0 = µ = = = g(X)
∂y y=0 L ∂Y Y =0 L ∂Y ∂Y Y =0 L X dη 2 η=0
r
µv∞ ReL
= 0.332 (9.30)
L X
153
Notas de clase: Francisco J. Valdés Parada
Figura 9.4: Perfiles de las componentes del vector velocidad en función de la variable de similaridad
η obtenidos numéricamente.
154
Notas de clase: Francisco J. Valdés Parada
9.3. Solución aproximada
Para llevar a cabo la solución aproximada, utilizaremos el método de von Kármán-Pohlhausen.
El primer paso de este método consiste en proponer expresiones para ux y uy . Por ello, proponemos
el siguiente polinomio de tercer orden para ux :
2 3
Y Y Y
ux = c0 + c1 + c2 + c3 (9.31)
∆ ∆ ∆
En Y = 0, ux = 0 (9.32a)
En Y = ∆, ux = 1 (9.32b)
Además, sabemos que después de la capa lı́mite no hay más transporte de momento y por tanto:
∂ux
En Y = ∆, =0 (9.32c)
∂Y
En el otro extremo, en Y = 0, el esfuerzo es máximo y por tanto:
∂ 2 ux
En Y = 0, =0 (9.32d)
∂Y 2
Utilizando estas condiciones de frontera obtenemos, c0 = 0, c1 = 3/2, c2 = 0, c3 = −1/2 y por
tanto:
1 Y 3
3 Y
ux = − (9.33)
2 ∆ 2 ∆
Para obtener uy utilizamos la ecuación de continuidad:
η=Y
∂ux
Z
uy = − dη (9.34)
∂X
η=0
donde
3 Y3
∂ux 3 Y d∆
= − 2
+ 4
(9.35)
∂X 2 ∆ 2 ∆ dX
El resultado de sustituir la expresión anterior en la integral de la Ec. (9.34) es:
η=Y " #
Y3 3 Y 2 3 Y 4 d∆
Z
3 Y 3 d∆
uy = − dη = − (9.36)
2 ∆2 2 ∆4 dX 4 ∆ 8 ∆ dX
η=0
Y2
∂ux 3 1 3
= − (9.37)
∂Y 2 ∆ 2 ∆3
155
Notas de clase: Francisco J. Valdés Parada
Estos resultados se usan en la siguiente ecuación, la cual se obtiene de integrar la Ec. (9.14b)
desde Y = 0 hasta Y = ∆ = δ/L:
YZ=∆ YZ=∆
∂ 2 ux
∂ux ∂ux 1 1 ∂ux
ux + uy dY = dY = − (9.38)
∂X ∂Y ReL ∂Y 2 ReL ∂Y Y =0
Y =0 Y =0
En X = 0, ∆=0 (9.44)
156
Notas de clase: Francisco J. Valdés Parada
donde se usó el resultado de la solución exacta. La fuerza que tenı́a el fluido antes de chocar con el
obstáculo puede calcularse como:
ρv 2
F∞ = 2W L ∞ (9.50)
2
La razón entre F y F∞ se conoce como Coeficiente de arrastre, CD , y, en este caso, es:
1.328
CD = √ (9.51)
ReL
Ejercicio 9.1.
Deduzca la expresión:
2
ρv∞
F∞ = 2W L
2
157
Notas de clase: Francisco J. Valdés Parada
Ejercicio 9.2.
Evalúe si se puede mejorar la solución aproximada desarrollada aquı́ utilizando un polinomio de
cuarto orden. Para ello, se requiere la condición de frontera adicional:
∂ 2 ux
En Y = ∆, =0 (9.52)
∂Y 2
Comente acerca del origen de esta condición de frontera y utilice el método de von Kármán para
desarrollar una solución aproximada del perfil de velocidades.
158
Chapter 9
Notas de clase: Francisco J. Valdés Parada
Turbulent Flow
In the preceding two chapters the equations of motion were solved for several types of laminar flow,
and in Sec. 7.6 we indicated how we could obtain dimensionless relations that are useful for the correlation
10 | Turbulencia
of experimental data. Experimental data are often essential for the analysis of turbulent flows in which all
three components of the velocity vector are non-zero and time-dependent. While turbulent flows are
extremely complex, the governing equations are relatively simple, i.e., the Navier-Stokes equations derived
in Chapter 6. Since many practical engineering problems are associated with turbulent flow, we must
develop some method of analyzing these flows. The first step in this analysis requires that we form the
time average of the equations of motion, and this will allow us to identify precisely what we do not know
and what we need to know. The second step requires the resolution of a closure problem and this leads to
Contenido
further analysis andPreliminares
10.1. the use of experimental
. . . . . . .data.
. . . Finally,
. . . . . it. .is. the
. . macroscopic
. . . . . . . . balances presented in
. . . . 160
Chapters 10 through 13 that provide
10.2. Suavizado en elus with the. most
tiempo . . . useful
. . . . tools
. . . for
. . the
. . analysis
. . . . . of. turbulent
. . . . . flows.
160
If we continuously measure
10.3. Problema deacerradura
component .of. the
. . .velocity
. . . . .at. some
. . . point
. . . .in. a. turbulent
. . . . . .flow,
162the result
will be similar to that illustrated in Figure 9-1. There we have shown the velocity, v x , as a rapidly
fluctuating function of time, the time average velocity, v x as a slowly changing function of time, and the
159
305
Notas de clase: Francisco J. Valdés Parada
10.1. Preliminares
El análisis hasta ahora se ha restringido a estudiar problemas en condiciones de flujo laminar.
El flujo turbulento se caracteriza por ser irregular, dependiente del tiempo y tridimensional.
Como se muestra en la Fig. 10.1, la velocidad, vx , fluctúa rápidamente con el tiempo, mientras
que la velocidad promedio en el tiempo, v x , cambia más lentamente con el tiempo. La diferencia
entre la velocidad local instantánea y su promedio en el tiempo se conoce como las desviaciones
turbulentas:
vx = v x + vx0 (10.1)
El promedio en el tiempo de una cantidad cualquiera ψ, se define como:
ζ=t+∆t
1
Z
ψ= ψdζ (10.2)
2∆t
ζ=t−∆t
El intervalo de tiempo ∆t debe ser lo suficientemente largo para suavizar en el tiempo las
fluctuaciones turbulentas, pero al mismo tiempo debe ser lo suficientemente corto que no elimine la
dependencia temporal de v x . Este tipo de descomposición es apropiada cuando existe una separación
de tiempos caracterı́sticos. Esto es, el tiempo caracterı́stico de la velocidad promedio es por lo
regular mayor que el de las desviaciones de la velocidad. Una consecuencia directa de esto es la
siguiente restricción de los campos de las desviaciones:
vx0 = 0 (10.3)
Ejercicio 10.1.
Aplique el operador de promediado a la Ec. (10.1) y explique qué debe suponerse para deducir la
Ec. (10.3)
Dado que los lı́mites de integración son independientes del tiempo, tenemos:
t+∆t
1
Z
∇· vdζ = 0 (10.6)
2∆t
ζ=t−∆t
160
Notas de clase: Francisco J. Valdés Parada
O bien,
∇·v =0 (10.7)
La cual es la forma final de la ecuación de continuidad promediada en el tiempo.
Dirigimos ahora la atención a la ecuación de Navier-Stokes. Antes de comenzar, es conveniente
expresarla, con ayuda de la ecuación continuidad, en la siguiente forma:
∂v
ρ + ∇ · (ρvv) = −∇p + ρg + µ∇2 v (10.8)
∂t
Aplicando el operador de promediado obtenemos:
t+∆t t+∆t t+∆t t+∆t
1 ∂v 1 1 1
Z Z Z Z
2
ρ dζ + ρ∇ · vvdζ = −∇ pdζ + ρg + µ∇ vdζ
2∆t ∂ζ 2∆t 2∆t 2∆t
ζ=t−∆t ζ=t−∆t ζ=t−∆t ζ=t−∆t
(10.9)
aquı́ tomamos en cuenta que la densidad es constante y además supusimos que la viscosidad no
cambia con el tiempo. Usando la definición de cantidades promedio, podemos expresar la ecuación
anterior como:
t+∆t
1 ∂v
Z
ρ dζ + ρ∇ · vv = −∇p + ρg + µ∇2 v (10.10)
2∆t ∂ζ
ζ=t−∆t
vv = (v + v0 )(v + v0 ) = v v + v0 v0 (10.14)
aquı́ se tomó en cuenta la restricción dada en la Ec. (10.3). Utilizando este resultado ası́ como la
forma integrada de la ecuación de continuidad dada en la Ec. (10.7) dan lugar a:
∂v
ρ + ρv · ∇v = −∇p + ρg + µ∇2 v − ∇ · (ρv0 v0 ) (10.15)
∂t
161
Notas de clase: Francisco J. Valdés Parada
La ecuación anterior es similar a la ecuación de Navier Stokes, excepto por el último término, el
cual se define como el tensor de esfuerzos turbulentos o esfuerzos de Reynolds:
τ = µ∇v (10.17)
τ (T ) = τ + τ (t) (10.19)
Evidentemente, para utilizar esta expresión es necesario conocer la relación entre τ (t) y v. Ante
este problema se han propuesto alternativas; por ejemplo Boussinesq razonó lo siguiente: dado que
el tensor de esfuerzos viscosos es:
τ = 2µD (10.20)
donde el tensor de deformación es D = ∇v + (∇v)T . Entonces, propuso que el tensor de esfuerzos
turbulentos sea:
τ (t) = 2µ(t) D (10.21)
donde µ(t) es la viscosidad turbulenta, la cual debe determinarse de manera experimental. Esta es
una propuesta más bien intuitiva y se ha probado que es limitada. Prandtl propuso que
(t) 2 ∂v x
µ =` (10.22)
dy
162
Notas de clase: Francisco J. Valdés Parada
ninguna de las dos cantidades, recurrimos a la deducción del problema de valor a la frontera para
las desviaciones de la velocidad.
Las versiones local-instantánea y promedio de la ecuación de continuidad son:
∇·v=0 (10.23a)
∇·v=0 (10.23b)
Notemos que:
la cual es una ecuación integrodiferencial no lineal que requiere conocer v y ∇v para poder ser
resuelta. Note que los términos que involucran cantidades promedio no están aislados, es decir, no
son términos fuente, sino términos homogéneos. En resumen, para calcular v, es necesario resolver
simultáneamente las ecs. (10.15) y (10.29), lo cual es claramente más complicado que resolver las
ecuaciones locales instantáneas. Estas dificultades han motivado el desarrollo de alternativas como
son: el desarrollo de modelos de la energı́a cinética turbulenta (K) que están acoplados con la
tasa de disipación viscosa turbulenta () o simplemente basados en la solución numérica directa
(DNS) de las ecuaciones locales instantáneas. Esta última alternativa ha motivado el desarrollo de
esquemas numéricos cada vez más robustos que permitan la solución de las ecuaciones de Navier
Stokes para valores del número de Reynolds que anteriormente eran imposibles. Si este campo
de estudio que actualmente ha cobrado especial importancia sigue progresando, es posible que en
un futuro no haya necesidad de resolver modelos turbulentos. Sin embargo, aún hace falta llevar
a cabo más investigación en este campo de la mecánica de fluidos para encontrar una respuesta
contundente.
163
Notas de clase: Francisco J. Valdés Parada
164
Notas de clase: Francisco J. Valdés Parada
11 | Balances macroscópicos
Contenido
11.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 166
11.2. Balance macroscópico de masa . . . . . . . . . . . . . . . . . . . . . 166
11.3. Balance macroscópico de cantidad de movimiento . . . . . . . . . . 170
11.3.1. Factor de fricción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 173
11.4. Ecuación de energı́a mecánica . . . . . . . . . . . . . . . . . . . . . . 175
11.5. Cálculo del coeficiente de pérdidas por disipación viscosa . . . . . 178
11.5.1. Tuberı́as lisas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 178
11.5.2. Accesorios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 178
11.6. Aplicaciones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179
11.6.1. Expansión brusca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 179
11.6.2. Flujo en tuberı́as con accesorios: Problema 1 . . . . . . . . . . . . . . 181
11.6.3. Flujo en tuberı́as con accesorios: Problema 2 . . . . . . . . . . . . . . 182
11.7. Análisis de aerogeneradores . . . . . . . . . . . . . . . . . . . . . . . 183
11.7.1. Conservación de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . 185
11.7.2. Cantidad de movimiento . . . . . . . . . . . . . . . . . . . . . . . . . . 185
11.7.3. Energı́a mecánica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 186
11.7.4. Potencia máxima de un aerogenerador . . . . . . . . . . . . . . . . . . 187
Averaging domain V
Aout
Ain
y
AL
x z
165
Notas de clase: Francisco J. Valdés Parada
11.1. Introducción
Considere un sistema cuya longitud caracterı́stica sea varios órdenes de magnitud mayor que la
longitud caracterı́stica en la cual se satisface la hipótesis del continuo. En este tipo de sistemas, a
menudo no es práctico resolver las ecuaciones vistas hasta el momento en todos los puntos. Por ello,
se recurre a técnicas de promediado para desarrollar modelos matemáticos macroscópicos a partir de
las ecuaciones gobernantes a la escala del continuo (escala microscópica). Las ecuaciones promedio,
a su vez, pueden entenderse como válidas en nuevo medio continuo donde se resolverán. Una de las
ventajas que presentan los modelos macroscópicos es que su complejidad es sustancialmente menor
a la de las ecuaciones microscópicas, además que en la práctica, la información disponible está a
menudo expresada en términos de cantidades promedio.
En este momento es pertinente recordar que el origen de las ecuaciones locales e instantáneas de
conservación (de masa, cantidad de movimiento, energı́a mecánica, etc.) pueden obtenerse a partir
de principios de conservación, los cuales ya son en sı́ macroscópicos. De hecho, en las primeras
partes del curso lo que se hizo fue plantear los principios de conservación en forma matemática y
posteriormente se usaron los teoremas del transporte y de la divergencia para obtener las ecuaciones
locales e instantáneas. Por tanto, es claro que el promediado de dichas ecuaciones lleva a las
formas integrales y macroscópicas de los principios de conservación. En otras palabras, no es
necesario promediar las ecuaciones locales e instantáneas para obtener las ecuaciones de trabajo
a escala macroscópica. Esto se mostrará a continuación, trabajando primero con el principio de
conservación de masa, seguido de la conservación de cantidad de movimiento y por último de
energı́a mecánica. Por último se presentarán problemas de aplicación donde se combinan los tres
principios de conservación.
d
Z Z
ρ dV + n · ρ(v − w) dA = 0 (11.2)
dt
Va Aa
166
Notas de clase: Francisco J. Valdés Parada
Cualquiera de estas tres versiones del principio de conservación expresan lo mismo, sólo difieren en
el tipo de región que se desee usar. Esta elección la determina la aplicación particular que se esté
estudiando.
En muchos problemas de aplicación es conveniente descomponer la superficie AA en tres: la
superficie en contacto con paredes fijas, Af , la superficie en contacto con paredes móviles, Am (t)
y las superficies de entradas y salidas Aes (t) (la cual a su vez se puede descomponer en Ae y As ).
En Af , tenemos que v = w = 0; mientras que en Am (t), v = w 6= 0. Bajo estas condiciones, la Ec.
(11.2) se expresa como:
d
Z Z Z
ρ dV = ne · ρ (v − w) dA − ns · ρ (v − w) dA (11.5)
dt
Va (t) Ae (t) As (t)
Note que el sistema coordenado se escogió de tal manera que n = −ne y n = ns . Para casos en
los que las fronteras de entrada y salida estén fijas (es decir w = 0), la ecuación anterior puede
expresarse como:
d
Z Z Z
ρ dV = ρ ne · v dA − ρ ns · v dA (11.6)
dt
Va (t) Ae (t) As (t)
donde además supusimos que la densidad es constante en las superficies de entrada y salida.
Para expresar el resultado anterior en una forma más sencilla es conveniente definir las siguientes
cantidades macroscópicas:
Z
mtot = ρ dV (11.7a)
Va (t)
1
Z
hvie = ne · v dA (11.7b)
Ae
Ae (t)
1
Z
hvis = ns · v dA (11.7c)
As
As (t)
dmtot
= ρhvie Ae − ρhvis As (11.8)
dt
Ejercicio 11.1.
167
216 Macroscopic Balances-lnertial Effects Chap. 7 S.C. 7.2 The ~acroscoplc
Notas de clase: Francisco J. Valdés Parada
•7.2 The Macrosco
.!!....
Dt
Representing the stress v
Reynolds transport theo
J ~a
j'".. <t> t
(pv) dV
Por lo tanto,
t See Prob. 4-15 for the de
d π
Z
dV = − D02 u0
dt 4
Va (t)
Como puede notarse, el sistema no tiene superficies de entradas, sólo de salidas que se localizan
en la superficie z = 0. Como esta superficie no se mueve, se tiene que w = 0, además, en esta
superficie, n = −k. Bajo estas condiciones, el balance macroscópico de masa es:
π
Z
− D02 u0 = vz dA
4
As
168
Notas de clase: Francisco J. Valdés Parada
Definiendo la velocidad promedio a la salida como:
1
Z
hvis = vz dA
TEM06 As
MODELOS TRANSITORIOS, ESTACIONARIOS Y As
CUASI-ESTACIONARIOS
donde As = π4 D12 − D02 . Tenemos entonces:
En el Tema TEM05 hemos analizado las etapas de los procesos de transporte, desde su
condición inicial hasta su comportamiento estacionario. Tanto en el flujoπ de2Couette como π
D12 − D02 hvis
en el de Hagen-Poiseuille hemos encontrado para flujo de agua y una−longitud
D0 u 0 =
característica
4 4
de 0.56 cm, que el intervalo de tiempo que podemos llamar tiempos cortos ( tc , tiempo de la
difusión de momentum
De donde en unapodemos
película de espesor igual a la
despejar lalongitud característica)
velocidad está en a la salida como:
promedio
el orden de 2 a 3 segundos, en tanto que la aproximación de la solución transitoria completa
indica que el tiempo necesario para estar cerca (a menos del 1% de diferencia) del régimen
2
estacionario, te va de 13 a 26 segundos, respectivamente. Las estimaciones anteriores D0
hvi s = −
pueden compararse con otro tiempo característico, el tiempo del transporte de masa en un 2
u0
D − D02
recipiente que se descarga td , es decir, el tiempo en el que se vacía que está en unos 90 1
segundos.
El signo negativo es de esperarse pues la velocidad a la salida está orientada en dirección contraria
Vaciado de un
altanque
eje z.
El proceso de vaciado de un recipiente a través de un tubo de descarga es transitorio, con
transporte de Ejercicio
masa y de momentum
11.2. lineal, que pueden analizarse partiendo de los balances
globales de masa y de energía mecánica. El modelo transitorio para este proceso se
construye considerando el Instrumento A, Tanque de descarga, reportado en la Parte I de
este texto, que bosquejamos aquí para su análisis:
Plano 1 Considere el sistema mostrado en la figura,
Depósito de diámetro DT donde se representa el vaciado de un tanque
H(t) z cilı́ndrico que tiene un tubo de descarga en la
Plano 2 parte inferior. Se desea determinar el tiempo de
descarga del tanque del fluido que contiene, el
cual se supone como newtoniano e incompresible.
L Solución:
Dado que el sistema no tiene entrada, sólo
Tubo de diámetro dT salida, el balance macroscópico de masa es, bajo
la suposición de flujo incompresible:
Plano 3
dVa
Z
w(t)
Figura T6.1 Diagrama del recipiente con un tubo de descarga
=− n · (v − w) dA
dt
As (t)
110 π
DT2 H(t) + d2t L
Note que, al igual que en el caso anterior, w = 0, además, n = ez y Va = 4
por lo tanto:
π 2 dH
Z
D = − vz dA = −hvz is As
4 T dt
As
DT2 dH
− = hvz is
d2T dt
169
Notas de clase: Francisco J. Valdés Parada
Donde pL = pA y p0 = pA + ρgH. De esta forma, el balance macroscópico de masa es:
ρgd2T d2T
H
ln +1 =− t+c
L 32µL DT2
32µL DT2
H0
tD = ln +1
ρgd2T d2T L
Por supuesto, este resultado está limitado por la suposición de que la ecuación de Hagen-Poiseuille
es aplicable a la salida del tubo, lo cual sabemos que no es del todo correcto.
d
Z Z Z Z
ρv dV + n · ρv(v − w) dA = ρb dV + n · (−Ip + τ ) dA (11.11)
dt
Va Aa Va Aa
170
Notas de clase: Francisco J. Valdés Parada
y para un dominio fijo
∂
Z Z Z Z
ρv dV + n · ρvv dA = ρb dV + n · (−Ip + τ ) dA (11.13)
∂t
Vf Af Vf Af
Para simplificar esta expresión, adoptaremos la suposición que en las superficies de entrada y salida
los esfuerzos superficiales normales son mucho mayores que los tangenciales, esto es:
en Aes , τ Ip (11.16)
Dirigiendo por el momento la atención al primer término del lado derecho de la ecuación anterior,
podemos descomponer a la velocidad en términos de su valor promedio y desviaciones espaciales:
v = hvie + ṽ (11.18)
donde
1
Z
hvie = v dA (11.19)
Ae
Ae
hṽie = 0 (11.20)
171
Notas de clase: Francisco J. Valdés Parada
Para situaciones en las que el primer término sea mucho mayor que el segundo (por ejemplo cuando
la velocidad es constante en la entrada), tenemos que
Z Z
ne · vv dA = ne dA · hvie hvie (11.22)
Ae (t) Ae (t)
Más aún, para casos en los que el vector unitario ne sea constante en la sección de la entrada,
tenemos: Z
ne · vv dA = Ae ne · hvie hvie (11.23)
Ae (t)
hvie
λe = (11.24)
hvie
Por último, si el flujo está orientado en la misma dirección que la entrada, es decir ne = λe ,
tendremos: Z
ne · vv dA = Ae λe hvi2e (11.26)
Ae (t)
Si aplicamos estas mismas suposiciones a la sección de salida y suponemos además que las presiones
a la entrada y salida del sistema son constates, podemos expresar la Ec. (11.17) como
donde se definió al vector de fuerza que ejerce el fluido sobre el sólido como:
Z Z
F= n · (Ip − τ ) dA = − n · T dA (11.28)
Af (t)+Am (t) Af (t)+Am (t)
Por lo tanto
172
Notas de clase: Francisco J. Valdés Parada
11.3.1. Factor de fricción
A menudo en las superficies fijas y móviles es razonable suponer que
pI τ (11.31)
A∗ = 2πRL (11.35a)
1
K ∗ = ρhvi2 (11.35b)
2
donde en el caso particular de flujo laminar en una tuberı́a, la ecuación de Hagen-Poiseuille lleva a
(P0 − PL )R2
hvz i = (11.36)
8µL
Del balance macroscópico de cantidad de movimiento tenemos que:
Por lo que el factor de fricción para condiciones de flujo laminar desarrollado es:
1 D P0 − PL 16
f= 1 2 = (11.39)
4 L 2 ρ hvz i Re
Note que la primera igualdad no está limitada a flujo laminar pero la segunda si lo está. En la
ecuación anterior, el número de Reynolds se define como:
ρhvz i2R
Re = (11.40)
µ
173
Notas de clase: Francisco J. Valdés Parada
La expresión anterior para el factor de fricción funciona bien para Re < 2100 en tubos lisos. Para
números de Reynolds mayores, se utiliza la fórmula de Blasius:
0.0791
f= , 2100 < Re < 105 (11.41)
Re1/4
Para tubos rugosos se utiliza la siguiente expresión empı́rica de Haaland
" #
k/D 10/9
1 6.9
√ = −3.6 log10 + (11.42)
f Re 3.7
Donde k/D es el factor de rugosidad (ver Figura 6.6-2 del texto de Bird y col. (2007)). La ecuación
anterior puede usarse para 4 × 104 < Re < 108 , 0 < k/D < 0.05.
Ejercicio 11.3.
El radio hidráulico se define como
Área de la sección transversal
Rh = (11.43)
Perı́metro de la sección transversal
Demuestre que para el flujo en una tuberı́a se tiene que D = 4Rh y que
Rh P0 − PL
f=
L 12 ρ hvz i2
Mientras que para valores aún mas grandes del número de Reynolds el resultado tiende a una
constante:
f = 0.44, 500 < Re < 105 (11.46)
Ejercicio 11.4.
Considere el flujo en estado estacionario de un fluido Newtoniano e incompresible a través de una
tuberı́a de D = 10 cm. El fluido es agua y tiene una densidad de 1000 kg/m3 y una viscosidad de
1 × 10−3 Pa s. Determine el gradiente de presión necesario para conseguir un flujo másico de 1 kg/s.
Solución: Primeramente se determina la velocidad del fluido como:
ṁ
hvz i = = 0.127m/s
ρA
Con los datos que aporta el problema podemos calcular el número de Reynolds como:
ρhvz iD
Re = = 1.27 × 104
µ
174
Notas de clase: Francisco J. Valdés Parada
Usando la fórmula de Blasius obtenemos que el factor de fricción es f = 7.45 × 10−3 . Por último,
usando la primera igualdad de la ecuación 11.39 obtenemos que
P0 − PL Pa
= 2.41
L m
d ρ 2 ρ
Z Z Z Z
v dV + n · v 2 (v − w) dA = ρb · v dV + n · (T · v) dA
dt 2 2
Va (t) Aa (t) Va (t) Aa (t)
Z
+ (p∇ · v − τ : ∇v) dV (11.48)
Va (t)
Aunque esta forma es útil para estudiar una gran variedad de problemas, es conveniente expresarla
en términos de la combinación de la energı́a cinética y potencial. Para ello, es conveniente restringir
el análisis a situaciones en las que la fuerza volumétrica se deba a la gravedad (esto es, b = g).
Además, el vector de gravedad puede expresarse en función de la energı́a potencial por unidad de
masa (φ̂ = −gh), por lo que
g = −∇φ̂ (11.49)
De esta forma, el término del trabajo debido a fuentes volumétricas se expresa como
Z Z
ρb · v dV = − ρ∇φ̂ · v dV (11.50)
Va (t) Va (t)
175
Notas de clase: Francisco J. Valdés Parada
Sustituyendo este resultado en la Ec. (11.50) se obtiene
∂ρφ̂
Z Z Z
ρb · v dV = − ∇ · ρφ̂v dV − dV (11.53)
∂t
Va (t) Va (t) Va (t)
d
Z Z Z h i
ρb · v dV = − ρφ̂ dV − n · ρφ̂ (v − w) dA (11.54)
dt
Va (t) Va (t) Aa (t)
Usando este resultado en la Ec. (11.48) se obtiene la siguiente forma de la ecuación macroscópica
de energı́a mecánica
Z
d 1 2 1 2
Z
ρv + ρφ̂ dV + n· ρv + ρφ̂ (v − w) dA =
dt 2 2
Va (t) Aa (t)
Z Z
n · (T · v)dA + (p∇ · v − τ : ∇v) dV (11.55)
Aa (t) Va (t)
Para simplificar la ecuación macroscópica de energı́a mecánica definimos las energı́as cinética y
potencial totales como:
Z
1 2
Ktot = ρv dV (11.56a)
2
Va (t)
Z
φtot = ρφ̂ dV (11.56b)
Va (t)
Al igual que en los casos anteriores podemos expresar el término relacionado con el flujo neto
de energı́a cinética a través de las fronteras del dominio como:
1 2 1 2
Z Z
n· ρv + ρφ̂ (v − w) dA = n· ρv + ρφ̂ vdA (11.57)
2 2
Aa (t) Aes (t)
Más aún, dado que la velocidad del fluido es nula en las superficies fijas, tenemos que
Z Z Z
n · (T · v)dA = n · (T · v)dA + n · (T · v)dA (11.58)
Aa (t) Aes (t) Am (t)
El último término del resultado anterior suele identificarse como la tasa a la cual las superficies
móviles ejercen trabajo sobre el fluido, esto es:
Z
Wm = n · (T · v)dA (11.59)
Am (t)
176
Notas de clase: Francisco J. Valdés Parada
Por último, las pérdidas de energı́a cinética se dan por dos mecanismos: por mecanismos reversibles
relacionados con la compresibilidad del fluido y por mecanismos irreversibles debidos a la disipación
viscosa. Esto es:
Z
Ec = − p(∇ · v)dV (11.60a)
Va (t)
Z
Eν = τ T : ∇vdV (11.60b)
Va (t)
Bajo estas condiciones, el balance macroscópico de energı́a mecánica puede expresarse como:
d 1 2
Z Z
(Ktot + φtot ) + n· ρv + ρφ̂ vdA = n · (T · v)dA + Wm − Ec − Eν (11.61)
dt 2
Aes (t) Aes (t)
Este resultado puede simplificarse aún más si se imponen las siguientes suposiciones:
i) Estado estacionario.
ii) En las superficies de entradas y salidas los esfuerzos viscosos son despreciables respecto a los
esfuerzos de presión y que la presión es constante en cada superficie.
iv) Es razonable despreciar las desviaciones espaciales de la velocidad en las superficies de entrada
y salida.
Por lo tanto,
1 2 ps 1 2 pe
hvis + ghs + − hvie + ghe + = Ŵm − Êν (11.64)
2 ρ 2 ρ
donde Ŵm = Wm /ṁ y Êν = Eν /ṁ. Para casos en los que no se ejerce trabajo por fuentes móviles
y las pérdidas por disipación viscosa sean despreciables, la ecuación anterior se reduce a la bien
conocida ecuación de Bernoulli:
1 1
ρhvi2e + ρghe + pe = ρhvi2s + ρghs + ps (11.65)
2 2
177
Notas de clase: Francisco J. Valdés Parada
11.5. Cálculo del coeficiente de pérdidas por disipación viscosa
11.5.1. Tuberı́as lisas
En el caso de tuberı́as de sección transversal constante, el balance macroscópico de cantidad de
movimiento mostramos que se reduce a:
11.5.2. Accesorios
Como puede notarse en la Ec. (11.68) hay un factor que corrige a la energı́a cinética debido a
la razón de aspecto L/R y el factor de fricción en tuberı́as. Este resultado ha motivado proponer
que, para porciones con accesorios como válvulas, codos, etc., se estimen las pérdidas por disipación
viscosa como sigue:
1
Êν = hvz i2 eν (11.70)
2
donde eν es el factor de pérdidas por fricción y adquiere valores especı́ficos para cada accesorio
como se mostrará más adelante. Por el momento, es conveniente escribir la ecuación de energı́a
mecánica tomando en cuenta estos resultados:
X 1 X
1 2 ps 1 2 pe 2 L 1 2
hvi + ghs + − hvi + ghe + = Ŵm − hvz i f − hvz i eν
2 s ρ 2 e ρ 2 Rh i 2 i
i i
| {z } | {z }
Contribución de tuberı́as Accesorios
(11.71)
Aquı́ usamos la definición del radio hidráulico en tuberı́as: Rh = R/2. Algunos autores, prefieren
expresar
2Le
eν = f (11.72)
R
donde Le es la longitud equivalente. Usando esta definición, la ecuación de energı́a mecánica da
como resultado:
X 1
1 2 ps 1 2 pe 2 L + Le
hvi + ghs + − hvi + ghe + = Ŵm − hvz i f (11.73)
2 s ρ 2 e ρ 2 Rh i
i
178
Notas de clase: Francisco J. Valdés Parada
11.6. Aplicaciones
Considere la situación mostrada en la Fig. 11.1 que consiste en el flujo en estado estacionario
de un fluido Newtoniano a través de una tuberı́a que súbitamente se ensancha. Como se puede
apreciar por las simulaciones numéricas, la velocidad incrementa los remolinos y zonas muertas
cuando se incrementa el número de Reynolds. Se busca obtener una expresión para las pérdidas por
disipación viscosa en este caso. Para ello, se analizarán los tres balances macroscópicos revisados
hasta el momento:
Balance de masa:
Bajo condiciones de flujo incompresible, la versión macroscópica de la ecuación de continuidad
es:
Z
Mechanical energy balance n · v dA = 0391 (11.74)
Aeys
v n dA
A
0 (11-63)
and for the control volume illustrated in Figure 11-4 this reduces to
v3 A3 v1 A1 (11-64)
Here we have used v1 to represent n v at the entrance of the control volume and v3 to represent n v at
the exit.
2
3
a) b)
pressure
z0 zL
Figure 11-4. Sudden expansion in a pipeline
t
(T)
v12 A1 v32 A3 k dA (11-66)
(n ) 179
A
in which we have assumed that k g 0 . The momentum flux is generally quite easy to evaluate since we
need only direct our attention to the entrances and exits; however, the area integral of the stress vector is
more complex and in this problem we need to consider the five distinct areas illustrated in Figure 11-4.
These are: (1) the entrance, (2) the upstream pipe wall, (3) the annular region at the sudden expansion, (4)
the downstream pipe wall, and (5) the exit.
Notas de clase: Francisco J. Valdés Parada
La cual da lugar de inmediato a:
hvz i1 A1 = hvz i3 A3 (11.75)
O bien:
hvz i1 D2
= 32 (11.76)
hvz i3 D1
Balance de cantidad de movimiento
La forma macroscópica de la ecuación de conservación de cantidad de movimiento es
F = Ae λe ρhvi2e + pe − As λs ρhvi2s + ps + mtot g
(11.77)
O bien, obteniendo su componente z y retomando la definición del vector de fuerza,
Z
2 2
ρhvz i3 A3 − ρhvz i1 A1 = hpi1 A1 − hpi3 A3 − n · (Ip − τ ) · ez dA (11.78)
Af (t)
Por último, dado que nos encontramos a la entrada de una nueva sección de flujo, es razonable pensar
que los esfuerzos viscosos no sean máximos aún, por ello, adoptaremos la siguiente suposición:
Z
τrz dA hpi1 (A3 − A1 ) (11.81)
Aparedes
180
08 Chapter 7 Macroscopic Balances for Isothermal Flow Systems
Notas de clase: Francisco J. Valdés Parada
Fig. 7.5-1. Pipeline flow
with friction losses be-
cause of fittings. Planes
1and 2 are just under
- - Plane 2 the surface of the liquid
Figura 11.2: Sistema de tuberı́as con accesorios. Las longitudes en metros son: 1.5 m, 91 m, 30.5
m, 36.6 m, 6.1 m.
A3 A23 A3 2
1 2 1 2
Êν = hvz i3 1 − 2 + = hvz i3 1 − (11.85)
2 A1 A21 2 A1
and
Este resultado seráthe
de Reynolds number
gran utilidad is
en el siguiente problema para estimar las pérdidas por expansión
brusca.
181
This is the work (per unit mass of fluid) done on the fluid in the pump. Hence the pump doe
2830 ft2/s2or 2830/32.2 = 88 ft lbf/lbrnof work on the fluid passing through the system. Th
mass rate of flow is
Notas de clase: Francisco J. Valdés Parada
Sin embargo, dado que no hay cambios en el diámetro de la tuberı́a, no hay cambios en la velocidad.
Además a la entrada y a la salida los recipientes se encuentran expuestos a la presión atmosférica.
Por lo tanto, la ecuación anterior se reduce a:
X 1 X 1
2 2L 2
Ŵm = g(hs − he ) + hvz i f + hvz i eν
2 R i 2 i
i i
Aquı́ usamos el hecho de que Rh = R/2. Además, del esquema de la figura tenemos que hs − he =
25.9m. Además,
X 1 m2
2 2L 1.5 + 91 + 30.5 + 36.6 + 6.1
hvz i f = (0.72)2 (0.0048) = 8.28 2
2 R i 0.05 s
i
m2
Ẇm = 263.13
s2
O bien,
Wm = ṁẆm = ρQ̇Ẇm = 1.49kW
X X 1
1 2 1 2 2L 2
hvi + ghs − hvi + ghe =− hvz i f − hvz i eν
2 s 2 e R i 2 i
i i
Q̇ m
hvie = = 0.233
Ae s
182
Notas de clase: Francisco J. Valdés Parada
Figura 11.3: Sistema de tuberı́as con accesorios. Los diámetros de las tuberı́as son 77.9 mm y 40.9
mm, respectivamente. Las válvulas se consideran de compuerta en el tanque y de globo a la salida.
y de la misma forma, la velocidad de salida es 0.846 m/s. Con esta información calculamos:
Ree = 1.82 × 104 ; Res = 3.46 × 104
Como el régimen es turbulento, usamos la fórmula de Blasius para obtener:
0.0791
fe =1/4
= 0.0068
Ree
y fs = 0.0058. Con esta información podemos calcular:
X m2
2L
hvz i f = 17.14 2
R i s
i
Además, tomando en cuenta las pérdidas por accesorios tenemos:
X 1 m2
2
hvz i eν = 5.92 2
2 i s
i
Con estos resultados, obtenemos que la altura del tanque debe ser de 0.38 m.
183
Notas de clase: Francisco J. Valdés Parada
Averaging domain V
Aout
Ain
y
AL
x z
a)
AM,F AM,B
AM,H
AF
x
b) z
Figura 11.4: Esquema del sistema y sus fronteras concibiendo al aerogenerador como un disco
actuador.
184
Notas de clase: Francisco J. Valdés Parada
11.7.1. Conservación de masa
Suponiendo condiciones de estado estacionario y considerando al aire como un fluido incompresible,
el balance macroscópico de masa se reduce a
Z
n · (v − w) dA = 0 (11.86)
Aa
donde Aa puede descomponerese en las superficies de entrada y salida (Ain y Aout ), las superficies
laterales (AL ), las superficies fijas (AF ), la nariz del aerogenerador (Am,H ) ası́ como las superficies
frontal y trasera del disco actuador (Am,F y AM,B ). Dado que en las superficies móviles v = w y
en las superficies fijas v = w = 0 y en el resto de las superficies w = 0, la ecuación anterior se
reduce a Z Z Z
nin · vin dA + nout · vout dA + nL · vL dA = 0 (11.87)
Ain Ain Ain
o bien,
− Ain hvin iin + Aout hvout iout + AL hvL iL = 0 (11.88)
A la entrada la velocidad del aire es máxima, por lo que suponiendo que
185
Notas de clase: Francisco J. Valdés Parada
y que las desviaciones espaciales de la velocidad son despreciables en las superficies de entrada y
salida, se tiene que Z
− ṁ(v1 − v2 ) = n · T · ez dA (11.95)
Aa
Dirigiendo la atención al lado derecho de esta ecuación, descomponiendo la integral en todas sus
contribuciones y suponiendo que en las superficies de entrada, salida y laterales, los esfuerzos
viscosos son despreciables respecto a los esfuerzos cortantes, el resultado es
Z Z Z Z Z
− ṁ(v1 − v2 ) = pin dA − pout dA − pL dA + n · T · ez dA + n · T · ez dA
Ain Aout AL AM,F AM,B
Z Z
+ n · T · ez dA + n · T · ez dA (11.96)
AM,H AF
De estos términos, aquellos asociados con las superficies frontal y trasera del disco actuador son los
mas grandes en magnitud, por lo que se puede suponer que
Z Z Z Z Z
pin dA − pout dA − pL dA + n · T · ez dA + n · T · ez dA
Ain Aout AL AM,H AF
Z Z
n · T · ez dA + n · T · ez dA (11.97)
AM,F AM,B
donde se tomó en cuenta que AM,F = AM,B = AM . Esta expresión se utiliza para calcular la fuerza
que imprime el viento sobre el aerogenerador.
Usando argumentos similares a los utilizados anteriormente la ecuación anterior se puede reformular
como
ṁ 2
Z
2
− (v1 − v2 ) = AM (hTn,z vz iF − hTn,z vz iB ) − τ : ∇v dV (11.100)
2
VA
Aquı́ se supuso que la única componente relevante de la velocidad alrededor del aerogenerador es la
componente z. En muchas aplicaciones es razonable suponer que el esfuerzo máximo se experimenta
186
Notas de clase: Francisco J. Valdés Parada
en la superficies AM,F por lo que se supone ahora que la disipación viscosa es despreciable respecto
a este término, esto es Z
τ : ∇v dV AM hTn,z vz iF (11.101)
VA
Más aún, despreciando los cambios espaciales de la velocidad del viento en las superficies de entrada
y salida del disco actuador, la ecuación (11.100) puede escribirse como
ṁ 2
P = vω AM (hTn,z iF − hTn,z iB ) = vω F = − (v − v22 ) (11.102)
2 1
la cual es una expresión de la potencia del viento sobre el aerogenerador. Este último resultado
puede trabajarse más sustituyendo la ecuación (11.98), para obtener
1
vω = (v1 + v2 ) (11.103)
2
Esta ecuación se conoce como teorema de Froude y es común en turbomaquinaria.
vω = v1 (1 − a) (11.105)
Usando este resultado en la ecuación (11.103) y despejando la velocidad de salida resulta que
v2 = v1 (1 − 2a) (11.106)
F = −ṁv1 2a (11.107)
Fa = ṁv1 2a (11.108)
Más aún, dado que no se pierde masa al pasar por el aerogenerador, se tiene que
187
Notas de clase: Francisco J. Valdés Parada
La potencia del viento antes de entrar al sistema es
1 2 ρ
Pin = v1 ṁ = AM v13 (11.112)
2 2
d2 CP
= 8(−2 + 3a) (11.116)
da2
al sustituir el valor crı́tico se obtiene
d2 CP
= −8 < 0, ∴ a = 1/3 corresponde a un máximo (11.117)
da2 a=1/3
188
Notas de clase: Francisco J. Valdés Parada
Referencias
Aris, R. (1976). How to get the most output of an equation without really trying. Chemical
Engineering Education 10, 114–124.
Basmadjian, D., Farnood, R. (2007). The Art of Modeling in Science and Engineering with
Mathematica. Chapman & Hall. segunda edición.
Bear, J., Cheng, A. (2010). Modeling Groundwater Flow and Contaminant Transport. Springer.
Bingham, E. (1916). An investigation of the laws of plastic flow. U.S. Bureau of Standards Bulletin
13, 309–353.
Bird, R., Stewart, W., Lightfoot, E. (2007). Transport Phenomena. Wiley. segunda edición.
Gibbs, J. (1957). The collected works of J. Willard Gibbs. Yale University Press.
Herrera, I., Pinder, G. (2012). Mathematical Modeling in Science and Engineering. Wiley.
Ochoa-Tapia, J., del Rı́o, J., Whitaker, S. (1993). Bulk and surface diffusion in porous media: an
application of the surface averaging theorem. Chemical Engineering Science 48, 2061–2082.
Ostwald, W. (1929). Ueber die rechnerische darstellung des strukturgebietes der viskosität. Kolloid
Zeitschrift 47, 176–187.
Pinder, G., Gray, W. (2008). Essentials of Multiphase Flow in Porous Media. Wiley.
Slattery, J., Oh, E., Fu, K. (2004). Extension of continuum mechanics to the nanoscale. Chemical
Engineering Science 59, 4621–4635.
Slattery, J., Sagis, L., Oh, E. (2006). Interfacial Transport Phenomena. Springer. segunda edición.
189
Notas de clase: Francisco J. Valdés Parada
Valdés-Parada, F. (2014). Breviario sobre modelado matemático. Universidad Autónoma
Metropolitana.
Whitaker, S. (1992b). The species mass jump condition at a singular surface. Chemical Engineering
Science 47, 1677–1685.
Wood, B. (2009). The role of scaling laws in upscaling. Advances in Water Resources 32, 723–736.
190