Thermodynamically consistent phase field model for hydrogen-assisted cracking

Este artigo apresenta um modelo de campo de fase termodinamicamente consistente que simula o fissuramento assistido por hidrogênio em materiais policristalinos ao acoplar a propagação de trincas com a segregação de hidrogênio e a redução da energia interfacial, capturando com sucesso a transição do fracasso transgranular para o intergranular sob mecanismos de decoesão potencializada por hidrogênio.

O essencial

Imagine que você tenha uma estrutura metálica muito forte, como o chassi de um carro ou uma ponte. Espera-se que ela suporte a pressão, mas às vezes, átomos invisíveis de hidrogênio entram furtivamente no metal e fazem com que ele se estilhace inesperadamente. Esse fenômeno é chamado de fragilização por hidrogênio. É como se o metal estivesse sendo secretamente "envenenado" por dentro, tornando-o quebradiço e propenso a romper.

Cientistas têm tentado construir modelos computacionais para prever exatamente como e onde esse metal irá quebrar. No entanto, modelos anteriores tinham uma falha importante: tratavam o comportamento do hidrogênio como uma regra simples e uniforme que se aplica em todos os lugares, embora a estrutura interna do metal seja, na verdade, uma colcha de retalhos complexa de diferentes grãos e contornos.

O Novo Modelo "Inteligente"
Os autores deste artigo criaram uma simulação computacional mais sofisticada (chamada de "modelo de campo de fase") que atua como um mapa de alta definição e termodinamicamente consistente. Veja como funciona, usando algumas analogias do cotidoso:

• O Metal como uma Multidão: Imagine que o metal é uma sala lotada de pessoas (átomos de metal). Os "contornos de grão" são as linhas invisíveis que separam diferentes grupos de pessoas. A "fenda" é uma lacuna crescente na multidão.
• O Hidrogênio como um Convidado Pegajoso: Átomos de hidrogênio são como convidados pegajosos que adoram se esconder nos espaços vazios entre as pessoas. Eles têm uma preferência especial: adoram se fixar nas bordas da fenda e nas linhas entre os grupos ainda mais do que gostam do meio da multidão.
• O Problema da "Cola": Em um metal saudável, a "cola" que mantém as bordas da fenda unidas é forte. Mas quando esses convidados pegajosos de hidrogênio se reúnem nas bordas da fenda, eles agem como um óleo escorregadio, enfraquecendo a cola. Isso torna a abertura da fenda muito mais fácil.
• A Abordagem Antiga vs. Nova:
  • Modelos Antigos: Usavam um livro de regras genérico (isoterma de Langmuir-McLean) que assumia que o hidrogênio estava distribuído uniformemente e em perfeito equilíbrio em todos os lugares. Isso é como assumir que todos na sala lotada estão parados e distribuídos uniformemente, o que não é verdade quando uma fenda está se formando.
  • Novo Modelo: Usa uma estrutura "variacional" flexível (baseada no formalismo de Kim-Kim-Suzuki). Em vez de forçar uma regra rígida, ele permite que o hidrogênio "migre" naturalmente para onde deseja ir (as bordas da fenda e os contornos de grão) com base nas condições locais. Ele calcula exatamente o quanto a "cola" enfraquece em tempo real conforme o hidrogênio se acumula.

O Que Eles Descobriram
A equipe testou seu novo modelo com dois cenários principais:

1. O Teste da Fenda Única: Eles simularam uma fenda em uma única peça de metal. Sem hidrogênio, a fenda cresceu exatamente como a física prevê (seguindo o "critério de Griffith"). Quando adicionaram hidrogênio, o modelo mostrou que a fenda cresceu muito mais facilmente porque o hidrogênio enfraqueceu a energia superficial. Os resultados coincidiram perfeitamente com as previsões teóricas, provando que o modelo funciona.

2. O Teste Policristalino (A Grande Descoberta): Eles simularam um metal feito de muitos cristais minúsculos (grãos) com contornos entre eles.

   • Sem Hidrogênio: A fenda preferia atravessar os grãos diretamente (fratura transgranular). Era como uma bola de demolição atravessando as paredes de uma casa porque as paredes eram mais fracas do que a argamassa entre elas.
   • Com Hidrogênio: O hidrogênio se acumulou pesadamente nos contornos entre os grãos, enfraquecendo a "argamassa" significativamente mais do que as "paredes". De repente, a fenda mudou seu caminho. Em vez de atravessar os grãos, ela começou a serpentear ao longo dos contornos (fratura intergranular). Foi como se o hidrogênio tivesse transformado a argamassa em areia molhada, fazendo a casa desmoronar pelas costuras em vez de através dos tijolos.

Por Que Isso Importa
Este novo modelo é uma atualização significativa porque não apenas adivinha para onde o hidrogênio vai; ele o calcula com base na termodinâmica real do sistema. Ele captura com sucesso a transição de um tipo de fratura para outro, o que é crucial para entender por que os materiais falham na presença de hidrogênio.

Os autores observam que, embora este modelo seja um grande passo à frente, ele atualmente foca em um mecanismo específico (o hidrogênio enfraquecendo a cola). Trabalhos futuros precisarão adicionar outros fatores complexos, como a forma como o metal dobra e torce (plasticidade) e como outros tipos de defeitos interagem com o hidrogênio. Mas, por enquanto, este modelo fornece uma maneira clara, consistente e precisa de ver como o hidrogênio transforma um metal forte em um frágil.

Resumo técnico

