Actividad 1 –Fluidos Computacionales

Nombres: Juan José Arenas Londoño – María Alejandra Pulgarín Muñoz.

ID: 000428010 - 000148177
Ecuaciones gobernantes y Solución analítica de ecuaciones de transporte (30%)

Objetivos:

• Identificar la importancia de las ecuaciones gobernantes en forma diferencial y su

importancia en el estudio de flujo de fluidos.
• Solucionar de forma analítica la ecuación de cantidad de movimiento y conservación de
masa en forma diferencial para determinar perfiles de velocidad en una geometría dada.

Instrucciones

1. Resuelva las preguntas y el problema que se presentan a continuación.

2. Envié el documento, en un archivo con extensión. pdf, a través del acceso creado en la
plataforma.

Preguntas:

1. Simplificar la ecuación general de continuidad:

para el caso de flujo bidimensional con densidad constante.

R=/ 0 0

𝜕𝑢 𝜕𝑣
𝜌 (𝜕𝑋 + 𝜕𝑦) = 0

𝜌=0

𝜕𝑢 𝜕𝑣
+ =0
𝜕𝑥 𝜕𝑦

1
2. ¿Cuál es la segunda ley del movimiento de Newton? (para fluidos en movimiento) Escribir
ecuación y significado de cada término.

R=/ La fórmula de la segunda ley del movimiento de Newton:

⃗
𝜕𝑉
𝜌( ⃗ ∗ ∇𝑉
+𝑉 ⃗ ) = −∇𝜌 + 𝜇∇2 𝑉
⃗ +𝐹
𝜕𝑡
Significado de cada termino:
⃗
𝜕𝑉
𝜌 ∗ ( ) : Aceleración temporal: es la tasa de cambio de la velocidad de una partícula del
𝜕𝑡
fluido en un punto específico con respecto al tiempo.

⃗ ∗ ∇𝑉
𝜌 ∗ (𝑉 ⃗ ) : Aceleración advección: es un componente de la aceleración temporal y es
un cambio en la velocidad de una partícula del fluido debido al movimiento del fluido.

−∇𝜌: Gradiente de presión: la variación de la presión del fluido por unidad de distancia.

𝜇∇2 𝑉⃗ : Viscosidad: es la dificultad que representa el fluido para deformarse y fluir cuando
se le aplica una fuerza externa.

𝐹 : Fuerza externa: este incluye cualquiera fuerza aplicada externamente al fluido, este
incluye la gravedad.

3. Un flujo laminar incompresible, no viscoso, unidimensional simplificado se define

mediantela siguiente ecuación de cantidad de movimiento en la dirección x:

Nombre cada término y discuta su contribución al flujo.

R=/ Esta fórmula describe cómo la velocidad del flujo cambia debido a la combinación de
variaciones temporales, movimiento a través del espacio.
𝜕𝑢
: Aceleración temporal
𝜕𝑡

𝜕𝑢
𝑢 : Advección
𝜕𝑥

Juntos, estos términos nos brindan información crucial sobre el comportamiento de la

velocidad del fluido: su valor en un instante determinado y cómo cambia a medida que la
partícula fluye.

Y por lo último en el lado derecho:

1 𝜕𝑝
− : Gradiente de presión
𝜌 𝜕𝑥
2
 4. ¿Qué propiedad del fluido cambia de constante a variable cuando el fluido se considera
un flujo compresible? ¿Cuál ecuación se usa complementariamente para evaluar esa
propiedadpara esa condición del flujo?

R=/ La densidad cambia de constante a variable en flujo compresible por los efectos de la
presión y la densidad.

𝜕𝜌 𝜕(𝜌𝑢) 𝜕(𝜌𝑣) 𝜕(𝜌𝑤)

+ + + =0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

5. Escriba la ecuación de continuidad para un flujo compresible e incompresible. ¿Qué

diferencias notas en las ecuaciones?

R=/ En la primera ecuación la densidad es constante:

𝜕𝑢 𝜕𝑣 𝜕𝑤
+ + =0
𝜕𝑥 𝜕𝑦 𝜕𝑧

