A Temperature-Coupled Cahn-Hilliard-Stokes-Heat Model for Thermally Driven Phase Separation

1. Problema e Motivação

O trabalho aborda a modelagem matemática e numérica da separação de fases em misturas viscosas incompressíveis, especificamente focado na condensação de ar úmido em geometrias confinadas, como as caixas de baterias de veículos elétricos.

• Contexto: A condensação de umidade em superfícies internas mais frias pode causar corrosão e degradação de isolamento. Testes experimentais abrangentes são caros, tornando a simulação numérica uma ferramenta essencial.
• Desafio: A maioria dos modelos existentes foca em cenários isotérmicos ou em configurações genéricas de ebulição/condensação, sem considerar adequadamente o acoplamento com escoamentos de Stokes de densidade variável e condições de contorno realistas em geometrias complexas.
• Objetivo: Desenvolver e analisar um modelo de fase difusa (diffuse-interface) que acople a equação de Cahn–Hilliard, as equações de Stokes e a equação de calor, permitindo a simulação de transições de fase induzidas termicamente (como a condensação) controladas por gradientes de temperatura.

2. Metodologia e Formulação Matemática

O modelo proposto é um sistema acoplado não linear que integra três componentes físicas principais:

A. Potencial Termodinâmico Dependente da Temperatura

Os autores modificam o potencial clássico de "duplo poço" (double-well) para torná-lo dependente da temperatura.

• Potencial de Landau: A densidade de energia livre de bulk é dada por $f(c, \Theta) = 2c^4 - (1 - \beta(\Theta))c^2$, onde $\beta(\Theta) = \Theta/\Theta_S$ é a razão entre a temperatura local e a temperatura de saturação.
• Mecanismo:
  • Se $\Theta > \Theta_S$ ($\beta > 1$), o potencial tem um único mínimo em $c=0$, favorecendo uma fase homogênea (vapor estável).
  • Se $\Theta < \Theta_S$ ($\beta < 1$), o potencial assume a forma de duplo poço, favorecendo a separação de fases (condensação).
• Parâmetros: O modelo não inclui umidade explícita ou calor latente como fontes separadas; esses efeitos são absorvidos na escolha da temperatura de saturação $\Theta_S$.

B. Sistema Acoplado (Cahn–Hilliard–Stokes–Calor)

O sistema governante em variáveis adimensionais consiste em:

1. Equação de Cahn–Hilliard Convectiva: Descreve a evolução do parâmetro de ordem $c$ (fração volumétrica), acoplado ao campo de velocidade $u$ e ao potencial químico $\mu$ que depende de $\Theta$.
2. Equações de Stokes: Descrevem o escoamento do fluido com densidade $\rho(c)$ e viscosidade $\eta(c)$ variáveis (leis afins). Inclui um termo de empuxo gravitacional ($G \rho(c) \hat{g}$) para modelar a sedimentação. O número de Reynolds é assumido muito baixo ($Re \ll 1$).
3. Equação de Calor: Descreve a evolução da temperatura $\Theta$ por convecção e difusão térmica.

C. Análise Matemática (Existência de Soluções)

Os autores provam a existência local no tempo de soluções fracas para o sistema regularizado em 2D e 3D.

• Abordagem: Utilização do Teorema do Ponto Fixo de Schauder.
• Estratégia: O sistema é decomposto em operadores:
  • Operador de Calor ($H$): Mapeia velocidade para temperatura.
  • Operador de Cahn–Hilliard ($C$): Mapeia (velocidade, temperatura) para o parâmetro de ordem.
  • Operador de Stokes ($S$): Mapeia concentração para velocidade.
• Regularização: Para fins analíticos, o sistema é regularizado com um termo de derivada temporal no Stokes ($-\alpha \partial_t u$) e uma suavização espacial (mollification) da temperatura no potencial químico.
• Desafios Analíticos: A ausência de limites globais $L^\infty$ para o parâmetro de ordem (devido ao potencial polinomial) e as condições de contorno mistas exigem estimativas a priori em topologias mais fracas ($H^2 \to L^2 \to H^{-2}$).

