ESCUELA SUPERIOR POLITÉCNICA DEL LITORAL

FACULTAD DE INGENIERÍA EN MECÁNICA Y CIENCIAS DE LA PRODUCCIÓN

Asignatura: Reología de fluidos Paralelo: 101
Año académico: 2020 Término: II-PAO
Estudiantes: Fecha de entrega: 31/01/2021
• Cárdenas Emanuel
• Jorgge Netti
• Mata María
Profesor: Eduardo Castillo Orozco

PRÁCTICA DE LABORATORIO

SIMULACIÓN CFD DEL FLUJO DE UN FLUIDO NO NEWTONIANO DENTRO DE UNA

TUBERÍA
PERFIL DE VELOCIDAD, Y CAÍDA DE PRESIÓN Y POTENCIA REQUERIDA PARA
SUPERAR LA CAÍDA DE PRESIÓN

1. RESUMEN
Esta práctica tuvo como objetivo el desarrollar un modelo de Dinámica de Fluidos
Computacional para simular el flujo de un fluido incompresible no Newtoniano Ley de Potencia
dentro de una tubería El fluido utilizado fue el chocolate y por medio del programa ANSYS se
obtuvo el campo de velocidad del fluido dentro de una tubería de 0.25 m de diámetro y los
perfiles de velocidad en ambos regímenes experimentalmente. Por medio de las ecuaciones,
se calculó el número de Reynolds en régimen laminar y turbulento con 44.5384 y 5236.63, el
coeficiente de consistencia y el índice de comportamiento del fluido fue de0.5 menor por lo
que se clasifica como pseudoplástico. Se concluye que la caída de presión es mayor cuando
el chocolate se encuentra en un régimen turbulento con un valor de 570.21 Pa en
comparación a la caída de presión obtenida en el régimen laminar de 116.38 Pa.
2. OBJETIVOS
2.1 OBJETIVO GENERAL
Desarrollar un modelo de Dinámica de Fluidos Computacional para simular el flujo
de un fluido incompresible no Newtoniano Ley de Potencia dentro de una tubería.
2.2 OBJETIVOS ESPECÍFICOS
- Observar la gráfica de velocidad vs radio del chocolate para la comparación
del perfil de velocidad en los regímenes laminar y turbulento.
- Determinar la potencia y caída de presión de este flujo mediante una
simulación CFD, así como los valores necesarios para superar esta caída de
presión para compararlos con los resultados de simulación con
aproximaciones teóricas.
- Interpretar los resultados del comportamiento del chocolate en la simulación
numérica y los cálculos teóricos para el análisis del fluido en cada régimen.

3. INTRODUCCIÓN

Cuando un fluido cuya viscosidad varía con la temperatura y la tensión cortante que se le
aplica es considerado un fluido no newtoniano. Es por eso por lo que, un fluido no newtoniano
no tiene un valor de viscosidad constante y definido como un fluido newtoniano. El término
no newtoniano, generalmente se utiliza para clasificar todos los fluidos en dónde el esfuerzo
de corte no es directamente proporcional a la relación de corte. La viscosidad es una
manifestación del movimiento molecular dentro del fluido, en el cual las moléculas de regiones
con alta velocidad global chocan con las moléculas que se mueven con una velocidad global
menor, y viceversa. Por medio de estos choques se transporta la cantidad de movimiento de
una región de fluido a otra debido a que los movimientos son aleatorios y se clasifican como
fluidos independientes y dependientes del tiempo (Sánchez,2016).

Figura 1. Fluidos no newtonianos independientes del tiempo.

Figura 2. Fluidos no newtonianos dependientes del tiempo.

La Ley de la Potencia es uno de los modelos más utilizados para calcular la viscosidad de
los fluidos no newtonianos por medio de la relación entre el esfuerzo de corte y la velocidad
de corte como se ilustra en la ecuación 1.

Ec. 1
Dónde:
K= coeficiente de consistencia.
n= índice de comportamiento del flujo.
Cuando n>1 hace referencia a un fluido dilatante. Si n= 1, el fluido tendrá un comportamiento
newtoniano, sin embargo, si n<1 el fluido es pseudoplástico. En la práctica, el fluido utilizado
es el chocolate y su n= 0.5 por lo que es un fluido pseudoplástico. Para el flujo de este tipo
de fluidos dentro de una tubería, se puede clasificar en régimen de flujo laminar y flujo
turbulento mediante el número de Reynolds. Cuando NGRe< 2200 pertenece al flujo laminar
mientras que cuando NGRe≥ 4000 el fluido se encuentra en un régimen completamente
turbulento (Sánchez,2016).

