Universidad Técnica Federico Santa María

de Ingeniería Química y Ambiental

ILQ-230 Fenómenos de Transporte 2° Semestre 2018

ILQ230- MC3: FLUJO TURBULENTO

Cecilia Lobos – Javiera Acuña
Resumen
El objetivo principal del siguiente informe era realizar una búsqueda bibliográfica sobre los modelos de
turbulencia que ofrece el software Comsol Multiphysics, reconociéndose en que situaciones es óptima la
utilización de ciertos modelos, sus ventajas y desventajas, junto con clasificaciones de estos. Detallándose
todo en la siguiente sección. Además, se buscó visualizar las diferencias existentes entre algunos modelos,
específicamente entre los modelos 𝑘-𝜀 y 𝑘-𝜔 y luego el v2-f con los otros dos ya mencionados, en base a
la simulación de una sección con baffles y el nado de un tiburón, lográndose concluir que k- 𝜔 describe
de forma más exacta el comportamiento del fluido al pasar por una sección con baffles puesto que es más
preciso para describir flujos que poseen una curvatura fuerte. Por otra parte, para el nado de un tiburón el
modelo k-e se adapta de mejor forma a la geometría compleja y dirección en la que fluye el fluido. Se
comprueba además que el modelo v2-f no es muy útil para la simulación de flujos abiertos.

1. Fenomenología del problema resolución de problemas de turbulencia, los

cuales son (4):
Enfocándonos en el flujo turbulento tenemos que
este presenta fluctuaciones de velocidad en todas L-VEL y y-Plus Algebraico: Los modelos de
direcciones, por lo que la velocidad instantánea, turbulencia L-VEL y y-Plus algebraico calculan
𝑣, es una función que oscila irregularmente. Se la viscosidad del remolino usando expresiones
puede considerar que esta velocidad se puede algebraicas basadas solo en la velocidad del
descomponer en un promedio en el tiempo (𝑣̅ ) y fluido local y la distancia a la pared más cercana.
una componente fluctuante de la velocidad (𝑣′), No resuelven ninguna ecuación de transporte
es decir en: adicional. Estos modelos resuelven el flujo en
todas partes y son los menos precisos de los 9
𝑣 = 𝑣̅ + 𝑣′ existentes. Aun así, proporcionan buenas
Desarrollando la ecuación de continuidad y aproximaciones para el flujo interno,
movimiento en base a esta igualdad tenemos que especialmente en aplicaciones de refrigeración
estas resultan para un tiempo ajustado: electrónica.

∇ ∙ 𝑣̅ = 0 k-ε: El modelo k-ε resuelve dos variables: k, la

energía cinética de turbulencia; y ε (épsilon), la
𝐷𝑣̅ velocidad de disipación de la energía cinética de
𝜌 = −∇𝑝̅ − (∇ ∙ 𝜏̅ (4) ) − (∇ ∙ 𝜏̅ (6) ) + 𝜌𝑔
𝐷𝑡 turbulencia. Se utilizan funciones de pared en
Es posible notar que términos adicionales este modelo, por lo que el flujo en la región del
amortiguada no se simula. El modelo k-ε ha sido
aparecen en la segunda ecuación, 𝜏̅ (4) y 𝜏̅ (6) . La
históricamente muy popular para aplicaciones
magnitud 𝜏̅ (4) es término adicional que aparece
industriales debido a su buena tasa de
como consecuencia del ajuste de tiempo,
convergencia y sus requisitos de memoria
normalmente tabulado, el segundo término 𝜏̅ (6) relativamente bajos. No calcula de forma muy
generalmente se denominan esfuerzo de precisa los campos de flujo que muestran
Reynolds (1). Estas ecuaciones solo pueden gradientes de presión adversos, curvatura fuerte
resolverse directamente para casos simples de del flujo o flujo de chorro. Se desempeña bien
flujo. Es aquí donde entra Comsol Multiphysics, para problemas de flujo externo alrededor de
el cual entrega nueve modelos distintos para la geometrías complejas. Por ejemplo, el modelo k-

1
 Universidad Técnica Federico Santa María
de Ingeniería Química y Ambiental
ILQ-230 Fenómenos de Transporte 2° Semestre 2018

