Este trabajo de investigación, realizado por Alikamp, explora la robustez de la frecuencia de desprendimiento de vórtices en simulaciones de flujo alrededor de un cilindro utilizando el método Lattice Boltzmann (LBM). El objetivo principal es reducir el costo computacional de estas simulaciones sin comprometer significativamente la precisión física. Tradicionalmente, simulaciones de alta resolución (DNS - Direct Numerical Simulation) son necesarias para capturar con precisión fenómenos como el desprendimiento de vórtices, pero esto es computacionalmente muy costoso. Este estudio demuestra que, sorprendentemente, la frecuencia de desprendimiento de vórtices es notablemente insensible a la resolución de la malla en simulaciones 2D a un número de Reynolds (Re) de 100.
En concreto, se observó que la frecuencia de Strouhal (St), un parámetro adimensional que caracteriza la frecuencia de desprendimiento, se mantuvo dentro del 2.5% de su valor de referencia (obtenido con una simulación de alta resolución) al reducir la resolución de la malla en un factor de 9 (de 320.000 a 35.511 celdas), lo que resultó en una aceleración de 37 veces en el tiempo de cálculo (wall time). El coeficiente de arrastre medio también se mantuvo dentro de los rangos esperados según la literatura. Esto implica que la física dominante del flujo en la estela (wake) está dictada por la geometría general y el número de Reynolds, y no por los detalles de la capa límite cercana a la superficie del cilindro.
Este descubrimiento tiene importantes implicaciones para diversas áreas. Permite el desarrollo de modelos de orden reducido más eficientes, facilita el uso de estrategias de mallado adaptativo (donde se concentra la resolución solo donde es necesaria para la precisión de la fuerza) y sugiere que los modelos sub-rejilla (sub-grid models) deberían enfocarse en la recuperación de la amplitud de la fuerza en lugar de la frecuencia, ya que esta última se conserva naturalmente con simulaciones de menor resolución. El trabajo se enmarca dentro del Kinetic-Pressure Ballooning Model (KPBM), un marco teórico que interpreta el desprendimiento de vórtices como una inestabilidad geométrica. El estudio proporciona una línea base DNS para evaluar el rendimiento del KPBM y su capacidad para mejorar la precisión de simulaciones de baja resolución.
Es importante tener en cuenta algunas limitaciones: el estudio se realizó en 2D y solo para un número de Reynolds de 100. La simulación presenta una relación de aspecto de 20:10 y una obstrucción del 10%, lo que afecta ligeramente el valor absoluto de St (18.7% por debajo del valor de la literatura). Además, no se abordó el estiramiento de vórtices en 3D ni se exploró la transición a un flujo turbulento a números de Reynolds más altos. El código fuente, los datos y los scripts de visualización están disponibles en GitHub, lo que facilita la reproducibilidad y la extensión de la investigación.