Ec.2
Dónde:
ρ= densidad del fluido.
V= velocidad promedio del flujo.
D= diámetro interno de la tubería
Por medio de la ecuación 3, se puede calcular el perfil de velocidad del flujo de un fluido
modelado por la Ley de Potencia encontrado en un régimen laminar:

Ec. 3

4. MATERIALES Y MÉTODOS
4.1 MATERIALES
Se utilizaron dos programas:
- ANSYS
- Microsoft Excel
El software ANSYS sirvió para crear la geometría y el mallado del cilindro para la
representación del dominio del flujo del chocolate en la tubería. Al momento de extraer los
resultados de las magnitudes de velocidad en función del radio del dominio, se exportaron los
datos al programa Microsoft Excel para hacer la gráfica correspondiente.
4.2 PROCEDIMIENTO

Tabla 1. Datos del fluido no newtoniano utilizado (chocolate)

 Nombre del fluido Chocolate
PROPIEDADES DEL FLUIDO
Densidad a 20°C y 1
1200 kg/m3
atm
Viscosidad
Temperatura (T) 303 K
Índice de
comportamiento del 0.5
flujo (𝜂)
Coeficiente de
0.7
consistencia (K)
Viscosidad del
2.8 n (Pa*s)
chocolate a 80°C
Punto de fusión 45-50 °C
Punto de ebullición 33°C
Tamaño de la tubería 1.18 m3
VELOCIDAD PROMEDIO DEL FLUIDO EN EL RÉGIMEN (V)
Laminar 0.075 m/s
Turbulento 1.8 m/s

Tabla 2. Datos del perfil de velocidad vs el radio del chocolate en un régimen turbulento
experimental
Velocidad [ m s^-1 ] Y[m]
3,24E+07 -5,00E+06
3,27E+07 -3,89E+06
3,27E+07 -2,78E+06
3,28E+07 -1,67E+06
3,28E+07 -5,56E+05
3,28E+07 5,56E+05
3,28E+07 1,67E+06
3,27E+07 2,78E+06
3,27E+07 3,89E+06
3,25E+07 5,00E+06

Tabla 3. Datos del perfilde velocidad vs el radio del chocolate en un régimen laminar
experimental

Velocidad [ m s^-1 ] Y[m]

1,12E+07 -5,00E+06
1,18E+07 -3,89E+06
1,20E+07 -2,78E+06
1,23E+07 -1,67E+06
1,25E+07 -5,56E+05
1,24E+07 5,56E+05
1,21E+07 1,67E+06
1,19E+07 2,78E+06
1,16E+07 3,89E+06
1,11E+07 5,00E+06
Tabla 4. Datos del perfil de velocidad vs el radio del chocolate en un régimen turbulento
teórico.

Velocity [ m s^-1 ] Y[m]

3,54E+07 -5,00E+06
3,57E+07 -3,89E+06
3,57E+07 -2,78E+06
3,58E+07 -1,67E+06
3,58E+07 -5,56E+05
3,58E+07 5,56E+05
3,58E+07 1,67E+06
3,57E+07 2,78E+06
3,57E+07 3,89E+06
3,55E+07 5,00E+06

Tabla 5. Datos del perfil de velocidad vs el radio del chocolate en un régimen laminar
teórico.

Velocity [ m s^-1] Y [ m]
1,32E+07 -5,00E+06
1,38E+07 -3,89E+06
1,40E+07 -2,78E+06
1,43E+07 -1,67E+06
1,45E+07 -5,56E+05
1,44E+07 5,56E+05
1,41E+07 1,67E+06
1,39E+07 2,78E+06
1,36E+07 3,89E+06
1,31E+07 5,00E+06

5. RESULTADOS Y DISCUSIÓN
5.1 RESULTADOS
Laminar
Figura 2. Perfil de velocidad vs el radio del chocolate en un régimen laminar experimental.

Figura 3. Perfil de velocidad vs el radio del chocolate en un régimen laminar teórico.