ε puede usarse para resolver el flujo de aire capa límite utilizando dos nuevas ecuaciones,
alrededor de un cuerpo. además de las dos ecuaciones para la energía
cinética de turbulencia (k) y la tasa de disipación
Reynolds bajo 𝑘-𝜀: La notación "bajo número de (ε). La primera ecuación describe el transporte de
Reynolds" no se refiere al flujo a escala global, fluctuaciones de velocidad turbulentas normales
sino a la región cercana a la pared donde a las líneas de corriente, mientras que la segunda
predominan los efectos viscosos, por lo que es un ecuación explica los efectos no locales, como la
modelo que reproduce correctamente el flujo a amortiguación inducida por la pared de la
medida que la distancia a la pared se acerca a redistribución de la energía cinética de
cero, es por esto, que si bien este modelo es turbulencia entre las direcciones normal y
similar al modelo k-ε, no necesita funciones de paralela. Las aplicaciones de este modelo van
pared, puede resolver el flujo en todas partes. Al más para flujos cerrados sobre superficies curvas,
ser una extensión del modelo k- ε comparte por ejemplo, para modelar ciclones.
muchas de sus ventajas y aplicaciones, pero
generalmente requiere una malla más densa, Realizable k-ε: La mayoría de los modelos de
puesto que al trabajar con bajo Reynolds tiende turbulencia incluyen restricciones de realización
como modelo a amortiguar la turbulencia. para garantizar la no negatividad de las tensiones
normales turbulentas, la desigualdad de Schwarz
El modelo Reynolds bajo k-ε puede calcular las entre cualquier cantidad fluctuante y para limitar
fuerzas de elevación y arrastre, y los flujos de la producción de turbulencia. Sin embargo, el
calor pueden modelarse con mayor precisión en modelo de turbulencia realizable k-ε aplica la
comparación con el modelo k-ε. realizabilidad al permitir que los coeficientes en
k-ω: El modelo k-ω es similar al modelo k-ε, pero las ecuaciones de transporte de turbulencia
resuelve para ω (omega), la tasa específica de varíen con respecto a la tasa de deformación del
disipación de la energía cinética. Es más no lineal flujo medio, así como k y ε. Esto da como
y, por lo tanto, más difícil de converger que el resultado un enfoque más suave y físico de los
modelo k-ε, y es bastante sensible a la estimación estados limitantes.
inicial de la solución. El modelo k-ω es útil en Los modelos anteriormente mencionados se
muchos casos en los que el modelo k-ε no es pueden agrupar las siguientes categorías según
preciso, como los flujos internos, flujos que las ecuaciones en las que se basan (3):
muestran una curvatura fuerte, flujos separados y
chorros. Un buen ejemplo de flujo interno es el Modelos de dos ecuaciones: k-ε, Realizable
flujo a través de una curva de tubería. k-ε, k-ω, SST, Bajo Re k-ε
SST: El modelo SST es una combinación del Modelos de ecuaciones de transporte: v2-f
modelo k-ε en el flujo libre y el modelo k-ω cerca Modelos de turbulencia algebraica: y-Plus
de las paredes. Es un modelo de número de algebraico y L-VEL
Reynolds bajo. Tiene requisitos de resolución Cabe destacar que hay ciertos modelos que ya no
similares al modelo k-ω y al modelo Reynolds se encuentran presentes en la versión 5.3 como el
bajo k-ε, pero su formulación elimina algunas modelo Spalart-Allmaras este es un modelo de
debilidades mostradas por los modelos puros k-ω número de Reynolds bajo que agrega una
y k-ε. Al ser una combinación de ambos modelos variable adicional para una viscosidad de
se podría considerar apto para aplicaciones remolinos cinemática no amortiguada. Este
industriales. modelo no calcula con precisión los campos que
v2-f: El modelo de turbulencia v2-f describe la exhiben flujo de corte o turbulencia en
anisotropía de la intensidad de turbulencia en la

2
 Universidad Técnica Federico Santa María
de Ingeniería Química y Ambiental
ILQ-230 Fenómenos de Transporte 2° Semestre 2018

descomposición. Su ventaja es que es bastante Tiburón: 𝑅𝑒 = 4,3 ⋅ 10C [−]

estable y muestra buena convergencia(2).
Todos los modelos nombrados corresponden a
modelos tipo RANS (Reynolds Averaged
Navier-Stokes), el software Comsol
Multiphysics no posee ningún modelo en base a
LES (Large Eddy Simulation) estos modelos
toman mucho tiempo incluso semanas para
obtener un resultado y son muchos más precisos.

2. Resultados:
Sección con Baffles: 𝑅𝑒 = 8,0 ⋅ 10> [−]

Figura 3 Líneas de Flujo del campo de velocidad para el modelo

k-e en el nado de un tiburón.

Figura 1 Líneas de Flujo del campo de velocidad para el modelo