Enunciado do Problema
A fragilização por hidrogênio (HE) representa um desafio crítico para as tecnologias de hidrogênio, causando degradação catastrófica de materiais e fraturas inesperadas em sistemas de energia descarbonizados. Embora mecanismos como a decoesão induzida por hidrogênio (HEDE), a plasticidade localizada induzida por hidrogênio (HELP) e as lacunas induzidas por tensão de hidrogênio (HESIV) tenham sido propostos, prever a HE permanece difícil devido à sua natureza multifísica que abrange as escalas atômica à macroscópica. Os modelos de campo de fase (PF) existentes para o fissuramento assistido por hidrogênio tipicamente acoplam a mecânica de fratura frágil com o mecanismo HEDE via difusão de hidrogênio. No entanto, esses modelos sofrem de duas limitações primárias: eles não reproduzem explicitamente a segregação de hidrogênio (o fenômeno responsável pela redução da energia superficial) e dependem da isoterma de Langmuir-McLean para calcular a cobertura de hidrogênio. Esta última é uma relação de equilíbrio macroscópica entre fases homogêneas de volume e superfície, tornando-a fisicamente injustificável como uma lei local para modelar microestruturas heterogêneas que envolvem interações elásticas.

Metodologia
Os autores propõem um modelo de campo de fase termodinamicamente consistente baseado no formalismo de Kim-Kim-Suzuki (KKS) para simular o fissuramento assistido por hidrogênio em materiais policristalinos. O modelo opera dentro de um arcabouço variacional onde a energia livre total inclui contribuições da energia de superfície de trinca, energia elástica, acoplamento de contorno de grão (GB) e energia livre química.

Principais características metodológicas incluem:

• Variáveis de Campo: O sistema é descrito por um campo de dano ($\eta$), uma probabilidade de ocupação de hidrogênio ($c$) e um conjunto de parâmetros de ordem ($\phi_i$) representando orientações cristalográficas.
• Formulação KKS: Em vez de impor uma lei local de Langmuir-Mclean, o modelo define fases distintas de volume ($c_b$) e de superfície de trinca ($c_{ck}$). A densidade de energia livre química é construída como uma soma ponderada das contribuições de volume e de trinca, restrita pela igualdade dos potenciais químicos ($\mu_b = \mu_{ck}$). Isso permite que as propriedades de interface sejam desacopladas dos campos de concentração.
• Consistência Termodinâmica: O modelo assume soluções sólidas ideais para a energia química e deriva as equações de equilíbrio governantes (equilíbrio mecânico e balanço de potencial químico) a partir do funcional de energia livre total.
• Implementação Numérica: As equações de equilíbrio são resolvidas utilizando métodos espectrais de Fourier e um solver FFT discreto para o equilíbrio mecânico.

Principais Contribuições

1. Modelagem Explícita de Segregação: O modelo contabiliza naturalmente a segregação de hidrogênio em superfícies de trinca e contornos de grão sem assumir uma isoterma macroscópica como uma lei local.
2. Derivação Variacional da Degradação da Energia de Superfície: A degradação efetiva da energia de superfície de trinca induzida pela segregação de hidrogênio é computada analiticamente dentro do arcabouço termodinâmico.
3. Captura de Transição: O modelo demonstra a capacidade de capturar a transição de fissuração transgranular para fissuração intergranular assistida por hidrogênio via mecanismo HEDE.

Resultos

• Validação em Monocristais: Simulações numéricas de uma trinca pré-existente em um monocristal sob carregamento de modo I mostraram que a tensão crítica para o crescimento da trinca ($\sigma_c$) segue o critério de Griffith. Quando o mecanismo HEDE foi ativado, $\sigma_c$ diminuiu devido a uma redução na energia de superfície efetiva ($\Delta \gamma_{ch} \approx -4,29$ J/m²), consistente com as previsões analíticas derivadas das equações do modelo.
• Aplicação Policristalina: Simulações em um material policristalino à base de níquel revelaram uma transição distinta nos caminhos de propagação de trinca. Em um ambiente livre de hidrogênio, a trinca propagou-se transgranularmente devido à alta tenacidade do contorno de grão (GB). Com a introdução de hidrogênio, o modelo previu uma mudança para fissuração intergranular. Isso foi atribuído a uma redução drástica na razão de energia de fratura efetiva GB-volume (de 0,5 para 0,07), causada por energias de segregação de hidrogênio mais altas em superfícies do tipo GB comparadas às superfícies do tipo volume. O mapa de concentração de hidrogênio confirmou a segregação significativa em GBs e superfícies de GB rompidas.

Significância
O artigo afirma que este trabalho fornece um método rigoroso e genérico para modelar o mecanismo HEDE na escala microestrutural. Ao utilizar o formalismo KKS, o modelo reproduz a segregação de hidrogênio e a resultante redução da energia de superfície sem depender da isoterma de Langmuir-McLean, que é termodinamicamente limitada para situações heterogêneas. O estudo valida que a isoterma emerge naturalmente apenas em limites específicos homogêneos e estáticos, justificando assim a nova abordagem para microestruturas complexas. O modelo captura com sucesso a transição observada experimentalmente de fissuração transgranular para intergranular utilizando dados de entrada realistas para o níquel, oferecendo uma ferramenta robusta para compreender a falha assistida por hidrogênio em metais policristalinos. Os autores observam que trabalhos futuros abordarão a anisotropia da energia de superfície e a incorporação de outros mecanismos de HE, como o HELP, através da modelagem de discordâncias.