Generalization of the viscous stress tensor to the case of non-small gradients of hydrodynamic velocity: a path to numerical modeling of turbulence non-locality

Este trabajo generaliza el método de Chapman-Enskog para obtener una representación integral del tensor de esfuerzos viscosos para grandes gradientes de velocidad, lo que permite la modelación numérica de la no localidad de la turbulencia y fenómenos como las discontinuidades tangenciales que las formulaciones estándar de Navier-Stokes tienen dificultades para capturar.

Lo esencial

El Gran Problema: La "Batidora" vs. La "Tormenta"

Imagina que intentas predecir cómo se mueve un fluido (como el agua o el aire). Durante más de un siglo, los científicos han utilizado un famoso conjunto de reglas llamado las ecuaciones de Navier-Stokes. Estas reglas dependen de un ingrediente específico llamado el tensor de esfuerzos viscosos.

Piensa en este tensor como una "calculadora de fricción". En la versión estándar, asume que si empujas un fluido, la resistencia (fricción) que siente depende únicamente de qué tan rápido se mueve el fluido justo al lado del punto que estás observando. Es como asumir que si remueves una taza de café, la resistencia que sientes está determinada solo por las moléculas de café que tocan tu cuchara.

El Defecto:
El autor, A.B. Kukushkin, señala que esta "calculadora de fricción" estándar falla cuando las cosas se vuelven caóticas, como en una tormenta turbulenta o cuando dos corrientes de fluido se deslizan una junto a la otra a altas velocidades (discontinuidades tangenciales).

• La Analogía: Imagina una multitud de personas caminando por un pasillo. El modelo estándar asume que todos solo chocan con la persona inmediatamente adyacente. Pero en una multitud real (turbulencia), una persona podría ser empujada por alguien que está tres filas atrás, o una ola de movimiento podría atravesar toda la habitación. El modelo estándar ignora estas interacciones de "larga distancia".
• La Paradoja: Las matemáticas estándar también conducen a un resultado extraño: sugieren que si las partículas colisionan más a menudo (como en una niebla densa), el fluido debería fluir más fácilmente (menor viscosidad). Esto va en contra de nuestra intuición.

La Solución: Observar el Cuadro Completo

Kukushkin propone una nueva forma de calcular esta fricción. En lugar de mirar solo el vecindario inmediato, su nueva fórmula examina la historia y la ubicación completas del movimiento del fluido.

El Nuevo Enfoque:

1. Abandonar la Regla de los "Pequeños Pasos": Las matemáticas antiguas (método de Chapman-Enskog) solo funcionan si el fluido cambia muy lenta y suavemente. Kukushkin elimina esta regla. Permite cambios repentinos y bruscos en la velocidad, lo cual ocurre en la turbulencia real.
2. La Analogía del "Mensajero": En lugar de mirar solo el fluido en un punto, imagina que el fluido está lleno de pequeños "mensajeros" (vórtices o perturbaciones).
   • En el modelo antiguo, un mensajero solo habla con su vecino.
   • En el nuevo modelo de Kukushkin, un mensajero nace en un punto, cruza la habitación y entrega su mensaje a un punto lejano antes de detenerse.
3. La Fórmula Integral: Las nuevas matemáticas son una integral (una suma sobre una gran área). Calcula el esfuerzo (fricción) en un punto específico sumando los efectos de todos estos mensajeros que viajan desde cualquier otro lugar del fluido hacia ese punto.

Por Qué Esto Importa

1. Arreglando la "Paradoja":
Al permitir que estos mensajeros viajen largas distancias, la nueva fórmula corrige la extraña paradoja sobre la viscosidad. Explica por qué los fluidos se comportan como lo hacen incluso cuando las partículas colisionan con frecuencia. Los "vuelos largos" de los mensajeros explican la resistencia de una manera que el antiguo modelo de "pasos cortos" no podía.

2. Conectando con el Caos del Mundo Real (La Ley de Richardson):
El artículo menciona una famosa observación llamada la ley $t^3$ de Richardson.

• La Analogía: Si sueltas dos hojas en un río turbulento, el modelo estándar predice que se separarán lentamente (como $t^2$). Pero en la realidad, vuelan separadas mucho más rápido (como $t^3$).
• La Conexión: Este nuevo modelo de "larga distancia" explica naturalmente por qué las partículas se separan tan rápidamente. Los mensajeros viajan lejos y rápido, llevando la perturbación a través del fluido, lo cual coincide con la observación del mundo real de cómo se propaga la turbulencia.