Y en la segunda ecuación la densidad es variable:

𝜕𝜌 𝜕(𝜌𝑢) 𝜕(𝜌𝑣) 𝜕(𝜌𝑤)
+ + + =0
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

6. En qué tipo de situación se puede simplificar la ecuación de energía 2-D general para
alcanzar la conocida ecuación de Laplace:

R=/ Para alcanzar la ecuación de Laplace las dos condiciones necesarias son: Flujo estacionario y flujo
en condición estática (sin movimiento en X o Y)
0 0 0

𝑘
=0
𝜌𝐶𝑝

Problema - Solución analítica de ecuaciones de transporte

Flujo entre dos cilindros en rotación

Existe una solución exacta para las ecuaciones de Navier-Stokes para el caso del flujo de un fluido

3
contenido entre dos cilindros concéntricos, incluso cuando ambos cilindros rotan a velocidades
constantes. Los cilindros son suficientemente largos comparados con el diámetro para asumir que
el flujo es bidimensional.

En la siguiente figura se observa la configuración, donde 𝑟𝑖 se refiere al radio del cilindro interno que
gira a una velocidad 𝜔𝑖 y 𝑟𝑒 representa el radio del cilindro externo que rota a velocidad 𝜔𝑒.

Figura 1. Geometría para el flujo entre dos cilindros concéntricos en rotación.

Para los valores indicados en la tabla a continuación, y empleado las ecuaciones de Navier-Stokes
expresadas en coordenada cartesianas, dar respuesta a las siguientes preguntas.

Tabla 1. Condiciones estudio

Variable Valor
𝜔𝑖 1 rad/s
𝜔𝑒 0 rad/s
𝑟𝑖 1 cm

1. Demuestre que las ecuaciones de Navier-Stokes se reducen a la siguiente forma:

R=/

𝜕𝑉𝑟 𝜕𝑉𝑟 𝑉𝜃 𝜕𝑉𝑟 𝜕𝑉𝑟 𝑉𝜃2

𝑟: 𝜌 ( + 𝑉 + + 𝑉𝑍 − )
𝜕𝑡 𝜕𝑟 𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝑟
𝜕𝑝 𝜕 1𝜕 1 𝜕 2 𝑉𝑟 𝜕 2 𝑉𝑟 2 𝜕𝑉𝜃
=− + 𝜌𝑔𝑟 + 𝜇[ ( (𝑟𝑉𝑟 )) + 2 + − 2 ]
𝜕𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 𝑟 𝜕𝜃

𝑉𝜃2 𝜕𝑝
𝑟: 𝜌 (− )=−
𝑟 𝜕𝑟

𝑉𝜃2 1 𝑑𝑝
𝑟: − =−
𝑟 𝜌 𝑑𝑟
4
 𝜕𝑉𝜃 𝜕𝑉𝜃 𝑉𝜃 𝜕𝑉𝜃 𝜕𝑉𝜃 𝑉𝑟 𝑉𝜃
𝜃: 𝜌 ( + 𝑉𝑟 + + 𝑉𝑍 − )
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝜕𝑧 𝑟
1 𝜕𝑝 𝜕 1𝜕 1 𝜕 2 𝑉𝜃 𝜕 2 𝑉𝑟 2 𝜕𝑉𝑟
=− + 𝜌𝑔𝜃 + 𝜇 [ ( (𝑟𝑉𝜃 )) + 2 + − ]
𝑟 𝜕𝜃 𝜕𝑟 𝑟 𝜕𝑟 𝑟 𝜕𝜃 2 𝜕𝑧 2 𝑟 2 𝜕𝜃

𝜃: 𝜇 = 0

𝑑 1𝑑
=[ (𝑟𝑣𝜃 )] = 0
𝑑𝑟 𝑟 𝑑𝑟

2. Para condiciones de no deslizamiento sobre las paredes de los cilindros y asumiendo que
únicamente gira el cilindro interior, demuestre que la distribución de velocidades del fluido
entre ambos cilindros es:

R=/
𝑑 1𝑑
[ ∗ (𝑟𝑢𝜃 )] = 0 𝑠𝑒 𝑖𝑛𝑡𝑒𝑔𝑟𝑎 𝑎 𝑎𝑚𝑏𝑜𝑠 𝑙𝑎𝑑𝑜𝑠
𝑑𝑟 𝑟 𝑑𝑟

1𝑑
(𝑟𝑢𝜃 ) = 𝐶1
𝑟 𝑑𝑟

𝑑
∫ (𝑟𝑢𝜃 ) = 𝐶1 ∫ 𝑟
𝑑𝑟

𝐶1 𝑟 2
𝑟𝑢𝜃 = + 𝐶2
2

𝐶1 𝑟 𝐶2
𝑢𝜃 = +
2 𝑟

Condiciones de frontera

1. 𝑟 = 𝑟𝑖 𝑢𝜃 = 𝑤𝑖 𝑟𝑖
2. 𝑟 = 𝑟𝑒 𝑢𝜃 = 0
𝐶1 𝑟 𝐶2
𝑢𝜃 = +
2 𝑟𝑖

1era condición.

5
 𝐶1 𝑟𝑖 𝐶2
𝑤𝑖 𝑟𝑖 = +
2 𝑟𝑖

𝐶1 𝑟𝑖2
𝑤𝑖 𝑟𝑖2 = + 𝐶2
2
2da condición.

𝐶1 𝑟𝑒 𝐶𝑒
0= +
2 𝑟𝑒

𝐶1 𝑟𝑒2
0= + 𝐶2
2

𝐶1 𝑟𝑒2
𝐶2 = −
2

Siguiente.
𝐶1 𝑟𝑖2 𝐶1 𝑟𝑒2
𝑤𝑖 𝑟𝑖2 = −
2 2

2𝑤𝑖 𝑟𝑖2
= 𝐶1
𝑟𝑖2 − 𝑟𝑒2

2𝑤𝑖 𝑟𝑖2 𝑟 𝑤𝑖 𝑟𝑖2 𝑟𝑒2

𝑢𝜃 = −
2(𝑟𝑖2 − 𝑟𝑒2 ) (𝑟𝑖2 − 𝑟𝑒2 )𝑟

𝑤𝑖 𝑟𝑖2 𝑟 𝑟𝑒2
𝑢𝜃 = (𝑟 − )
(𝑟𝑖2 − 𝑟𝑒2 ) 𝑟

𝑤𝑖 𝑟 1
𝑢𝜃 = 𝑟𝑖2 ∗ 𝑟𝑒2 ( )( − )
1 1 𝑟𝑒2 𝑟
−
𝑟𝑒2 𝑟𝑒2

𝑤𝑖 𝑟𝑖2 ∗ 𝑟𝑒2
𝑢𝜃 =
1 1
𝑟𝑖2 ∗ 𝑟𝑒2 ( 2 ) − 𝑟𝑖2 ∗ 𝑟𝑒2 ( 2 )
𝑟𝑒 𝑟𝑖

𝑤𝑖 𝑟𝑖2 ∗ 𝑟𝑒2 𝑟 1
𝑢𝜃 = 2 2 ∗ ( 2 − )
𝑟𝑖 ∗ 𝑟𝑒 𝑟𝑒 𝑟

6
 𝑤𝑖 𝑟𝑖2 ∗ 𝑟𝑒2 1 𝑟 1
𝑢𝜃 = 2 2 ∗ 2 2 ( 2− )
𝑟𝑖 ∗ 𝑟𝑒 𝑟𝑖 ∗ 𝑟𝑒 𝑟𝑒 𝑟
1
2
𝑟𝑖 ∗ 𝑟𝑒2

𝑤𝑖 𝑟 1
𝑢𝜃 = ( − )
1 𝑟𝑒2 𝑟𝑒2 𝑟
𝑟𝑖2 ∗ −
𝑟𝑖2 ∗ 𝑟𝑒2 𝑟𝑖2 ∗ 𝑟𝑒2

𝑤𝑖 𝑟 1
𝑢𝜃 = ( − )
1 1 𝑟𝑒2 𝑟
−
𝑟𝑒2 𝑟𝑖2