Figura 4. Perfil de velocidad vs el radio del chocolate en un régimen laminar ANSYS-
teórico

Figura 5. Campo de velocidad del chocolate dentro de la tubería en régimen laminar.

Turbulento

Figura 6. Perfil de velocidad del chocolate en un régimen turbulento experimental.

Figura 7. Perfil de velocidad vs el radio del chocolate en un régimen turbulento teórico.

Figura 8. Perfil de velocidad vs el radio del chocolate en un régimen turbulento ANSYS-
teórico
 Figura 9. Campo de velocidad del chocolate dentro de la tubería en régimen turbulento.

5.2 RESULTADOS TEÓRICOS

Considerar el flujo de chocolate dentro de una tubería de 0.25 m de diámetro interior
y 6 m de longitud. El fluido fluye a una velocidad de 0.075 m/s. Determinar la caída de
presión y la potencia requerida para superar esta caída de presión.
Datos
D= 0.25 m
L= 6 m
V laminar = 0.075 m/s
V turbulento= 1.8 m/s
Fluido: chocolate (comportamiento pseudoplástico)
ρ = 1200 Kg/m3
n = 0.5
K= 0.7
Solución
Régimen laminar

𝐷 𝑛 ∗ 𝑉 2−𝑛 ∗ ρ 4𝑛 𝑛
𝑁𝐺𝑅𝐸 = ( ) = 8.24
K8n−1 3𝑛 + 1

0.5
(0.25)0.5 (0.075)2−0.5 (1200) 4(0.5)
𝑁𝐺𝑅𝐸 = 0.5−1
( )
(0.7)8 3(0.5) + 1
𝑁𝐺𝑅𝐸 = 44.5384
a) Factor de fricción de Fanning
16 16
fF = = = 0.3592
NGR 44.5384

b) Caída de presión y potencia requerida

2L 2
∆p = fF ρV
D
2(6)
∆p = 0.3592 (1200)(0.075)2 = 116.38 𝑃𝑎
0.25

D2 π
m = ρV
4
(0.25)2 π
m = (1200)(0.075) = 4.4178 𝐾𝑔/𝑠
4

m (116.38)(4.4178)
W = ∆P = = 0.4284 𝑊𝑎𝑡𝑡
ρ 1200
 Régimen turbulento

𝐷 𝑛 ∗ 𝑉 2−𝑛 ∗ ρ 4𝑛 𝑛
𝑁𝐺𝑅𝐸 = n−1 ( ) = 8.24
K8 3𝑛 + 1

0.5
(0.25)0.5 (1.8)2−0.5 (1200) 4(0.5)
𝑁𝐺𝑅𝐸 = 0.5−1 ( )
(0.7)8 3(0.5) + 1

𝑁𝐺𝑅𝐸 = 5236.63
c) Factor de fricción de Fanning
16 16
fF = = = 3.05𝑥 10−3
NGR 5236.63

d) Caída de presión y potencia requerida

2L 2
∆p = fF ρV
D
2(6)
∆p = 3.05𝑥 10−3 (1200)(1.8)2 = 570.21 𝑃𝑎
0.25

D2 π
m = ρV
4
(0.25)2 π
m = (1200)(1.8) = 106.03 𝐾𝑔/𝑠
4

m (570.21)(106.03)
W = ∆P = = 50.38 𝑊𝑎𝑡𝑡
ρ 1200

Cálculo de error de la velocidad en el régimen laminar

∣ Valor experimental − valor teórico ∣
𝑬𝒓𝒓𝒐𝒓 𝒂𝒃𝒔𝒐𝒍𝒖𝒕𝒐 =
valor teórico
∣ 11174967,1 − 13274967,1 ∣
𝑬𝒓𝒓𝒐𝒓 𝒂𝒃𝒔𝒐𝒍𝒖𝒕𝒐 = = 0.15819
13274967,1
𝑬𝒓𝒓𝒐𝒓 𝒑𝒐𝒓𝒄𝒆𝒏𝒕𝒖𝒂𝒍 = 100% ∗ 0.15819 = 15.89 %

Error
Y [ m]
porcentual
16% -5,00E+06
14% -3,89E+06
91% -2,78E+06
91% -1,67E+06
14% -5,56E+05
99% 5,56E+05
 99% 1,67E+06
