300px|derecha|miniatura|Simulación de grandes remolinos de un campo de velocidad de gas turbulento. La simulación de grandes remolinos (LES, por sus siglas en inglés) es un modelo matemático para la turbulencia utilizado en la dinámica de fluidos computacional. Fue propuesto inicialmente en 1963 por Joseph Smagorinsky para simular las corrientes de aire atmosférico, y explorado por primera vez por Deardorff (1970). Actualmente, la LES se aplica en una amplia variedad de aplicaciones de ingeniería, incluida la combustión, acústica, y simulaciones de la capa límite atmosférica.

La simulación de flujos turbulentos mediante la resolución numérica de las ecuaciones de Navier-Stokes requiere resolver una amplia gama de escalas de tiempo y longitud, todas las cuales afectan al campo de flujo. Esta resolución se puede lograr con la simulación numérica directa (DNS), pero la DNS es computacionalmente costosa y su coste impide la simulación de sistemas de ingeniería prácticos con geometrías o configuraciones de flujo complejas, como chorros turbulentos, bombas, vehículos y trenes de aterrizaje.

La idea principal detrás de la LES es reducir el coste computacional ignorando las escalas de longitud más pequeñas, que son las más costosas de resolver, mediante el filtrado de paso bajo de las ecuaciones de Navier-Stokes. Este filtrado de paso bajo, que puede considerarse como un promedio temporal y espacial, elimina eficazmente la información a pequeña escala de la solución numérica. Sin embargo, esta información no es irrelevante y su efecto sobre el campo de flujo debe modelarse, una tarea que es un área activa de investigación para problemas en los que las pequeñas escalas pueden desempeñar un papel importante, como los flujos cercanos a las paredes, flujos reactivos, y flujos multifásicos. Definición y propiedades del filtro

miniatura|derecha|250px|Un campo de velocidad producido por una simulación numérica directa (DNS) de turbulencia homogénea en descomposición. El tamaño del dominio es . miniatura|derecha|250px|El mismo campo de velocidad DNS filtrado con un filtro de caja y . miniatura|derecha|200px|El mismo campo de velocidad DNS filtrado con un Filtro (simulación de grandes remolinos)#Filtro de caja y . Se puede aplicar un filtro LES a un campo espacial y temporal y realizar una operación de filtrado espacial, una operación de filtrado temporal o ambas. El campo filtrado, denotado con una barra, se define como:

donde es el núcleo de convolución del filtro. Esto también se puede escribir como:

El núcleo del filtro tiene una escala de corte asociada y una escala de tiempo de corte . Las escalas menores que estas se eliminan de . Utilizando la definición de filtro anterior, cualquier campo puede dividirse en una parte filtrada y otra subfiltrada (denotada con un apóstrofo), como

Es importante señalar que la operación de filtrado de simulación de grandes remolinos no satisface las propiedades de un Operador de Reynolds.

Ecuaciones determinantes filtradas

Las ecuaciones determinantes de LES se obtienen filtrando las ecuaciones diferenciales parciales que determinan el campo de flujo . Existen diferencias entre las ecuaciones determinantes LES incompresibles y compresibles, lo que lleva a la definición de una nueva operación de filtrado.

Flujo incompresible

Para el flujo incompresible, se filtran la ecuación de continuidad y las ecuaciones de Navier-Stokes, lo que da como resultado la ecuación de continuidad incompresible filtrada,

y las ecuaciones de Navier-Stokes filtradas,

donde es el campo de presión filtrado y es el tensor de velocidad de deformación evaluado utilizando la velocidad filtrada. El término de advección filtrado no lineal es la principal causa de dificultad en el modelado LES. Requiere conocer el campo de velocidad sin filtrar, que se desconoce, por lo que debe modelarse. El análisis que sigue ilustra la dificultad causada por la no linealidad, es decir, que provoca una interacción entre escalas grandes y pequeñas, lo que impide la separación de escalas.

El término de advección filtrado puede dividirse, siguiendo a Leonard (1975), como:

donde es el tensor de tensión residual, de modo que las ecuaciones de Navier-Stokes filtradas se convierten en

con el tensor de tensión residual agrupando todos los términos no cerrados. Leonard descompuso este tensor de tensión como y proporcionó interpretaciones físicas para cada término. , el tensor de Leonard, representa las interacciones entre grandes escalas, , el término similar a la tensión de Reynolds, representa las interacciones entre las subescalas del filtro (SFS), y , el tensor de Clark, representa las interacciones entre escalas grandes y pequeñas. Modelar el término no cerrado es la tarea de los modelos de escala sub-grid (SGS). Esto resulta complicado debido a que el tensor de tensión subgrid debe tener en cuenta las interacciones entre todas las escalas, incluidas las escalas filtradas con escalas no filtradas.

La ecuación determinante filtrada para un escalar pasivo , como la fracción de mezcla o la temperatura, se puede escribir como

donde es el flujo difusivo de , y es el flujo subfiltrado para el escalar . El flujo difusivo filtrado no está cerrado, a menos que se le asuma una forma particular, como un modelo de difusión gradiente . se define de forma análoga a ,

y puede dividirse de forma similar en contribuciones de interacciones entre varias escalas. Este flujo de subfiltro también requiere un modelo de subfiltro.

Derivación Utilizando la notación de Einstein, las ecuaciones de Navier-Stokes para un fluido incompresible en coordenadas cartesianas son

Filtrando la ecuación del momento se obtiene

Si suponemos que el filtrado y la diferenciación son conmutativos, entonces

Referencias

Categoría:Dinámica de fluidos Categoría:Mecánica de fluidos Categoría:Ecuaciones en derivadas parciales