k-e en una sección de baffles.

Figura 4 Líneas de Flujo del campo de velocidad para el modelo

k-w en el nado de un tiburón.

Figura 2 Líneas de Flujo del campo de velocidad para el modelo

k-w en una sección de baffles.

Figura 5 Líneas de Flujo del campo de velocidad para el modelo

v2-f en el nado de un tiburón.

3
 Universidad Técnica Federico Santa María
de Ingeniería Química y Ambiental
ILQ-230 Fenómenos de Transporte 2° Semestre 2018

3. Discusión de resultados figuras complejas. Lo antes mencionado se puede

comprobar claramente en las Figuras 3, 4 y 5, en
Sección con baffles:
donde se aprecia que el modelo k-e , figura 3,
Para este caso se utilizaron dos modelos de
describe notoriamente las turbulencias presentes
turbulencia para ejemplificar el comportamiento
en la aleta dorsal y la cola del tiburón. En cambio,
del fluido, en este caso agua, con un valor de
en el modelo k-w los vórtices son más sutiles y
Reynolds aproximado de 8,0 ⋅ 10> [−]. Los
en el modelo v2-f no se logra apreciar vórtices o
modelos utilizados fueron k-e, que considera la
flujos turbulentos.
energía cinética de turbulencia y la velocidad de
Aunque existen diferencias entre los modelos, se
disipación de la energía cinética en la
puede observar, tanto en la sección con baffles
turbulencia, y k-w, que considera la energía
como en el tiburón, que los modelos k-e y k-w
cinética de turbulencia y la tasa especifica de
muestran una buena aproximación al
disipación de la energía cinética.
comportamiento turbulento del fluido.
Se puede observar en la figura 2 claramente los
vórtices formados por el flujo turbulento, en
4. Conclusiones y recomendaciones
cambio, en la figura 1 los vórtices formados son
más sutiles. Esto es debido a que el primer Se pudo concluir, con ayuda de la herramienta
modelo mencionado no se adapta adecuadamente Comsol Multiphysics, que para una sección con
a los flujos de curvatura fuerte que son causados baffles un modelo k-w permite describir una
por los baffles. Por el contrario, se tiene que el buena aproximación sobre el comportamiento del
segundo modelo se adecua mejor a los flujos fluido frente a curvaturas fuertes y que para el
internos con desviaciones, comprobando así que nado de un tiburón el modelo k-e se adapta de
el modelo k-w describe de forma más exacta el mejor forma a la geometría compleja y dirección
comportamiento del fluido al pasar por una en que fluye el fluido. Se comprueba además que
sección con baffles o con cambio de diámetro. el modelo v2-f no es muy útil para la simulación
de flujos abiertos. Lo antes mencionado,
Tiburón: corrobora el comportamiento descrito de estos
Se aplicaron tres modelos distintos para analizar modelos en fenomenología.
el comportamiento del nado de un tiburón. Al Se recomienda que antes de simular alguna
igual que el caso anterior se utilizaron los determinada situación en el software realizar una
modelos k-e y k-w, agregando el modelo v2-f. búsqueda bibliográfica de los modelos
Para el análisis se utilizo un Reynolds turbulentos existentes, de tal manera conseguir
aproximado de 4,3 ⋅ 10C [−]. los mejores resultados.
A diferencia del diseño anterior, sección con
baffles, el tiburón se caracteriza por ser una 5. Referencias
figura compleja, por lo cual, conociendo las (1) Bird, R. B., Stewart, W. E., & Lightfoot, E. N. (2006).
BIRD: Fenómenos de transporte. (Vol. Segunda
características mencionadas en fenomenología edición). México: Limusa.
de cada modelo se puede prever que tanto el (2) Comsol. (s.f.). CFD Module. Obtenido de Comsol:
modelo k-w como el modelo v2-f no van a https://www.comsol.com/release/5.2/cfd-module
(3) Comsol. (s.f.). Simulate Fluid Flow Applications with
generar una buena aproximación para la the CFD Module. Obtenido de Comsol:
situación, debido a que el primero se caracteriza https://www.comsol.com/cfd-module
por modelar flujos internos o curvas fuertes, y el (4) Frei, W. (2017). Which Turbulence Model Should I
Choose for My CFD Application? Obtenido de Comsol:
segundo por modelar flujos cerrados sobre https://www.comsol.com/blogs/which-turbulence-
superficies curvas, situaciones que no van al model-should-choose-cfd-application/
caso. Por otro lado, el modelo k-e se caracteriza
por ser utilizado en describir flujos alrededor de