D. Esquema Numérico

Foi desenvolvido um esquema semi-implícito totalmente discreto baseado em Elementos Finitos:

• Discretização Espacial: Elementos de Taylor-Hood ($P2-P1$) para Stokes e elementos $P1$ contínuos para $c$, $\mu$ e $\Theta$.
• Discretização Temporal:
  • Equação de Calor: Euler para trás (implícito).
  • Stokes: Resolver estacionário (quase-estático) com coeficientes atualizados.
  • Cahn–Hilliard: Convex Splitting (divisão convexa-concava) estabilizada. A parte convexa da energia é tratada implicitamente e a parte côncava explicitamente, garantindo estabilidade energética incondicional para o subproblema difusivo isotérmico.
• Software: Implementado em FreeFEM++.

3. Principais Contribuições

1. Modelo Termodinamicamente Motivado: Introdução de um potencial de Landau onde a transição de fase é controlada diretamente pelo parâmetro termodinâmico $\beta(\Theta)$, permitindo simular a transição entre mistura homogênea e separação de fases sem necessidade de fontes de calor latente ad hoc.
2. Acoplamento com Escoamento de Stokes Variável: Integração rigorosa de leis de densidade e viscosidade dependentes da concentração e do termo de empuxo, essencial para modelar a sedimentação de gotículas em baterias.
3. Teoria de Existência Sob Condições de Contorno Mistas: Prova de existência de soluções fracas para um sistema complexo com condições de contorno de deslizamento (no-slip) e tensão livre (traction-free), além de condições térmicas mistas, em um contexto não isotérmico.
4. Esquema Numérico Robusto: Desenvolvimento de um algoritmo de splitting que preserva a conservação de massa e a estabilidade energética (no caso isotérmico), validado contra benchmarks clássicos.

4. Resultados Numéricos

Os autores realizaram quatro experimentos numéricos em 2D para validar o modelo:

1. Benchmark de Cavidade Acionada por Tampo (Lid-Driven Cavity): Comparação com resultados da literatura para validar o acoplamento hidrodinâmico e o transporte da interface. O modelo reproduziu com precisão os perfis de velocidade e a deformação da interface.
2. Separação de Fases Isotérmica:
   • $\Theta > \Theta_S$: O sistema relaxa para um estado homogêneo (potencial de poço único).
   • $\Theta < \Theta_S$: Ocorre separação de fases rápida seguida de coarsening (amadurecimento), validando o comportamento do potencial de duplo poço.
3. Separação de Fases Induzida por Gravidade: Simulação com densidades e viscosidades diferentes. Mostrou a migração de domínios densos para a parte inferior do domínio devido ao empuxo, formando uma camada pesada na base, demonstrando o acoplamento correto entre fase e gravidade.
4. Mudança de Fase Térmica em Fronteiras:
   • Fronteira Aquecida: Um gradiente de temperatura que aquece a parede localmente suprime a separação de fases, fazendo com que as estruturas existentes se dissolvam e o sistema se torne homogêneo perto da parede.
   • Fronteira Resfriada: O resfriamento de uma parede induz a nucleação de fases separadas (condensação) na vizinhança da parede fria, que crescem ao longo do tempo.

5. Significado e Conclusão

O trabalho fornece uma ferramenta matemática e computacional robusta para estudar a separação de fases induzida termicamente em geometrias confinadas, com aplicações diretas na engenharia de baterias e sistemas térmicos.

• O modelo captura qualitativamente fenômenos críticos como nucleação de gotas, coalescência, sedimentação e mudança de fase induzida por paredes.
• A análise de sensibilidade aos parâmetros adimensionais (números de Péclet, número de Cahn, razões de densidade/viscosidade e gravidade) oferece insights sobre como esses fatores influenciam a formação de padrões e a estrutura do escoamento.
• A abordagem proposta preenche uma lacuna entre modelos termodinâmicos ideais e simulações de engenharia realistas, permitindo a investigação de cenários de condensação em condições operacionais complexas sem a necessidade de testes experimentais extensivos.