Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

Resumen Técnico: Reconstrucción Compatible con la Entropía para Flujos Viscoelásticos de Alto Número de Weissenberg

Planteamiento del Problema
El problema del alto número de Weissenberg sigue siendo un obstáculo central en la simulación de flujos viscoelásticos. En los modelos de tensor de conformación, el estado polimérico se representa mediante un tensor simétrico definido positivo (SPD) $A$. Si bien la positividad es esencial para el significado físico de la energía libre elástica ($\log \det A$) y el logaritmo de la matriz, el artículo sostiene que la positividad por sí sola es insuficiente para la estabilidad numérica en números de Weissenberg elevados.

Los métodos existentes, como la conformación logarítmica ($A = \exp \Psi$) y las reconstrucciones de raíz cuadrada ($A = BB^T$), logran imponer la restricción SPD con éxito. Sin embargo, los experimentos numéricos muestran que estos métodos aún pueden producir sobreoscilaciones de tensión no físicas y un crecimiento no físico de la energía elástica. El problema central identificado es la desconexión entre la restricción de cono (positividad) y la restricción de balance acoplado (compatibilidad de la energía libre). Específicamente, si el tensor utilizado en el trabajo de tensión de la ecuación de cantidad de movimiento difiere del tensor utilizado en el trabajo de estiramiento o en las variables de entropía de la ecuación de conformación, surge un defecto de trabajo con signo indefinido. Además, las reconstrucciones estrictamente positivas pueden introducir tres mecanismos no físicos específicos:

1. Sesgo de entropía de tipo Jensen: Los mapas no lineales (como la exponencial) aplicados a perturbaciones de media cero en el espacio logarítmico generan un sesgo de media positiva en el tensor físico y en la entropía.
2. Amplificación exponencial: En estados altamente estirados, las pequeñas perturbaciones logarítmicas son amplificadas por la exponencial de la matriz, escalando con el autovalor más grande $\lambda_{\max}(A)$.
3. Defectos de trabajo de signo indefinido: El uso de tensores reconstruidos discordantes para los términos de tensión y entropía impide la cancelación exacta del trabajo polimérico, lo que conduce a una inyección o eliminación incontrolada de energía.

Metodología
El artículo propone una reconstrucción logarítmica corregida y compatible con la entropía. En lugar de reemplazar los métodos de preservación de la positividad existentes, este enfoque actúa como un filtro de admisibilidad aplicado a una reconstrucción logarítmica de alto orden bruta antes de que entre en el balance acoplado de tensión-trabajo.

La metodología consiste en:

• Definición de la Trayectoria Logarítmica: Dado un predictor físico $\hat{A}$ y una reconstrucción de alto orden bruta $\tilde{A}$, se define una trayectoria en el espacio logarítmico: $A(\theta) = \exp(\hat{\Psi} + \theta(\tilde{\Psi} - \hat{\Psi}))$ para $\theta \in [0, 1]$.
• Restricción del Presupuesto de Entropía: El método busca el $\theta$ más grande tal que la energía libre discreta total (entropía) no exceda la entropía del predictor físico más un presupuesto local $\tau_h$. El presupuesto se escala para coincidir con el error de truncamiento esperado del esquema de alto orden (por ejemplo, $O(h^{2k+2})$).
• Algoritmo de Bisección: El $\theta$ óptimo se encuentra mediante una búsqueda de bisección local, por celda. Si la reconstrucción bruta satisface el presupuesto, $\theta=1$ (sin cambios). De lo contrario, el método retrocede a lo largo de la trayectoria hasta el punto máximo admisible.
• Imposición de Consistencia: El tensor seleccionado se utiliza de manera consistente en todos los términos de acoplamiento: trabajo de tensión de cantidad de movimiento, estiramiento de conformación, variables de entropía y cuadratura.

Contribuciones Clave y Resultados Teóricos
El artículo proporciona un análisis riguroso de este principio de reconstrucción, demostrando varias propiedades clave:

• Existencia y Convexidad: El funcional de entropía a lo largo de la trayectoria logarítmica es convexo, lo que garantiza la existencia de un parámetro admisible máximo $\theta^*$.
• Mínimo Amortiguamiento: Entre todos los estados admisibles en la trayectoria, el $\theta^*$ seleccionado minimiza la distancia a la reconstrucción bruta, haciendo que la corrección sea de "mínimo amortiguamiento".
• Control Espectral: La corrección preserva la independencia de marco y mantiene los límites espectrales derivados de los logaritmos de los extremos.
• Inactividad Asintótica: Si la reconstrucción de alto orden bruta ya es entropía-admisible (es decir, su defecto está dentro de la escala del error de discretización), la corrección se vuelve inactiva ($\theta^* = 1$), preservando la precisión formal de alto orden.
• Ventaja de Resolución en Alto Estiramiento: En regímenes de alto estiramiento, el método mantiene una alta precisión en métricas de log-relativo y entropía, mientras que las reconstrucciones directas en el espacio físico sufren problemas de escalamiento de derivadas proporcionales a $\lambda_{\max}(A)$.
• Estimación de Energía: Los autores demuestran una estimación de energía libre discreta para el sistema acoplado, mostrando que el presupuesto de reconstrucción $\tau_h$ actúa como una fuente de entropía controlada y acotada en lugar de una viscosidad artificial.

Resultados Numéricos
El artículo valida la teoría mediante una serie de diagnósticos:

• Pruebas Escalares: Se demostró que las perturbaciones logarítmicas de media cero crean un sesgo de entropía positivo (sesgo de Jensen) y que las reconstrucciones lineales, aunque de bajo sesgo, no preservan la positividad.
• Pruebas de Amplificación: Se confirmó que los defectos de fuerza de tensión en las reconstrucciones logarítmicas brutas crecen linealmente con el estiramiento de fondo $\lambda_{\max}(A)$.
• Verificación del Defecto de Trabajo: Se mostró que el uso de tensores discordantes para la tensión y la entropía conduce a defectos de energía de signo indefinido.
• Bucles de Alto Número de Weissenberg Dinámicos: En simulaciones con $Wi=80$ y alto estiramiento, las reconstrucciones logarítmicas estándar acumularon errores de entropía y de fuerza de tensión con el tiempo, mientras que la reconstrucción corregida controló estos defectos, reduciendo los errores de fuerza en órdenes de magnitud.
• Benchmarks de Resolución: El método corregido mantuvo la compatibilidad de la entropía y controló los defectos de fuerza a través de refinamientos de malla y distintos niveles de estiramiento, mientras que los métodos brutos mostraron errores crecientes al aumentar el estiramiento.

Significado y Reivindicaciones
El artículo posiciona este trabajo como el llenado de un vacío crítico entre la preservación de la positividad (gestionada por métodos de conformación logarítmica/raíz cuadrada) y la compatibilidad termodinámica (gestionada por métodos disipativos de energía libre). Los autores afirman que:

1. La positividad es una condición necesaria pero insuficiente para la robustez ante altos números de Weissenberg.
2. La corrección propuesta es una capa local, computable y controlada espectralmente que puede añadirse a los esquemas de alto orden existentes sin alterar su estructura fundamental.
3. El método evita la inyección no física de energía al asegurar que el tensor que entra en la ecuación de cantidad de movimiento sea compatible con las variables de entropía.
4. Ofrece una ventaja de resolución condicional en regímenes de alto estiramiento al operar en el espacio de variables logarítmicas suaves mientras respeta el presupuesto de entropía.

Los autores concluyen que, para computaciones de alto número de Weissenberg, monitorear el autovalor mínimo no es suficiente; también se debe rastrear el exceso de entropía de la reconstrucción y los defectos de tensión-fuerza para asegurar la fidelidad física. La reconstrucción logarítmica corregida proporciona un mecanismo para imponer esta compatibilidad permaneciendo asintóticamente inactiva para flujos suficientemente resueltos.