IGA-LBM: Isogeometric lattice Boltzmann method

Resumen Técnico: IGA-LBM: Método de Lattice Boltzmann Isogeométrico

Planteamiento del Problema
El Método de Lattice Boltzmann (LBM) es un marco mesoscópico potente para la dinámica de fluidos computacional, valorado por su paralelismo intrínseco y su sencillo tratamiento de fronteras. Sin embargo, el LBM tradicional depende de redes cartesianas uniformes, que aproximan las fronteras curvas mediante una representación de "escalonamiento" (stair-step). Como ilustra el artículo, esta discretización geométrica introduce errores significativos: las facetas dentadas de la frontera distorsionan las normales a la pared y la forma efectiva del cuerpo, propagando fuerzas espurias y degradando las predicciones de la capa límite. Aunque existen diversas estrategias de mitigación (p. ej., tratamientos basados en interpolación, métodos de celda cortada, técnicas de frontera inmersa), estas suelen comprometer la estabilidad numérica, introducir artefactos de interpolación o incurrir en costes computacionales prohibitivos. El desafío central es desarrollar un marco que preserve la simplicidad y eficiencia del LBM eliminando al mismo tiempo los errores de discretización geométrica en dominios curvos.

Metodología
Los autores proponen el Método de Lattice Boltzmann de Colocación Isogeométrica (IGA-LBM), el cual integra el Análisis Isogeométrico (IGA) con el sistema BGK de velocidades discretas. La metodología se aparta del modelo clásico de "transporte-colisión" (stream-collide) de dos maneras fundamentales:

1. Representación Geométrica Exacta: En lugar de una rejilla cartesiana, el dominio físico se parametriza utilizando funciones de base B-spline o NURBS (B-splines no uniformes racionales). Esto permite una representación analíticamente exacta de fronteras curvas (p. ej., círculos, aerofolios) directamente desde las definiciones CAD. Las funciones de distribución se aproximan utilizando las mismas funciones de base NURBS que la geometría (formulación isoparamétrica).
2. Colocación de Forma Fuerte: En lugar de transportar poblaciones a lo largo de enlaces de red, el método resuelve la ecuación Boltzmann-BGK de tiempo continuo en forma fuerte sobre la rejilla paramétrica ajustada al cuerpo (body-fitted).
   • Discretización Espacial: Las funciones de distribución se colocan en puntos de Greville. El término de advección transformado se evalúa utilizando operadores de diferenciación de alto orden. En el interior de rejillas uniformes, el operador de colocación de spline cúbico es algebraicamente idéntico a un esquema de diferencias finitas compactas (Padé) de cuarto orden.
   • Tratamiento de Fronteras: En las rejillas agrupadas hacia la pared, necesarias para resolver capas límite, la colocación de spline ingenua pierde precisión en la pared. Los autores reemplazan la derivada de colocación cerca de las fronteras con un operador de diferencias finitas de orden 5, consistente con la frontera, construido directamente sobre los nodos físicos. Esto restaura una precisión uniforme de cuarto orden hasta la pared.
   • Términos Métricos: Las velocidades de red contravariantes se computan una sola vez a partir de los términos métricos analíticamente exactos del mapeo NURBS. Esto asegura que los términos métricos estén libres de errores de diferenciación de rejilla.
   • Integración Temporal: El sistema semidiscreto se avanza utilizando un integrador Runge–Kutta de cuatro etapas (RK4) explícito.
   • Estabilización: Dado que el operador de colocación centrado no es disipativo, se aplica un filtro de paso bajo implícito de sexto orden, selectivo de escala, una vez por cada paso de tiempo para amortiguar las inestabilidades de escala de rejilla sin afectar el rango resuelto.
   • Condiciones de Frontera: El método emplea un enfoque de consistencia de momentos para las fronteras (paredes de no deslizamiento, campo lejano) en lugar de sobrescribir poblaciones, ya que la regla de rebote (bounce-back) por enlace es incompatible con la formulación de colocación.

Contribuciones Clave y Resultados Teóricos
El artículo proporciona un análisis teórico completo que sustenta el esquema:

• Límite Hidrodinámico: El esquema recupera las ecuaciones de Navier–Stokes incompresibles con viscosidad $\nu = c_s^2 \tau$. Notablemente, debido a que el método integra directamente la ecuación de tiempo continuo (en lugar de usar un transporte exacto), evita el artefacto clásico de corrección de viscosidad de $-\Delta t/2$.
• Preservación Exacta de Flujo Libre: El Teorema 7 demuestra que la discretización advectiva preserva un flujo uniforme con precisión de máquina en rejillas arbitrarias ajustadas al cuerpo sin requerir una ley de conservación geométrica discreta (GCL). Esto se debe a que los operadores de diferenciación anulan las constantes, independientemente de la precisión métrica.
• Orden de Precisión: Se demuestra que el operador de colocación en el interior es de cuarto orden de precisión (Lema 4). El Teorema 6 establece que el cierre de diferencias finitas en nodos físicos restaura la precisión uniforme de cuarto orden hasta la pared en rejillas agrupadas, corrigiendo la degradación del orden de la colocación ingenua.
• Exactitud Métrica: La Proposición 9 demuestra que el uso de términos métricos analíticamente exactos preserva el orden formal de la discretización de la derivada direccional a través del mapeo. Un estudio de ablación muestra que el uso de métricas diferenciadas de menor orden (incluso de cuarto orden) puede degradar silenciosamente el rendimiento del esquema si no coinciden perfectamente, mientras que la métrica NURBS exacta elimina esta fuente de error.
• Restricciones de Estabilidad: El Teorema 10 establece que el paso de tiempo explícito está limitado por el término de colisión estirado (stiff) a $\Delta t \lesssim 2.785 \tau$, lo que implica un escalamiento de $\Delta t = \mathcal{O}(\text{Re}^{-1})$ para parámetros de flujo fijos. Esta es una restrición necesaria para el sistema BGK no dividido (unsplit).

Resultados Numéricos y Validación
El método se valida contra varios casos de prueba:

• Pruebas de Flujo Libre: En rejillas curvas (cilindro, elipse, NACA0012), el método preserva un flujo libre uniforme hasta la precisión de máquina ($\sim 10^{-15}$), mientras que los métodos tradicionales de escalonamiento exhiben errores geométricos de $\mathcal{O}(1)$ en las normales a la pared.
• Vórtice de Taylor–Green: Los estudios de convergencia en 2D y 3D confirman los órdenes de precisión diseñados (2, 4 y 6) para el operador de advección. El método simula con éxito la transición a la turbulencia en el vórtice de Taylor–Green 3D.
• Cavidad Impulsada por Tapa (Re=100): El método reproduce los perfiles de velocidad de referencia (Ghia et al.) con alta precisión, limitado únicamente por el suelo de compresibilidad $\mathcal{O}(\text{Ma}^2)$ del modelo atérmico.
• Flujo alrededor de un Cilindro Circular:
  • Estacionario (Re=40): La longitud de recirculación y el coeficiente de arrastre coinciden con los datos de referencia dentro de un $\sim 9\%$ y $6\text{--}12\%$, respectivamente. Los autores atribuyen la ligera sobrepredicción del arrastre a la evaluación del esfuerzo en la pared y a los efectos de dominio finito, más que a la precisión en el interior.
  • No estacionario (Re=100): El número de Strouhal se predice con una precisión de $\sim 2.5\%$ respecto a los valores de referencia. La estela permanece estable hasta Re=1000.
• Aerofolio NACA0012 (Re=500, $\alpha=10^\circ$): El método reproduce las distribuciones de presión superficial y fricción superficial de referencia (Hafez et al.) con un error de pocos puntos porcentuales. Los coeficientes de sustentación y arrastre han convergido con la rejilla.

Significancia y Reivindicaciones
El artículo enmarca explícitamente la significancia de IGA-LBM con matices específicos:

1. Fidelidad Geométrica, No Superioridad Interior: Los autores enfatizan que la principal ventaja de la ruta isogeométrica es la geometría exacta y las métricas exactas, no un estencil interior superior. En rejillas suaves, el operador interior es equivalente a un esquema de diferencias finitas compactas estándar de cuarto orden. La ganancia es la eliminación de los errores de escalonamiento y de los errores de diferenciación métrica, lo que permite preservar el orden formal de precisión hasta la pared en fronteras curvas.
2. Preservación de Flujo Libre: Aunque la preservación exacta del flujo libre es una característica clave, el artículo señala que es una propiedad de la forma advectiva en cualquier rejilla suave y no requiere una GCL discreta. Es una característica de robustez más que una exclusividad isogeométrica.
3. Compromisos (Trade-offs): El método no es intrínsecamente conservativo (es un esquema de valores puntuales de forma fuerte) y requiere un parámetro de filtro ajustable para la estabilidad. La restricción del paso de tiempo explícito escala con $\text{Re}^{-1}$, lo que puede ser computacionalmente costoso para números de Reynolds elevados.
4. Alcance: El método se presenta como un LBM de alto orden y preservación de flujo libre para flujos con fronteras curvas construidos sobre geometría analítica exacta. No pretende resolver las limitaciones de los flujos separados no estacionarios de alto número de Reynolds sin un desarrollo adicional (p. ej., discretizaciones de forma débil conservativas).

En conclusión, el artículo presenta un LBM de alto orden y rigurosamente analizado que logra transferir los beneficios del análisis isogeométrico —específicamente la geometría exacta y la evaluación métrica— al solver cinético, resolviendo el problema de larga data de los errores de discretización geométrica en simulaciones de fluidos basadas en redes.