Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

Resumo Técnico: Reconstrução Compatível com a Entropia para Fluxo Viscoelástico de Alto Número de Weissenberg

Definição do Problema
O problema do alto número de Weissenberg permanece como um obstáculo central na simulação de fluxos viscoelásticos. Em modelos de tensor de conformação, o estado polimérico é representado por um tensor simétrico definido positivo (SPD) $A$. Embora a positividade seja essencial para o significado físico da energia livre elástica ($\log \det A$) e do logaritmo de matriz, o artigo argumenta que a positividade por si só é insuficiente para a estabilidade numérica em números de Weissenberg elevados.

Métodos existentes, como a conformação logarítmica ($A = \exp \Psi$) e reconstruções de raiz quadrada ($A = BB^T$), conseguem impor a restrição SPD. No entanto, experimentos numéricos mostram que esses métodos ainda podem produzir sobre-sinais (overshoots) de tensão não físicos e crescimento não físico da energia elástica. O problema central identificado é uma desconexão entre a restrição de cone (positividade) e a restrição de balanço acoplado (compatibilidade da energia livre). Especificamente, se o tensor usado no trabalho de tensão da equação do momento diferir do tensor usado no trabalho de estiramento ou nas variáveis de entropia da equação de conformação, surge um defeito de trabalho de sinal indefinido. Além que, reconstruções estritamente positivas podem introduzir três mecanismos não físicos específicos:

1. Viés de entropia do tipo Jensen: Mapeamentos não lineares (como o exponencial) aplicados a perturbações de média zero no espaço logarítmico geram um viés médio positivo no tensor físico e na entropia.
2. Amplificação exponencial: Em estados altamente estirados, pequenas perturbações logarítmicas são amplificadas pelo exponencial de matriz, escalando com o maior autovalor $\lambda_{\max}(A)$.
3. Defeitos de trabalho de sinal indefinido: O uso de tensores reconstruídos incompatíveis para os termos de tensão e entropia impede o cancelamento exato do trabalho polimérico, levando à injeção ou remoção descontrolada de energia.

Metodologia
O artigo propõe uma reconstrução logarítmica corrigida e compatível com a entropia. Em vez de substituir os métodos de preservação de positividade existentes, esta abordagem atua como um filtro de admissibilidade aplicado a uma reconstrução logarítmica de alta ordem bruta antes que ela entre no balanço acoplado de tensão-trabalho.

A metodologia envolve:

• Definição do Caminho Logarítmico: Dado um preditor físico $\hat{A}$ e uma reconstrução de alta ordem bruta $\tilde{A}$, um caminho é definido no espaço logarítmico: $A(\theta) = \exp(\hat{\Psi} + \theta(\tilde{\Psi} - \hat{\Psi}))$ para $\theta \in [0, 1]$.
• Restrição do Orçamento de Entropia: O método busca o maior $\theta$ tal que a energia livre discreta total (entropia) não exceda a entropia do preditor físico mais um orçamento local $\tau_h$. O orçamento é escalonado para corresponder ao erro de truncamento esperado do esquema de alta ordem (ex: $O(h^{2k+2})$).
• Algoritmo de Bisseção: O $\theta$ ótimo é encontrado via uma busca de bisseção local por célula. Se a reconstrução bruta satisfizer o orçamento, $\theta=1$ (sem alteração). Caso contrário, o método recua ao longo do caminho até o ponto admissível máximo.
• Imposição de Consistência: O tensor selecionado é usado de forma consistente em todos os termos de acoplamento: trabalho de tensão do momento, estiramento de conformação, variáveis de entropia e quadratura.

Contribuições Principais e Resultados Teóricos
O artigo fornece uma análise rigorosa deste princípio de reconstrução, provando várias propriedades fundamentais:

• Existência e Convexidade: O funcional de entropia ao longo do caminho logarítmico é convexo, garantindo a existência de um parâmetro admissível máximo $\theta^*$.
• Mínimo Amortecimento: Entre todos os estados admissíveis no caminho, o $\theta^*$ selecionado minimiza a distância em relação à reconstrução bruta, tornando a correção de "mínimo amortecimento".
• Controle Espectral: A correção preserva a independência de referencial (frame indifference) e mantém limites espectrais derivados dos logaritmos dos pontos finais.
• Inatividade Assintótica: Se a reconstrução de alta ordem bruta já for entropicamente admissível (ou seja, seu defeito está dentro da escala do erro de discretização), a correção torna-se inativa ($\theta^* = 1$), preservando a precisão formal de alta ordem.
• Vantagem de Resolução em Alto Estiramento: Em regimes de alto estiramento, o método retém alta precisão em métricas log-relativas e de entropia, enquanto reconstruções diretas no espaço físico sofrem de problemas de escala de derivada proporcionais a $\lambda_{\max}(A)$.
• Estimativa de Energia: Os autores provam uma estimativa de energia livre discreta para o sistema acoplado, mostrando que o orçamento de reconstrução $\tau_h$ atua como uma fonte controlada e limitada de entropia, em vez de viscosidade artificial.

Resultados Numéricos
O artigo valida a teoria através de uma série de diagnósticos:

• Testes Escalares: Demonstrou que perturbações logarítmicas de média zero criam um viés de entropia positivo (viés de Jensen) e que reconstruções lineares, embora de baixo viés, falham em preservar a positividade.
• Testes de Amplificação: Confirmou que defeitos de força de tensão em reconstruções logarítmicas brutas crescem linearmente com o estiramento de fundo $\lambda_{\max}(A)$.
• Verificação de Defeito de Trabalho: Mostrou que o uso de tensores incompatíveis para tensão e entropia leva a defeitos de energia de sinal indefinido.
• Loops de Alto Número de Weissenberg Dinâmicos: Em simulações com $Wi=80$ e alto estiramento, as reconstruções logarítmicas padrão acumularam erros de entropia e de tensão-força ao longo do tempo, enquanto a reconstrução corrigida controlou esses defeitos, reduzindo os erros de força em ordens de magnitude.
• Benchmarks de Resolução: O método corrigido manteve a compatibilidade de entropia e controlou os defeitos de força através de refinamentos de malha e diferentes níveis de estiramento, enquanto os métodos brutos mostraram erros crescentes com o aumento do estiramento.

Significância e Alegações
O artigo posiciona este trabalho como o preenchimento de uma lacuna crítica entre a preservação de positividade (tratada por métodos de conformação logarítmica/raiz quadrada) e a compatibilidade termodinâmica (tratada por métodos dissipativos de energia livre). Os autores afirmam que:

1. A positividade é uma condição necessária, mas insuficiente, para a robustez em altos números de Weissenberg.
2. A correção proposta é uma camada local, computável e controlada espectralmente que pode ser adicionada a esquemas de alta ordem existentes sem alterar sua estrutura fundamental.
3. O método evita a injeção não física de energia ao garantir que o tensor que entra na equação do momento seja compatível com as variáveis de entropia.
4. Ele oferece uma vantagem de resolução condicional em regimes de alto estiramento ao operar no espaço de variáveis logarítmicas suaves, respeitando o orçamento de entropia.

Os autores concluem que, para computações de alto número de Weissenberg, monitorar o autovalor mínimo não é suficiente; deve-se também rastrear o excesso de entropia da reconstrução e os defeitos de tensão-força para garantir a fidelidade física. A reconstrução logarítmica corrigida fornece um mecanismo para impor essa compatibilidade, permanecendo assintoticamente inativa para fluxos suficientemente resolvidos.