A hybrid sharp-diffuse interface approach to accurately model melt pool dynamics with rapid evaporation in laser-based processing of metals

Resumo Técnico: Uma Abordagem de Interface Híbrida Aguda–Difusa para a Dinâmica de Poças de Fusão

Definição do Problema
A simulação preditiva da dinâmica de poças de fusão em processos de processamento de metais baseados em laser (ex: soldagem por feixe de laser, fusão em leito de pó) requer a resolução de interações termohidrodinâmicas complexas na interface metal–gás. Forças motrizes fundamentais, como a pressão de recuo induzida pela evaporação e a tensão superficial dependente da temperatura, dependem sensivelmente — frequentemente de forma exponencial — da temperatura da interface. As formulações padrão de interface difusa (DI), embora robustas no tratamento de mudanças complexas de topologia como ruptura e coalescência, suavizam gradientes térmicos extremos através de uma região de transição finita. Esse suavizamento introduz erros críticos na previsão da temperatura da interface, o que, consequentemente, degrada a precisão das previsões das forças de interface e da resultante dinâmica da poça de fusão. Por outro lado, os métodos de interface aguda (SI) oferecem alta fidelidade térmica, mas frequentemente enfrentam a complexidade numérica e a robustez necessária para fluxos multifásicos altamente dinâmicos envolvendo grandes razões de densidade e mudanças topológicas. Os modelos existentes tipicamente impõem um compromisso entre fidelidade física e viabilidade computacional, exigindo resoluções de malha proibitivamente finas ou calibrações extensivas de parâmetros para compensar as imprecisões térmicas.

Metodologia
Os autores propõem uma abordagem de Interface Híbrida Aguda–Difusa (HSDI) que acopla um modelo térmico de interface aguda com um modelo de fluxo de interface difusa para alcançar uma modelagem de alta fidelidade sem incorrer nos custos computacionais de uma simulação multifásica totalmente de interface aguda.

1. Modelo Térmico de Interface Aguda (CutFEM): O problema de transferência de calor é resolvido utilizando um método de elementos finitos não ajustado de interface aguda (CutFEM). Esta formulação resolve com precisão as descontinuidades de temperatura e fluxo de calor na interface metal–gás. Os autores apresentam duas variantes:

   • SI-2P (Acoplado de Duas Fases): Resolve a equação do calor tanto no metal quanto na fase de gás ambiente, impondo fracamente a continuidade de temperatura através da interface usando o método de Nitsche.
   • SI-1P (Fase Única): Resolve a equação do calor apenas no domínio metálico, tratando o gás ambiente como adiabático, exceto para fluxos de calor impostos. Esta variante é motivada pela capacidade térmica negligenciável do gás ambiente em modelos de evaporação não resolvidos, oferecendo maior eficiência computacional.
     Ambas as variantes empregam estabilização de penalidade fantasma baseada em faces para lidar com elementos cortados pequenos.

2. Modelo de Fluxo de Interface Difusa: O problema de fluxo multifásico, caracterizado por grandes razões de densidade e dinâmicas de interface complexas, é modelado utilizando uma formulação de elementos finitos de um único fluido de interface difusa baseada em nível de conjunto (level-set). Esta abordagem utiliza um modelo de Força de Superfície Contínua (CSF) para tensão superficial e pressão de recuo, permitindo o tratamento natural de ruptura e coalescência de interface.

3. Estratégia de Acoplamento Consistente: Para unir as duas estruturas, o campo de temperatura de interface aguda é estendido para uma banda estreita ao redor da interface metal–gás usando um algoritmo de projeção de ponto mais próximo. Este campo de temperatura estendido ($\tilde{T}_\Gamma$) é então utilizado para avaliar as forças de interface dependentes da temperatura dentro do framework de fluxo de interface difusa. Isso garante que os dados térmicos altamente precisos do modelo SI conduzam o solver de fluxo DI de forma consistente.

