Revisiting creeping viscoelastic cross-slot flow: Global linear stability and structural sensitivity analyses

Resumo Técnico: Revisitando o Fluxo de Cross-Slot Viscoelástico em Regime de Creeping

Definição do Problema
O artigo aborda a origem física não resolvida há muito tempo da instabilidade elástica no fluxo de cross-slot bidimensional em regime de creeping para fluidos Oldroyd-B. Observada experimentalmente há quase meio século e reproduzida numericamente há duas décadas, este fluxo exibe uma transição de quebra de simetria de um estado simétrico estacionário para um estado assimétrico estacionário conforme o número de Weissenberg ($Wi$) excede um limiar crítico. Apesar dos extensos estudos fenomenológicos, o mecanismo preciso que impulsiona esta instabilidade permanece em debate. Hipóteses anteriores atribuíram o início do fenômeno a diversos fatores, incluindo o fluxo compressivo a montante do ponto de estagnação, o fluxo extensional a jusante, singularidades geométricas em cantos agudos, alto cisalhamento próximo aos cantos ou a interação entre elasticidade e pseudoplasticidade (shear-thinning). Os autores observam que carecia uma previsão teórica de primeira ordem sobre o início da instabilidade e a estrutura do modo, em parte devido às dificuldades numéricas associadas à simulação de fluxos viscoelásticos (o "problema do alto número de Weissenberg") em geometrias complexas que carecem de estados de base analíticos.

Metodologia
Os autores empregam uma abordagem abrangente e multifacetada combinando quatro técnicas analíticas e numéricas para um fluido Oldroyd-B difusivo em um canal de cross-slot 2D padrão com número de Reynolds nulo:

1. Simulações Numéricas Diretas (DNS): Utilizando o método de elementos finitos (FEM) via FreeFem++, os autores resolvem as equações de Stokes incompressíveis acopladas à lei constitutiva difusiva de Oldroyd-B. Eles utilizam a formulação de gradiente de tensão dividida elástico-viscosa discreta (DEVSS-G) e o esquema streamline upwind Petrov–Galerkin (SUPG) para estabilização numérica.
2. Análise de Estabilidade Linear Global (LSA): Os autores linearizam as equações governantes em torno do estado de base estacionário para resolver os autovalores e automodos complexos. Isso permite a previsão do $Wi$ crítico e das taxas de crescimento temporal das perturbações.
3. Análise de Sensibilidade Estrutural: Ao computar tanto os automodos diretos quanto os adjuntos, os autores definem campos de sensibilidade estrutural (baseados tanto em velocidade quanto em tensão) para identificar a região do "gerador de ondas" (wavemaker) — o núcleo espacial onde a instabilidade é mais sensosa a um feedback localizado.
4. Análise de Orçamento de Energia: Os autores decompõem a equação de evolução da energia da perturbação para identificar os mecanismos físicos específicos que contribuem para a produção ou dissipação de energia. Eles focam no termo de acoplamento entre o estado de base e as perturbações, analisando especificamente a densidade de potência polimérica para localizar a fonte energética da instabilidade.

O estudo também investiga duas variantes geométricas do cross-slot base: uma com cantos arredondados (removendo singularidades geométricas) e outra com um cilindro centralizado (eliminando o ponto de estagnação livre), para testar a robustez dos mecanismos identificados.

Resultos Principais

• Bifurcação e Validação: A DNS reproduz com sucesso a bifurcação de garfo (pitchfork) supercrítica canônica, com o fluxo transitando de simétrico para assimétrico em um número de Weissenberg crítico $Wi_{cr} \approx 0,36$. A LSA prevê consistentemente este limiar e as taxas de crescimento temporal, que se alinham com os dados de DNS e os resultados de Decomposição de Modos Dinâmicos (DMD).
• Estrutura do Modo de Instabilidade: O automodo direto principal consiste em uma perturbação quiral de velocidade-tensão. A perturbação de velocidade forma uma estrutura vortical em formato de diamante centrada no ponto de estagnação. Crucialmente, a perturbação da tensão normal polimérica ($\hat{\tau}_{yy}$) apresenta paridade antissimétrica, o que inclina e rotaciona a fibra birrefringente gerada pelo fluxo extensional. Esta rotação desvia as separatrizes do fluxo, levando a uma distribuição desigual de fluxo nos ramos de saída.
• Núcleo Espacial (Wavemaker): A análise de sensibilidade estrutural revela que a instabilidade não é impulsionada pelos cantos ou pelo ponto de estagnação em si. Em vez disso, o wavemaker baseado em tensão está localizado em cristas estreitas de alta taxa de extensão dentro do fluxo extensional a jusante do ponto de estagnação. Estas cristas exibem um padrão de "espinha de peixe" (herringbone) de perturbação de conformação polimérica, caracterizado por estiramento e compressão lateral simultâneos.
• Fonte Energética: A análise do orçamento de energia identifica o acoplamento entre o estado de base e a perturbação como a fonte de energia dominante. Especificamente, a densidade de potência polimérica de perturbação ($P$) é positiva e máxima ao longo das mesmas cristas verticais identificadas pela análise do wavemaker. Isso indica uma transferência localizada de energia elástica armazenada das cadeias poliméricas deformadas de volta para a energia cinética da perturbação, impulsionando a instabilidade.
• Variantes Geométricas: Em ambas as variantes — cantos arredondados e cilindro centralizado — a instabilidade persiste. O pico de densidade de potência polimérica permanece localizado ao longo de cristas verticais nas regiões de fluxo extensional. Isso demonstra que nem os cantos agudos, nem o ponto de estagnação central livre são o ingrediente desestabilizador essencial.

Significância e Alegações
O artigo afirma resolver a controvérsia sobre a origem física da instabilidade de cross-slot elástica. Os autores concluem que a instabilidade não decorre de singularidades geométricas (cantos) ou do ponto de estagnação propriamente dito, mas sim da liberação de energia elástica em regiões dominadas por extensão, características do fluxo de cross-slot.

Especificamente, o mecanismo envolve:

1. A geração de uma fibra birrefringente com tensões normais e gradientes substanciais na extensão do fluxo a jusante do ponto de estagnação.
2. Uma perturbação quiral que inclina e rotaciona esta fibra.
3. A liberação de energia elástica armazenada em cristas estreitas de alta taxa de extensão (a região do wavemaker), que alimenta o crescimento da instabilidade.

Os autores afirmam que suas descobertas reconciliam interpretações conflitantes anteriores: enquanto os cantos agudos não são necessários (como mostrado pelos experimentos OSCER), e o ponto de estagnação não é o único motor (como mostrado pelos experimentos de inserção de cilindro), o ingrediente essencial é o fluxo de alta taxa de extensão e o consequente estiramento polimérico. O estudo valida ainda a utilidade da integração de estabilidade global, sensibilidade estrutural e análise de orçamento de energia para compreender instabilidades em fluxos viscoelásticos complexos onde estados de base analíticos não estão disponíveis.