𝑛 91% 2,78E+06
1
𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒍𝒂𝒎𝒊𝒏𝒂𝒓 = ∑ 𝑎1 91% 3,89E+06
𝑛
𝑖=1 92% 5,00E+06
16 + 14 + 91 + 91 + 14 + 99 + 99 + 91 + 91 + 92
𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒍𝒂𝒎𝒊𝒏𝒂𝒓 =
10
𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒍𝒂𝒎𝒊𝒏𝒂𝒓 = 70%

Cálculo de error de la velocidad en el régimen turbulento

∣ Valor experimental − valor teórico ∣
𝑬𝒓𝒓𝒐𝒓 𝒂𝒃𝒔𝒐𝒍𝒖𝒕𝒐 =
valor teórico
∣ 3,24E + 07 − 3,54E + 07 ∣
𝑬𝒓𝒓𝒐𝒓 𝒂𝒃𝒔𝒐𝒍𝒖𝒕𝒐 = = 0.08469
3,54E + 07
𝑬𝒓𝒓𝒐𝒓 𝒑𝒐𝒓𝒄𝒆𝒏𝒕𝒖𝒂𝒍 = 100% ∗ 0.08469 = 8.469 %

Error Porcentual Y [ m]
8,469% -5,00E+06
8,410% -3,89E+06
8,398% -2,78E+06
8,386% -1,67E+06
8,375% -5,56E+05
8,380% 5,56E+05
8,390% 1,67E+06
8,401% 2,78E+06
8,411% 3,89E+06
8,460% 5,00E+06

𝑛
1
𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒕𝒖𝒓𝒃𝒖𝒍𝒆𝒏𝒕𝒐 = ∑ 𝑎1
𝑛
𝑖=1

𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒕𝒖𝒓𝒃𝒖𝒍𝒆𝒏𝒕𝒐
8.469 + 8.410 + 8.398 + 8.386 + 8.375 + 8.380 + 8.390 + 8.401 + 8.411 + 8.460
=
10
𝑬𝒓𝒓𝒐𝒓𝒑𝒓𝒐𝒎 𝒕𝒖𝒓𝒃𝒖𝒍𝒆𝒏𝒕𝒐 = 8.408%

5.3 DISCUSIÓN
Para poder analizar y estudiar el comportamiento del chocolate dentro de una tubería se hizo
uso del software ANSYS con la que se pudo representar el perfil de velocidad y el campo de
velocidad del flujo dentro de la tubería. En la tabla 1, se detallan las propiedades de este
fluido, el tamaño de la tubería que fue de 1.18 m3. En las tablas posteriores, se especifican y
reconocen como varia la velocidad de este fluido a distintas velocidades que presenta la
tubería según el flujo del chocolate que teóricamente fue clasificado como laminar y turbulento
con 44.5384 y 5236.63 al calcular el número de Reynolds. En la figura 4 y 8 se observa el
comportamiento de este fluido dentro de la tubería en el cual tiene mayor velocidad al inicio
de la tubería, pero al pasar el fluido dentro de la misma, su velocidad va decreciendo. El
coeficiente de consistencia en el régimen laminar es de 106.03 Kg/s, valor mucho mayor en
comparación al coeficiente del régimen turbulento con un valor de 4.4178 Kg/s ya que se
incrementa la velocidad. El chocolate es considerado un fluido pseudoplástico ya que su
índice de comportamiento es de 0.5 menor al valor 1. Así mismo, la caída de presión es mayor
cuando el chocolate se encuentra en un régimen turbulento con un valor de 570.21 Pa en
comparación a la caída de presión obtenida en el régimen laminar de 116.38 Pa.

En la figura 3 y 7, se observa claramente como difiere experimentalmente el régimen laminar

