A Kinetic Phase-Field Model of Diffusion Bonding: A Nonlocal Approach to Interface Coalescence

Resumo Técnico: Um Modelo de Campo de Fase Cinético para Soldagem por Difusão: Uma Abordagem Não Local para Coalescência de Interfaces

1. O Problema

A soldagem por difusão no estado sólido é uma tecnologia crítica para a fabricação avançada (nuclear, aeroespacial, armazenamento de energia), onde a microestrutura da interface determina o desempenho e a segurança. No entanto, a modelagem computacional desse processo enfrenta um desafio fundamental:

• Limitação dos Modelos Convencionais: Os modelos de campo de fase tradicionais (baseados na equação de Allen-Cahn e na minimização de energia livre) tendem a promover a coalescência (fusão) de interfaces sólido-sólido sempre que elas estão próximas.
• Consequência Indesejada: Em processos de soldagem por difusão, é frequentemente desejável que certas camadas intermediárias (interlayers) ou interfaces permaneçam distintas e não desapareçam, mesmo sob condições termodinâmicas que, em modelos padrão, levariam à sua homogeneização.
• Falha Atual: As abordagens existentes dependem de ajustes numéricos delicados ou de paisagens de energia complexas para tentar impedir a coalescência, o que muitas vezes falha em capturar a física real de interfaces que devem permanecer estáveis sob condições cinéticas ou termodinâmicas específicas.

2. Metodologia

Os autores desenvolveram um modelo de campo de fase cinético inovador que substitui a evolução baseada em gradiente de energia (típica de Allen-Cahn) por uma equação de evolução baseada em uma lei de conservação geométrica para interfaces.

• Equação de Evolução: O modelo utiliza uma formulação onde a velocidade da interface é controlada por uma função cinética que depende da geometria local do campo de fase ($\phi$), e não apenas da força motriz termodinâmica.
  • A equação de evolução é dada por: $\frac{\partial \phi}{\partial t} = -\kappa [1 - g(\phi, \nabla \phi, \dots)] |\nabla \phi| \delta_\phi \mathcal{P}$.
  • O termo $g$ é uma função de coalescência que atua como um fator de frenagem.
• Critério Geométrico Não Local ($g$): A função $g$ é projetada para detectar configurações geométricas específicas que indicam uma camada fina persistente (como uma camada intermediária de titânio entre cerâmicas).
  • Detecção de Mínimos Locais: A função $g$ ativa-se (valor $\approx 1$) quando o campo de fase exibe um mínimo espacial local (indicando uma camada fina) com curvaturas positivas (convexidade).
  • Condições de Ativação: A ativação depende de limiares para o valor de $\phi$, o gradiente $|\nabla \phi|$ e as derivadas de segunda ordem (curvaturas principais ou invariantes do Hessian).
  • Mecanismo de Controle: Quando $g \approx 1$, a mobilidade efetiva da interface $(1-g)$ torna-se zero, parando a evolução da interface e preservando a camada, mesmo que a energia livre ainda "deseje" a coalescência.
• Implementação Numérica:
  • Utiliza-se uma formulação de Elementos Finitos Mistas (Mixed FEM) para lidar com as derivadas de segunda ordem necessárias para calcular as curvaturas.
  • Um campo vetorial auxiliar $\psi = \nabla \phi$ é introduzido para garantir a consistência na avaliação dos termos geométricos.
  • O modelo é implementado na plataforma de código aberto FEniCS.
• Consistência Termodinâmica: O modelo foi provado satisfazer a Segunda Lei da Termodinâmica, garantindo que a taxa de dissipação de energia seja não-negativa, pois a função de frenagem $g$ está limitada ao intervalo $[0, 1]$.

3. Contribuições Principais

1. Mecanismo de Frenagem Geométrica: Introdução de um critério geométrico não local explícito na lei cinética para controlar a coalescência, permitindo que interfaces permaneçam distintas sob condições onde modelos convencionais as eliminariam.
2. Flexibilidade Física: A capacidade de acoplar os limiares de ativação da função $g$ a variáveis macroscópicas (como temperatura ou tensão), permitindo modelar a transição entre "soldagem completa" e "preservação de interface" baseada em condições experimentais.
3. Redução de Complexidade Computacional: Diferente de modelos que acoplam explicitamente equações de transporte de soluto (como Cahn-Hilliard) para controlar a espessura da camada, esta abordagem captura o resultado macroscópico sem resolver equações de difusão adicionais, reduzindo o custo computacional.
4. Validação Multidimensional: Demonstração robusta da detecção de interfaces em 1D, 2D e 3D, incluindo casos complexos como inclusões circulares e interfaces diagonais.

4. Resultados

Os autores realizaram simulações numéricas para validar o modelo:

• Simulações 1D e 2D:
  • Em cenários convencionais ($g=0$), as interfaces coalescem e a camada intermediária desaparece.
  • Com a função $g$ ativa, as interfaces aproximam-se até que os critérios geométricos (mínimo local, curvatura positiva) sejam atendidos, momento em que a evolução é arrestada, preservando uma camada estável de espessura controlada.
  • O modelo demonstrou robustez frente a variações de malha e passo de tempo, e a sensibilidade aos parâmetros de limiar ($\alpha, \beta, \gamma$) permite calibrar a espessura final da camada preservada.
• Aplicação ao Sistema Ti-ZrC (Soldagem por Difusão):
  • O modelo foi aplicado ao sistema de soldagem de compósitos ZrC-SiC com uma camada interlayer de Titânio (Ti).
  • Espessura Crítica: Simulações mostraram que, para uma camada inicial de 200 µm, o modelo prevê a estabilização de uma camada residual de Ti, enquanto para camadas muito finas (50 µm), a coalescência pode ser controlada para atingir uma espessura específica desejada.
  • Acoplamento Mecânico: A inclusão de tensões elásticas e deformações eigenstrain (devido à transformação de fase $\beta$-Ti para $\gamma$-ZrC) mostrou que, embora as tensões alterem a forma local da interface, a cinética de frenagem geométrica mantém a topologia global da interface preservada.

5. Significância e Impacto

Este trabalho oferece uma ferramenta computacional poderosa para a engenharia de materiais avançados:

• Precisão no Projeto de Processos: Permite prever e controlar a microestrutura final de juntas soldadas por difusão, essencial para aplicações onde a presença de uma camada intermediária é crítica para a resistência mecânica ou estabilidade térmica.
• Superação de Limitações Teóricas: Resolve a contradição entre a tendência natural de minimização de energia (que funde interfaces) e a necessidade física de manter interfaces distintas em certas condições cinéticas.
• Versatilidade: A abordagem não se limita à soldagem por difusão; pode ser aplicada a outros problemas de transformação de fase onde a geometria da interface dita o comportamento, como o "pinning" de contornos de grão ou a preservação de trincas sob certas condições.

Em suma, o modelo proposto preenche uma lacuna crítica na modelagem de materiais, oferecendo um mecanismo físico e matematicamente consistente para controlar a coalescência de interfaces, alinhando as simulações de campo de fase com as realidades observadas experimentalmente em processos de união de materiais avançados.