3. Considere las siguientes condiciones de frontera con deslizamiento sobre las fronteras,
dadopor la longitud de deslizamiento 𝐿𝑠.

Obtenga una solución similar a la presentada en el numeral anterior, en la que se incluya el

efecto de 𝐿𝑠 sobre el perfil de velocidades 𝑢𝜃. Demuestre que para 𝐿𝑠 = 0, la solución
corresponde al caso desarrollado previamente para no deslizamiento.

R=/
𝐶2
𝑢(𝑟) = 𝐶1 𝑟 + ( )
𝑟

𝜕𝑢𝜃 𝐶2
= 𝐶1 − ( 2 )
𝜕𝑟 𝑟

𝑢𝜃 𝐶2
= 𝐶1 + ( 2 )
𝑟 𝑟

𝜕𝑢𝜃 𝑢𝜃 𝐶2 𝐶2 𝐶2
− = 𝐶1 − ( 2 ) − 𝐶1 − ( 2 ) = −2 ( 2 )
𝜕𝑟 𝑟 𝑟 𝑟 𝑟
𝐶2 2𝐿𝑠 𝐶2
𝐶1 𝑅𝑖 + = 𝑤𝑖 𝑅𝑖 −
𝑅𝑖 𝑅𝑖2

7
 𝐶2 2𝐿𝑠 𝐶2
𝐶1 𝑅𝑒 + = 𝑤𝑒 𝑅𝑒 −
𝑅𝑒 𝑅𝑒2

1 2𝐿𝑠
𝐶1 = 𝑤𝑖 − 𝐶2 ( + )
𝑅𝑖2 𝑅𝑖3

1 2𝐿𝑠
𝐶1 = 𝑤𝑒 − 𝐶2 ( + )
𝑅𝑒2 𝑅𝑒3

𝑅𝑒 − 2𝐿𝑠 𝑅𝑖 − 2𝐿𝑠
𝑤𝑖 − 𝑤𝑒 + 𝐶2 ( − )=0
𝑅𝑒3 𝑅𝑖3

𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 − 2𝐿𝑠 (𝑅𝑖3 + 𝑅𝑒3 )

𝑤𝑖 − 𝑤𝑒 + 𝐶2 ( )=0
𝑅𝑖3 𝑅𝑒3
Se define C2

𝑅𝑖3 𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 )

𝐶2 = ( )=0
𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 − 2𝐿𝑠 (𝑅𝑖3 + 𝑅𝑒3 )
Se define C1
𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 )(𝑅𝑖 + 2𝐿𝑠 )
𝐶1 = 𝑤𝑖 −
𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 − 2𝐿𝑠 (𝑅𝑖3 + 𝑅𝑒3 )

𝑅𝑖3 𝑅𝑒 𝑤𝑖 − 𝑅𝑖 𝑅𝑒3 𝑤𝑒 − 2𝐿𝑠 (𝑤𝑖 𝑅𝑖3 + 𝑤𝑖 𝑅𝑒3 )

𝐶1 =
𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 − 2𝐿𝑠 (𝑅𝑖3 + 𝑅𝑒3 )

𝑅𝑖3 𝑅𝑒 𝑤𝑖 − 𝑅𝑖 𝑅𝑒3 𝑤𝑒 − 2𝐿𝑠 (𝑤𝑖 𝑅𝑖3 + 𝑤𝑖 𝑅𝑒3 ) 𝑅𝑖3 𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 ) 1

𝐶1 = 3 3 3 3 + 3 3 3 3
𝑅𝑒 𝑅𝑖 − 𝑅𝑖 𝑅𝑒 − 2𝐿𝑠 (𝑅𝑖 + 𝑅𝑒 ) 𝑅𝑒 𝑅𝑖 − 𝑅𝑖 𝑅𝑒 − 2𝐿𝑠 (𝑅𝑖 + 𝑅𝑒 ) 𝑟

Ls = 0

𝑅𝑖3 𝑅𝑒 𝑤𝑖 − 𝑅𝑖 𝑅𝑒3 𝑤𝑒 𝑅𝑖3 𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 ) 1