3. Un Puente hacia Mejores Simulaciones por Computadora:
Actualmente, las simulaciones por computadora de la turbulencia a menudo deben usar "trucos" o números inventados porque las matemáticas estándar fallan en los bordes afilados (como donde un ala se separa del aire).

• La nueva fórmula de Kukushkin proporciona un puente matemático. Convierte los "trucos" en un cálculo riguroso basado en primeros principios. Permite que las computadoras modelen la turbulencia sumando estas interacciones de larga distancia, en lugar de simplemente adivinar.

Resumen en Poca Cosa

El artículo argumenta que la antigua forma de calcular la fricción de los fluidos es como intentar entender una conversación escuchando solo a la persona que está de pie junto a ti. Se pierde la imagen general.

Kukushkin ha escrito un nuevo reglamento que escucha a la habitación completa. Al tener en cuenta cómo las perturbaciones viajan a través de todo el fluido (no localidad), estas nuevas matemáticas:

• Arreglan una paradoja lógica sobre qué tan grueso o delgado debería ser un fluido.
• Explican por qué las partículas en una tormenta vuelan separadas tan rápidamente.
• Ofrecen un camino para que las computadoras simulen flujos complejos y caóticos (como el viento alrededor de un avión o el agua en una tubería) con mucha más precisión, sin necesidad de depender de conjeturas.

El autor señala que esta misma lógica podría eventualmente aplicarse a la transferencia de calor e incluso a la física del plasma, pero el logro central aquí es reescribir las reglas de la fricción de los fluidos para manejar la realidad "desordenada" de la turbulencia.

Resumen técnico

Resumen Técnico: Generalización del Tensor de Esfuerzos Viscosos para Gradientes No Pequeños

Planteamiento del Problema
El artículo aborda limitaciones fundamentales en la descripción estándar de la turbulencia hidrodinámica mediante las ecuaciones de Navier-Stokes, específicamente en lo que respecta al tensor de esfuerzos viscosos ($\sigma_{ik}$). El tensor convencional, derivado mediante el método de Chapman-Enskog a partir de la ecuación de Boltzmann, se basa en la suposición de gradientes pequeños de la velocidad hidrodinámica. El autor argumenta que esta suposición hace que el tensor estándar sea inadecuado para describir fenómenos de flujo complejos, como discontinuidades tangenciales, flujos separados y la no localidad inherente a la turbulencia.

Además, el artículo destaca una paradoja teórica en la derivación estándar: el coeficiente de viscosidad dinámica ($\eta$) resulta ser inversamente proporcional a la sección transversal de colisión de las partículas ($\sigma_{coll}$), en contra de las expectativas intuitivas. El autor atribuye esta paradoja a la restrictiva condición de "gradiente pequeño" del método de Chapman-Enskog, la cual se vuelve cada vez más factible a medida que aumenta la sección transversal de colisión, enmascarando así el comportamiento de vórtices localizados de larga duración.

Los modelos de turbulencia existentes (por ejemplo, Spalart-Allmaras, RANS) a menudo recurren a correcciones fenomenológicas o ecuaciones diferenciales locales que no logran capturar la naturaleza no local de la turbulencia, evidenciada empíricamente por la ley $t^3$ de Richardson para las correlaciones de pares en medios turbulentos. Si bien se han propuesto modelos de transporte no local (como vuelos/paseos de Lévy), a menudo carecen de una derivación rigurosa a partir de primeros principios dentro de la hidrogasdinámica.

Metodología
El autor propone una generalización del tensor de esfuerzos viscosos relajando la condición de gradientes pequeños de velocidad dentro del marco de Chapman-Enskog. La metodología procede de la siguiente manera:

1. Fundamento Cinético: La derivación comienza con la ecuación cinética de Boltzmann para un gas de un solo tipo. La función de distribución $f(\mathbf{p}, \mathbf{r}, t)$ se descompone en una parte de equilibrio termodinámico local (distribución de Maxwell $f^{(0)}$) y un término de corrección $\varphi$.
2. Relajación de las Restricciones de Gradiente: A diferencia del enfoque estándar, donde se desprecian las derivadas espaciales del término de corrección, el autor conserva las derivadas de $\varphi$ en el lado izquierdo de la ecuación de transporte. Esto reconoce que en flujos de corte y vórtices localizados, los gradientes del término de corrección pueden ser significativamente mayores que los de la distribución de equilibrio suave.
3. Formulación de la Ecuación de Transporte: Se utiliza un integral de colisión simplificado (tipo BGK), lo que conduce a una ecuación cinética de transferencia para $\varphi$. Esta ecuación presenta un término fuente dependiente del gradiente de la velocidad de masa promedio y un término sumidero que representa la absorción de portadores (inverso del camino libre).
4. Solución Integral: La ecuación de transporte se resuelve asumiendo que el movimiento de los portadores es significativamente más rápido que el movimiento del medio masivo. La solución para $\varphi$ se expresa como una integral sobre las trayectorias de los portadores desde el límite del medio hasta el punto de observación, incorporando un factor de decaimiento exponencial basado en la longitud del camino libre.
5. Derivación del Tensor de Esfuerzos: El tensor de esfuerzos viscosos se calcula integrando el producto del término de corrección, la distribución de equilibrio y el flujo de momento sobre el espacio de fases. Al transformar la integración de trayectorias en una integración de volumen sobre el medio, el autor deriva una representación integral no local del tensor de esfuerzos.

Contribuciones y Resultados Clave
El resultado principal es la derivación de un tensor de esfuerzos viscosos generalizado $\sigma_{ik}$ que es una integral sobre coordenadas espaciales en lugar de un operador diferencial local.

• Representación No Local: La nueva expresión para $\sigma_{ik}$ (Ecuación 2.13) depende explícitamente de los gradientes de velocidad en todos los puntos dentro del medio, ponderados por un núcleo que involucra la distancia entre puntos y la longitud del camino libre. Esta estructura refleja las teorías de transporte no local, como el modelo de Biberman-Holstein para la transferencia de radiación resonante.
• Recuperación de la Teoría Estándar: En el límite de longitudes de camino pequeñas (gradientes pequeños), la expresión integral se reduce al tensor de esfuerzos viscosos local estándar (Ecuación 1.2), con el coeficiente de viscosidad $\eta = P\tau$. Esto asegura la consistencia con la teoría establecida en el régimen apropiado.
• Resolución de la Paradoja: La forma generalizada resuelve la dependencia inversa de la viscosidad respecto a la sección transversal de colisión. El autor postula que la paradoja surge porque la derivación estándar asume gradientes pequeños, una condición incompatible con los vórtices localizados de larga vida que actúan como portadores de perturbaciones en flujos turbulentos.
• Leyes de Conservación: La derivación demuestra que la corrección a la función de distribución conserva el número total de partículas en el medio, representando un intercambio dinámico entre el medio principal y el conjunto de portadores de perturbaciones.

Significado y Afirmaciones
El artículo afirma que este trabajo constituye un paso hacia una teoría no local de la turbulencia hidrodinámica derivada de primeros principios.

• Modelado de Discontinuidades: El tensor generalizado proporciona un operador matemático capaz de describir discontinuidades tangenciales y flujos de separación donde las derivadas de la velocidad son ilimitadas, abordando una brecha de larga data en el modelado numérico.
• Conexión con la Ley de Richardson: La estructura no local del tensor derivado se presenta como una base necesaria para modelar numéricamente la no localidad de la turbulencia, específicamente para reproducir la ley empírica $t^3$ de Richardson para las correlaciones de pares.
• Validación del Enfoque de Kolmogorov: El autor sugiere que esta generalización puede explicar el éxito del razonamiento dimensional de Kolmogorov en la predicción de fenómenos de turbulencia, a pesar de que históricamente ha sido difícil derivar estos resultados directamente de la ecuación estándar de Navier-Stokes.
• Aplicabilidad Más Amplia: Se señala que la metodología es potencialmente aplicable a la transferencia de calor en gases y líquidos, así como a medios más complejos como el plasma, ofreciendo un enfoque unificado para fenómenos de transporte no local.

El artículo concluye que, si bien la expresión derivada es un avance teórico significativo, la demostración completa de su capacidad para reproducir la ley de Richardson requiere un modelado numérico adecuado posterior.