Informe Modelos de Turbulencia
Informe Modelos de Turbulencia
Informe Modelos de Turbulencia
Turbomáquinas I
Integrantes
Lima,2019
“Año de la lucha contra la corrupción y la impunidad”
INTRODUCCIÓN
DEFINICIONES
1. TURBULENCIA
Desde el punto de vista de la teoría de CFD, la turbulencia tiene que ver con tres
cosas:
Nuestra visión de nuestro mundo exterior,
Nuestra máquina digital, y
La naturaleza de nuestras propias facultades mentales.
“Año de la lucha contra la corrupción y la impunidad”
a) Escalas de turbulencia:
de la hipótesis del medio continuo), incluidos los flujos turbulentos. Sin embargo,
para un estudiante promedio de mecánica de fluidos esta premisa podría parecer
una paradoja porque las ecuaciones de N-S que él conoce son para flujo en
régimen laminar; lo que no es cierto. Es más, desde siempre, el estudio del flujo
turbulento se ha enfocado a partir de las ecuaciones del flujo laminar, tal como se
describe en las siguientes Secciones.
Por lo anterior, se sabe que las ecuaciones exactas que describen el movimiento
turbulento son las ecuaciones de Navier-Stokes8. Entonces, previamente
recordemos que con el término “ecuaciones de Navier-Stokes” (o “ecuaciones de
N-S”) nos referimos a la ecuación de conservación de momentum lineal, llamado
también ecuación de movimiento, y a la ecuación de conservación de masa.
Asimismo, conviene recordar que, aun cuando las ecuaciones de N-S son
bastante genéricas, pueden simplificarse notoriamente de acuerdo las condiciones
particulares del problema. Por ejemplo, la ecuación vectorial de continuidad para
flujo incompresible toma la siguiente forma sencilla:
𝛻⃗. 𝑉
⃗ =0
𝜕𝑢 𝜕𝑣 𝜕𝑤
+ 𝜕𝑦 + =0
𝜕𝑥 𝜕𝑧
⃗
𝐷𝑉
𝜌 𝐷𝑡 = −𝛻⃗𝑝 + 𝜇𝛻 2 𝑉
⃗ + 𝜌𝑔
Componente en x:
𝝏𝒖 𝝏𝒖 𝝏𝒖 𝝏𝒖 𝜹𝒑 𝝏𝟐 𝒖 𝝏𝟐 𝒖 𝝏𝟐 𝒖
𝝆 ( 𝝏𝒕 + 𝒖 𝝏𝒙 + 𝒗 𝝏𝒚 + 𝒘 𝝏𝒛 ) = − 𝜹𝒙 + 𝝁 (𝝏𝒙𝟐 + 𝝏𝒚𝟐 + 𝝏𝒛𝟐 ) + 𝝆𝒈𝒙
Componente en y:
𝝏𝒗 𝝏𝒗 𝝏𝒗 𝝏𝒗 𝜹𝒑 𝝏𝟐 𝒗 𝝏𝟐 𝒗 𝝏𝟐 𝒗
𝝆 ( 𝝏𝒕 + 𝒖 𝝏𝒙 + 𝒗 𝝏𝒚 + 𝒘 𝝏𝒛 ) = − 𝜹𝒚 + 𝝁 (𝝏𝒙𝟐 + 𝝏𝒚𝟐 + 𝝏𝒛𝟐 ) + 𝝆𝒈𝒚
“Año de la lucha contra la corrupción y la impunidad”
Componente en z:
𝝏𝒘 𝝏𝒘 𝝏𝒘 𝝏𝒘 𝜹𝒑 𝝏𝟐 𝒘 𝝏𝟐 𝒘 𝝏𝟐 𝒘
𝝆 ( 𝝏𝒕 + 𝒖 𝝏𝒙 + 𝒗 𝝏𝒚 + 𝒘 𝝏𝒛 ) = − 𝜹𝒛 + 𝝁 ( 𝝏𝒙𝟐 + 𝝏𝒚𝟐 + ) + 𝝆𝒈𝒛
𝝏𝒛𝟐
𝝏𝝓
⃗⃗⃗ 𝝓𝑽
+ (𝜵. ⃗⃗ ) = 𝝓𝑮 + 𝜹𝜵𝟐 𝝓
𝝏𝒕
Donde
𝜱
𝝓 = ∀ : concentración de la propiedad escalar 𝛷, siendo el volumen de un volumen
de control elemental,
⃗𝑽
⃗ : campo de velocidad (velocidad promedio de masa),
𝜕𝜙
⃗⃗⃗ 𝑉
+ 𝜙(𝛻. ⃗ . 𝛻⃗)𝜙 = 𝜙𝐺 + 𝛿𝛻 2 𝜙
⃗ ) + (𝑉
𝜕𝑡
𝜕𝜙
⃗ . 𝛻⃗)𝜙 = 𝛿𝛻 2 𝜙
+ (𝑉
𝜕𝑡
“Año de la lucha contra la corrupción y la impunidad”
𝜕𝜙 𝜕𝜙 𝜕𝜙 𝜕𝜙 𝜕 2𝜙 𝜕 2 𝜙 𝜕 2 𝜙
+𝑢 +𝑣 +𝑤 = 𝛿( 2 + 2 + 2 )
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧
Modelos de turbulencia
“Año de la lucha contra la corrupción y la impunidad”
Imagen con simulación con método RANS (izquierda) y con método LES (derecha)
de un flujo no estacionario en la estela de un cilindro.
que se emplea el método LES se observa cómo un vórtice de gran escala se
rompe en vórtices de menor escala. Este método de simulación numérica
permite conocer con mayor precisión cómo se entremezcla el fluido en la
estela porque utiliza esquemas numéricos menos difusivos. El método
RANS, por el contrario, obtiene una solución representativa del flujo
medio.
Pese a las claras ventajas que tiene el método LES frente al método RANS, el
método RANS en más barato desde el punto de vista de la computación, es
un sistema más robusto y más fácil de utilizar. Por estas razones RANS
es más utilizado en el ámbito de la industria y está disponible en
prácticamente todos los sistemas de CFD.
Ventajas:
“Año de la lucha contra la corrupción y la impunidad”
Es fácil de implementar
Proporciona una buena predicción de chorros, capas de mezcla, estelas y
capas límite
Está suficientemente validado
Inconvenientes:
Es incapaz de describir flujos con separación o recirculación
Sólo calcula propiedades medias y tensiones turbulenta
b) Modelo k-ε.
pretende corregir los defectos del método anterior, y permitir calcular flujos con
recirculación o separación. Se define la energía cinética turbulenta instantánea
como kI(t) = K + k, siendo
1
Energía cinética turbulenta media: k=2 ( 𝑈 2 +𝑉 2 +𝑊 2 )
1 ̅̅̅̅ ̅̅̅̅ ̅̅̅̅̅
Energía cinética turbulenta: k=2 (𝑈 ´2 +𝑉 ´2 +𝑊 ´2 )
𝑘 3/2 𝑘2
v= 𝑘1/2 ; l= ; 𝑢𝑡 = 𝑝. 𝑐𝑢 .
𝜀 𝜀
Ventajas
Inconvenientes:
“Año de la lucha contra la corrupción y la impunidad”
Para incluir las mejores características de ambos modelos Menter combina k-ε y k-
ω para formar un nuevo modelo de turbulencia de dos ecuaciones. Incorpora el
modelo k-ω en la región interior de la capa límite (hasta aproximadamente δ/2),
cambiando al modelo k-ε en la región exterior de la capa límite y la zona de la
estela, sin la necesidad de que el usuario intervenga. Las dos versiones de este
modelo son: El modelo base (Baseline model, BSL) y una versión modificada de
éste, el modelo del transporte de los esfuerzos de Reynolds (The shearstress
transport model, SST), desarrollado en 1994. Este último modelo incluye el
transporte de los esfuerzos de Reynolds en capas límites con gradientes de
presión adversos. En su inicio era un modelo muy utilizado en el ámbito de la
aeronáutica, pero actualmente se utiliza en la mayoría de los códigos del ámbito
industrial e incluso en la investigación.
Según los estudios realizados por Menter (12), ambos modelos proporcionan
resultados sin la dependencia que tenía con los valores especificados en la
corriente libre22 . Además, el modelo SST predice de forma adecuada las capas
límites con gradientes de presión adversos.
A continuación, se ofrece un estudio detallado de la mejora que introduce
Menter23 en el modelo de BSL para posteriormente explicar el modelo SST
presente el ANSYS Fluent®.
MODELO LES
La simulación completa de un flujo turbulento, dependiendo del tiempo y a
elevados número de Reynold, es prácticamente imposible por la finura de la
desratización que necesita. Dentro del amplio rango de escalas espaciales y
temporales de los flujos turbulentos. Las escalas más grandes Son las más
“Año de la lucha contra la corrupción y la impunidad”
Submodelo de Smagorinsky
𝜇𝐿𝐸𝑆 𝑐𝑝 𝜇𝐿𝐸𝑆
𝜇𝐿𝐸𝑆 = 𝜌(𝐶𝑆 ∆)2 |𝑆̅|; 𝑘𝐿𝐸𝑆 = ; (𝜌𝐷)𝑙,𝐿𝐸𝑆 =
𝑃𝑟 𝑆𝑐
“Año de la lucha contra la corrupción y la impunidad”
Este modelo está basado en el anterior, asumiendo que los coeficientes de las
escalas
modeladas son idénticos a los de las calculadas. Este submodelo es poco
disipativo y por
ello no se emplea casi nunca de forma independiente. No obstante, su empleo
conjunto con el submodelo de Smagorinsky reduce las inexactitudes provocadas
por la incertidumbre del coef. 𝐶𝑠 .
Características
“Año de la lucha contra la corrupción y la impunidad”
MODELO DNS
La "simulación numérica directa", a menudo abreviada DNS, es un término
utilizado principalmente en el contexto de la dinámica de fluidos computacional
(CFD).
Realizar DNS solo implica una discretización de estas ecuaciones y una solución
del sistema algebraico resultante de ecuaciones con una computadora. Nada más.
Sin embargo, estas ecuaciones son muy difíciles de resolver, debido a las
diferencias de escala extremas, tanto en tiempo como en espacio, que aparecen
en aplicaciones típicas. De hecho, las computadoras no son lo suficientemente
potentes como para poder simular con precisión con el flujo de DNS, por ejemplo,
un automóvil que se mueve rápidamente.
DNS EN SIMULACIÓN:
la malla se genera con vistas a asegurar que todas las escalas turbulentas se
puedan calcular directamente, sin la necesidad del modelado. Esta estrategia
implica en mallas con un número muy elevado de elementos, ya que es
impracticable para los casos reales.
“Año de la lucha contra la corrupción y la impunidad”
grandes para la
mayoría de las
aplicaciones
prácticas
CONCLUSIÓN
• Los errores que se podrán generar al aplicar este método sería al momento de
realizar las discretizaciones numéricas.
• Otro inconveniente de este método es que su utilización práctica queda limitada
a geometrías sencillas y a bajos número de Re.
BIBLIOGRAFÍA
https://upcommons.upc.edu/bitstream/handle/2099/10445/III%20-%20V24N3.pdf
https://www.sciencedirect.com/topics/engineering/direct-numerical-simulation
https://www.quora.com/What-is-the-main-theory-behind-direct-numerical-simulations
http://bibing.us.es/proyectos/abreproy/3718/fichero/Parte+I%252FCapitulo+2.pdf
“Año de la lucha contra la corrupción y la impunidad”
file:///C:/Users/ASUS/Downloads/Turbulencia%20(1).pdf
http://bibing.us.es/proyectos/abreproy/3718/fichero/Parte+I%252FCapitulo+2.pdf
https://matheo.uliege.be/bitstream/2268.2/3118/7/TFE_RubenGomez.pdf
http://www.bakker.org/dartmouth06/engs150/12-les.pdf
https://www.researchgate.net/figure/Different-approaches-of-turbulence-models-
regarding-the-model-scope-and-fields-of_fig3_326257568
“Año de la lucha contra la corrupción y la impunidad”