𝑢(𝑟) = 𝑟 +
𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 𝑅𝑒 𝑅𝑖3 − 𝑅𝑖 𝑅𝑒3 𝑟

8
 𝑅𝑖 𝑅𝑒 (𝑅𝑖2 𝑤𝑖 − 𝑅𝑒2 𝑤𝑒 ) 𝑅𝑖3 𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 ) 1
𝑢(𝑟) = 𝑟 +
𝑅𝑖 𝑅2 (𝑅𝑖2 − 𝑅𝑒2 ) 𝑅𝑖 𝑅𝑒 (𝑅𝑖2 − 𝑅𝑒2 ) 𝑟

𝑅𝑖2 𝑤𝑖 − 𝑅𝑒2 𝑤𝑒 ) 𝑅𝑖3 𝑅𝑒3 (𝑤𝑒 − 𝑤𝑖 ) 1

𝑢(𝑟) = 𝑟 +
𝑅𝑖 𝑅2 (𝑅𝑖2 − 𝑅𝑒2 ) 𝑅𝑖 𝑅𝑒 (𝑅𝑖2 − 𝑅𝑒2 ) 𝑟

𝑅𝑒2 𝑤𝑒 − 𝑅𝑖2 𝑤𝑖 𝑅𝑖2 𝑅𝑒2 (𝑤𝑖 − 𝑤𝑒 ) 1

𝑢(𝑟) = 𝑟+
𝑅𝑒2 − 𝑅𝑖2 𝑅𝑖2 − 𝑅𝑒2 𝑟

𝑅𝑒2 𝑤𝑒 − 𝑅𝑖2 𝑤1 𝑅𝑖2 𝑅𝑒2 (𝑤𝑖 − 𝑤𝑒 ) 1

𝑢(𝑟) = 𝑟 +
𝑅𝑒2 − 𝑅𝑖2 𝑅𝑖2 − 𝑅𝑒2 𝑟

4. Grafique el perfil de velocidad 𝑢𝜃. en el espacio entre los cilindros para valores de 𝐿𝑠
=0, 0.1, 1, 10 y 𝑟𝑒⁄𝑟𝑖 = 1.5, 2, 3.

5. Analice los resultados obtenidos y concluya respecto al comportamiento que ocurre al

variarambos parámetros. Cómo justifica el valor de velocidad mayor en el cilindro externo
estáticorespecto al cilindro interno en movimiento para determinados casos de 𝐿𝑠.
R/: En cuanto al efecto de ambos parámetros, ambos impactan en el comportamiento del
perfil de velocidades, empezando con el parámetro 𝐿𝑠 , el cual, siendo 0, no tiene efecto
alguno y el perfil de velocidades es similar a los ejemplos clásicos, sin embargo, e medida
que este parámetro aumenta, el perfil de velocidades comienza a volverse más uniforme,
9
 puesto que permite mayor movilidad del fluido cerca de las paredes, llegando al punto en
que incluso en la pared externa, la velocidad del flujo puede ser mayor que la velocidad en
la pared interna. En cuanto a la separación de los dos cilindros, cuando esta es pequeña, el
gradiente de velocidad es más pronunciado y el deslizamiento es muy notable cerca de las
paredes, sin embargo, a medida que esta separación incrementa, el perfil de velocidad
tiene una distribución más homogénea, y el deslizamiento se suaviza.
En cuanto a la velocidad en el cilindro externo, Con Ls grande, el deslizamiento permite
mayor movilidad del fluido cerca del cilindro externo, resultando en una velocidad mayor
en esta región. Esto suaviza el gradiente de velocidad y permite un efecto de "arrastre" que
incrementa la velocidad en el cilindro externo.

Referencias

J. Tu, G.-H. Yeoh, and C. Liu, “Chapter 3 - Governing Equations for CFD: Fundamentals,” in
Computational Fluid Dynamics (Third Edition), Third Edit., J. Tu, G.-H. Yeoh, and C. Liu, Eds.
Butterworth-Heinemann, 2018, pp. 65–124.

[Link]
3k8GsMtzOk&t=2358s&ab_channel=MiguelZoritaCalvo

[Link]

10