Unstable Interface Dynamics for Gravity Stokes Flow

Resumo Técnico: Dinâmica de Interface Instável para o Fluxo de Stokes por Gravidade

Definição do Problema
Este artigo investiga a dinâmica de longo prazo e os mecanismos de instabilidade do problema da interface nítida para o fluxo de Stokes bidimensional sob ação da gravidade. O sistema modela dois fluidos incompressíveis e homogêneos com densidades constantes distintas ($\rho_+$ e $\rho_-$) evoluindo sob a lei de Stokes com uma força gravitacional. A interface $\Gamma(t)$ separa a região superior $\Omega_+(t)$ da região inferior $\Omega_-(t)$ em um domínio periodicamente horizontal ($\mathbb{T} \times \mathbb{R}$).

O estudo foca em dois regimes primários:

1. Regime Instável (Rayleigh-Taylor): O fluido mais denso está acima do fluido mais leve ($\rho_+ < \rho_-$). Embora a existência global para dados iniciais suaves seja conhecida (ex: [12]), o comportamento de longo prazo de quantidades geométricas (perímetro, curvatura) permanece amplamente inexplorado.
2. Regime Estável: O fluido mais denso está abaixo do fluido mais leve ($\rho_+ > \rho_-$). O artigo examina se configurações iniciais estáveis podem evoluir para estados instáveis, investigando especificamente "instabilidades de tombamento" (turning instabilities), onde a interface deixa de ser um gráfico em tempo finito.

O sistema é analisado via a Equação de Dinâmica de Contorno (CDE) derivada em [11]:
$$\partial_t z(\alpha, t) = (\rho_- - \rho_+) \int_{\mathbb{T}} S(z(\alpha, t) - z(\beta, t)) \dot{z}^\perp(\beta, t) z^2(\beta, t) d\beta$$
onde $S$ é o Stokeslet $\text{x}_1$-periódico.

Metodologia
Os autores empregam uma combinação de técnicas analíticas e simulações numéricas:

• Análise de Energia Potencial: Uma ferramenta central é um funcional de energia potencial modulada $E(t) = \int_{\mathbb{T}\times\mathbb{R}} x_2(\rho_s(x_2) - \rho(x,t)) dx$, onde $\rho_s$ representa o estado estratificado instável. Os autores provam que $E(t)$ é monotonicamente crescente e sua derivada temporal relaciona-se com a norma $L^2$ de um operador não local aplicado à densidade.
• Exploração de Simetrias: A análise baseia-se fortemente em simetrias específicas dos dados iniciais: simetria central ($z(\alpha) = -z(-\alpha)$) e simetria par em relação a $x_1 = \pm \pi/2$. Essas simetrias são provadas como conservadas pelo fluxo e permitem a redução do domínio para uma faixa restrita para estimativas geométricas.
• Dicotomia Geométrica: Para provar resultados de crescimento, os autores estabelecem uma dicotomia baseada na quantidade $\delta(t) = \|( -\Delta)^{-1} \partial_1 \rho(t) \|_{L^2}^2$. Eles demonstram que, se o perímetro $L(t)$ não crescer suficientemente rápido, a quantidade $\delta(t)$ deve tornar-se pequena, o que, por sua vez, força a curvatura máxima $K(t)$ a crescer rapidamente.
• Analiticidade e Fingering (Formação de Dedos): Para a análise de fingering, os autores utilizam a analiticidade das soluções em uma faixa complexa. Eles aplicam um lema sobre a distribuição de pontos críticos de funções analíticas para estimar o número de "dedos" (extremos locais) que podem se desenvolver.
• Construção da Instabilidade de Tombamento: Para provar a quebra em tempo finito no regime estável, os autores constroem dados iniciais específicos que satisfazem condições de simetria ímpar/par. Eles analisam a derivada temporal da inclinação na origem para mostrar que uma interface estável pode desenvolver uma tangente vertical e tombar.
• Simulações Numéricas: Os autores implementam um esquema numérico de tipo espectral usando diferenças centrais e a regra de Simpson para as integrais singulares, regularizadas por viscosidade artificial, para visualizar o crescimento do comprimento da interface e da curvatura.

Contribuições Principais e Resultados

1. Crescimento em Tempo Infinito no Regime Instável (Teorema 1.1):
   O artigo prova que, para dados iniciais instáveis satisfazendo condições de simetria específicas, a interface exibe crescimento ilimitado em pelo menos uma quantidade geométrica (perímetro $L(t)$ ou curvatura máxima $K(t)$) quando $t \to \infty$. Especificamente, para qualquer $\epsilon > 0$:
   $$\limsup_{t \to \infty} t^{-\frac{1}{17} + \epsilon} (L(t) + K(t)) = \infty$$
   Este resultado estabelece que as soluções não se estabilizam para um estado estacionário, mas exibem crescimento de longo prazo genérico, contrastando com o resultado de existência global que se aplica independentemente da estratificação de densidade.

2. Estimativa sobre o Fenômeno de Fingering (Teorema 1.3):
   Os autores fornecem um limite inferior para o crescimento do número de dedos (pontos críticos da inclinação da interface) para uma classe de soluções analíticas. Eles mostram que o número de dedos $F(t)$ cresce pelo menos como $C(1 + \sqrt{t})$. Isso conecta o crescimento analítico da solução ao fenômeno físico de viscous fingering.

3. Instabilidade de Tombamento em Tempo Finito no Regime Estável (Teorema 1.4):
   Ao contrário da intuição de que a estratificação estável previne a instabilidade, os autores provam a existência de dados iniciais suaves no regime estável ($\rho_+ > \rho_-$) que evoluem para uma instabilidade de tombamento em tempo finito. A interface deixa de ser um gráfico e a condição de Rayleigh-Taylor falha. O artigo constrói uma família de tais curvas iniciais e prova que o tombamento ocorre em pontos simétricos específicos.

4. Validação Numérica:
   Simulações numéricas confirmam as previsões teóricas, mostrando crescimento significativo no comprimento da interface e curvatura limitada no curto prazo para configurações instáveis, bem como o desenvolvimento de dedos na evolução da interface.

Significância e Alegações
O artigo afirma abordar uma lacuna na compreensão da dinâmica de longo prazo do sistema Stokes-Transporte. Embora a existência global tenha sido estabelecida em trabalhos anteriores [12], o comportamento dessas soluções globais era desconhecido. Os autores argumentam que seus resultados demonstram que:

• A existência global não implica estabilidade; configurações instáveis levam ao crescimento geométrico ilimitado.
• Configurações estáveis não estão imunes à instabilidade; dados iniciais suaves específicos podem evoluir para singularidades de tombamento, quebrando a propriedade de gráfico.
• O crescimento de quantidades geométricas (perímetro e curvatura) é uma característica genérica do regime instável, apoiando a conexão com os fenômenos de viscous fingering observados em fluxos de meios porosos.

Os autores explicitamente observam que a taxa de crescimento derivada no Teorema 1.1 ($t^{1/17}$) provavelmente não é aguda (sharp), com evidência numérica sugerindo uma taxa de crescimento linear, mas a prova analítica fornece um limite inferior rigoroso para a instabilidade. O trabalho destaca a distinção da dinâmica do fluxo de Stokes em relação ao problema de Muskat, que é mais singular, particularmente quanto à natureza do amortecimento linear e aos mecanismos específicos de instabilidade.