4. Estrutura Numérica: A implementação utiliza a biblioteca deal.II com avaliação de operador livre de matriz (matrix-free), refinamento de malha adaptativo e paralelização baseada em MPI. Um esquema de divisão de operador particionado é empregado para integração temporal, tratando termos de acoplamento explicitamente enquanto resolve campos individuais implicitamente.

Principais Contribuições

• Nova Formulação Híbrida: Este trabalho apresenta a primeira abordagem HSDI para termohidrodinâmica de duas fases que acopla a modelagem de temperatura CutFEM de interface aguda com a modelagem de fluxo FEM de interface difusa baseada em um modelo de evaporação não resolvido.
• Precisão Térmica: O desenvolvimento de formulações SI baseadas em CutFEM (SI-2P e SI-1P) que permitem prever com alta precisão o campo de temperatura, particularmente na superfície da poça de fusão, onde os modelos DI padrão falham.
• Algoritmo de Acoplamento: Uma estratégia consistente para estender a temperatura de interface aguda para a região de transição difusa, permitindo a avaliação precisa de termos de fonte dependentes da temperatura no solver de fluxo DI.
• Benchmarking: A proposta de um novo exemplo de benchmark 2D representando uma depressão de vapor induzida por laser para avaliar de forma abrangente a precisão dos modelos termohidrodinâmicos acoplados.

Resultos

• Convergência Térmica: Em benchmarks puramente térmicos mimetizando interações laser–metal, os modelos SI exibiram convergência espacial de segunda ordem ($O(h^2)$), enquanto os modelos DI padrão mostraram apenas convergência de primeira ordem ($O(h^1)$).
• Precisão vs. Resolução de Malha: Para uma precisão alvo de 1% na temperatura da interface e na pressão de recuo, os modelos SI permitiram tamanhos de elementos finitos duas ordens de magnitude maiores do que os exigidos pelos modelos DI. Especificamente, a abordagem SI atingiu a tolerância do benchmark com uma malha aproximadamente 50 vezes mais grosseira que a abordagem DI.
• Termohidrodinâmica Acoplada: No novo benchmark de depressão de vapor 2D, a abordagem HSDI mostrou-se um ordem de magnitude mais precisa do que um modelo puramente DI na mesma malha. A abordagem HSDI previu com sucesso as dinâmicas estáveis e instáveis (keyhole) da depressão de vapor com erro significativamente menor na geometria da interface e energia cinética em comparação ao modelo DI.
• Robustez 3D: Uma simulação de prova de conceito 3D de fusão estacionária induzida por laser de uma placa de metal nua demonstrou a robustez da abordagem HSDI em lidar com fluxos de duas fases complexos e altamente dinâmicos, incluindo formação de keyhole, ruptura e coalescência, sem instabilidade numérica.
• Eficiência do SI-1P: A variante térmica de fase única (SI-1P) mostrou preservar uma precisão comparável à variante SI-2P totalmente acoplada, ao mesmo tempo em que oferece eficiência computacional, validando a suposição de que a transferência de calor na fase gasosa é negligenciável neste contexto.

Significância
O artigo afirma que a abordagem HSDI proporciona um equilíbrio sem precedentes entre precisão preditiva e eficiência computacional para simulações de poças de fusão. Ao desacoplar os problemas térmicos e de fluxo em suas respectivas formulações ideais (aguda para térmica, difusa para fluxo), o método supera as limitações dos modelos puramente DI, que exigem malhas impaticavelmente finas para resolver gradientes térmicos, e dos modelos puramente SI, que têm dificuldade em lidar com mudanças topológicas complexas. A metodologia permite previsões de alta fidelidade da dinâmica da poça de fusão em resoluções de malha viáveis para escalas de tempo e comprimento industrialmente relevantes. Os autores enfatizam que esta abordagem é aplicável a uma ampla gama de problemas de fluxo dependentes de temperatura envolvendo fontes de calor fortemente localizadas, como soldagem por feixe de laser e fusão em leito de pó por laser, e serve como base para trabalhos futuros que incorporem físicas mais detalhadas, como ray-tracing e absortividade dependente da temperatura.