y turbulento entre el perfil de velocidad vs. Radio del chocolate. El perfil de velocidad en flujo
laminar totalmente desarrollado en una tubería es parabólico con un máximo en línea central
y el mínimo -1,67E+06 en la pared de la tubería como se indica en la figura 3. Cuando este
fluido presenta un comportamiento laminar, el perfil de velocidad es representado con una
parábola o curva no simétrica debido que las fuerzas al ser aplicadas pueden hacer que se
vuelva el fluido más grueso por la viscosidad en el cual las moléculas de regiones con alta
velocidad global chocan con las moléculas que se mueven con una velocidad global menor,
y viceversa como se muestra en la figura 4. En la tabla 5 por medio de la ecuación 3, se
obtuvieron los valores teóricos de la velocidad y se graficó el perfil de velocidad del chocolate
en el régimen laminar. Se puede visualizar que la curva hecha por ANSYS en el primer punto
del radio interno, la velocidad respectiva es 3,24E+07. Sin embargo, cuando se calculan los
valores teóricos en el primer punto la velocidad obtenida fue de 1,32E+07 como se ilustra en
la figura 5. Por lo que se tiene un porcentaje de error de 70% del comportamiento del
chocolate en el régimen laminar.

Cuando se trata de un flujo turbulento las fluctuaciones son aleatorias y rápidas en las
regiones giratorias en este caso tenemos el chocolate, y estos giros transportan masa,
cantidad de movimiento y energía a otras regiones del flujo con mayor rapidez. El perfil de
velocidad cuando se encuentra en un régimen turbulento se basó en el análisis y mediciones
empíricas como se ilustra en la figura 7. Por otra parte, cuando el flujo es turbulento
predominan las fuerzas inerciales ante las viscosas por lo que el movimiento de las partículas
es más desordenado y como se observa en la figura 7 y 8 la forma de la parábola es más
aplanada. Se puede observar que la curva hecha por ANSYS en el primer punto del radio
interno, la velocidad respectiva es 3,24E+07. Sin embargo, cuando se calculan los valores
teóricos en el primer punto la velocidad obtenida fue de 3,54E+07 como se visualiza en la
figura 9. Por lo que se tiene un porcentaje de error de 8.408 % del comportamiento del
chocolate en el régimen turbulento.

6. CONCLUSIONES Y RECOMENDACIONES
6.1 CONCLUSIONES
− Se simulo el flujo del chocolate dentro de una tubería de 0.25m de diámetro y
6 m de longitud. Con los datos obtenidos en ANSYS y los valores teóricos
calculados, se observó que el régimen turbulento presenta una velocidad
teórica y experimental mayor que la del flujo laminar.
− La caída de presión para el régimen laminar fue menor que la del régimen
turbulento. Esto se debe a que la velocidad del régimen turbulento es mayor y
por tanto la resistencia de fricción al flujo del chocolate cuando presenta este
 comportamiento también debe ser mayor. Por lo que si se desea superar la
caída de presión del régimen laminar debe ser mayor a 570.21 𝑃𝑎.
− La potencia del fluido calculada teóricamente es menor que la potencia que la
obtenida a partir de la simulación realizada en el software ANSYS, esto indica
que para que el fluido pueda superar la caída de presión del sistema debe
tener una potencia mayor o igual a la potencia calculada teóricamente, para
que así no exista ningún error en el proceso.

6.2 RECOMENDACIONES
- Se recomienda ubicar los nombres a las superficies de la estructura ya que se evitarán
errores en el momento en que las condiciones establecidas estén mezcladas al
momento de realizar el diseño estructural.
- Ubicar la entrada del flujo como punto de origen para evitar confusiones al momento
de obtener resultados cuando se realice la estructura de la tubería.
- Las mallas encontradas alrededor de las paredes laterales deben ser finas ya que se
evita la existencia de gradientes altos impidiendo el deslizamiento del fluido interno.

7. BIBLIOGRAFÍA
Carrero, J. A & Rengel, J. E. (2012). Influencia del número de Reynolds sobre el campo
de flujo alrededor de dos cilindros. Universidad de los Andes, 33, 1-11.

Sánchez, L. (2016). Clasificación de los fluidos no newtonianos. Universidad Nacional

Autónoma de México, 1,1-5.

8. ANEXOS
ANEXO A. TABLA DE APORTACIONES
Integrante Aportación
Introducción
Cálculo de errores
Emanuel Cárdenas Recomendaciones
Simulación en ANSYS
Discusión
Cálculos teóricos de régimen turbulento
Graficas de régimen turbulento
Netti Jorgge Discusión
Conclusiones
Resultados teóricos
Cálculos teóricos de régimen laminar
Resumen
Graficas de régimen laminar
María Mata
Discusión
Objetivos específicos
Materiales y métodos