Thermodynamically consistent phase field model for hydrogen-assisted cracking

Este artículo presenta un modelo de campo de fase termodinámicamente consistente que simula el agrietamiento asistido por hidrógeno en materiales policristalinos mediante el acoplamiento de la propagación de grietas con la segregación de hidrógeno y la reducción de la energía interfacial, capturando con éxito la transición del fallo transgranular al intergranular bajo mecanismos de decohesión mejorada por hidrógeno.

Lo esencial

Imagina que tienes una estructura metálica muy fuerte, como el chasis de un coche o un puente. Esperarías que soporte la presión, pero a veces, átomos de hidrógeno invisibles se cuelan en el metal y hacen que se rompa inesperadamente. Este fenómeno se llama fragilización por hidrógeno. Es como si el metal estuviera siendo "envenenado" secretamente desde el interior, volviéndolo quebradizo y propenso a partirse.

Los científicos han estado intentando construir modelos informáticos para predecir exactamente cómo y dónde se romperá este metal. Sin embargo, los modelos anteriores tenían un fallo importante: trataban el comportamiento del hidrógeno como una regla simple y uniforme que se aplica en todas partes, a pesar de que la estructura interna del metal es en realidad un mosaico complejo de diferentes granos y límites.

El nuevo modelo "inteligente"
Los autores de este artículo han creado una simulación informática más sofisticada (llamada "modelo de campo de fase") que actúa como un mapa de alta definición y termodinámicamente consistente. Así es como funciona, utilizando algunas analogías de la vida cotidiana:

• El metal como una multitud: Imagina que el metal es una habitación llena de gente (átomos de metal). Los "límites de grano" son las líneas invisibles que separan a los distintos grupos de personas. La "grieta" es un hueco que crece en la multitud.
• El hidrógeno como un invitado pegajoso: Los átomos de hidrógeno son como invitados pegajosos a los que les encanta esconderse en los espacios vacíos entre la gente. Tienen una preferencia especial: les encanta pegarse a los bordes de la grieta y a las líneas entre los grupos (límites de grano) incluso más de lo que les gusta el medio de la multitud.
• El problema del "pegamento": En un metal sano, el "pegamento" que mantiene unidos los bordes de la grieta es fuerte. Pero cuando estos invitados pegajosos de hidrógeno se reúnen en los bordes de la grieta, actúan como un aceite resbaladizo, debilitando el pegamento. Esto hace que la grieta sea mucho más fácil de abrir.
• El enfoque antiguo frente al nuevo:
  • Modelos antiguos: Utilizaban un libro de reglas genérico (la isoterma de Langmuir-McLean) que asumía que el hidrógeno estaba distribuido uniformemente y en perfecto equilibrio en todas partes. Esto es como asumir que todo el mundo en la habitación llena de gente está quieto y distribuido uniformemente, lo cual no es cierto cuando se está formando una grieta.
  • Nuevo modelo: Utiliza un marco "variacional" flexible (basado en el formalismo de Kim-Kim-Suzuki). En lugar de imponer una regla rígida, permite que el hidrógeno "migre" naturalmente hacia donde quiere ir (los bordes de la grieta y los límites de grano) basándose en las condiciones locales. Calcula exactamente cuánto se debilita el "pegamento" en tiempo real a medida que el hidrógeno se reúne.

Lo que descubrieron
El equipo probó su nuevo modelo con dos escenarios principales:

1. La prueba de la grieta única: Simularon una grieta en una sola pieza de metal. Sin hidrógeno, la grieta creció exactamente como predice la física (siguiendo el criterio de Griffith). Al añadir hidrógeno, el modelo mostró que la grieta crecía con mucha más facilidad porque el hidrógeno había debilitado la energía superficial. Los resultados coincidieron perfectamente con las predicciones teóricas, demostrando que el modelo funciona.

2. La prueba policristalina (El gran descubrimiento): Simularon un metal compuesto por muchos cristales diminutos (granos) con límites entre ellos.

   • Sin hidrógeno: La grieta prefería atravesar directamente los granos (fractura transgranular). Era como una bola de demolición atravesando las paredes de una casa porque las paredes eran más débiles que el mortero entre ellas.
   • Con hidrógeno: El hidrógeno se acumuló fuertemente en los límites entre los granos, debilitando el "mortero" significativamente más que las "paredes". De repente, la grieta cambió su trayectoria. En lugar de atravesar los granos, empezó a serpentear a lo largo de los límites (fractura intergranular). Fue como si el hidrógeno convirtiera el mortero en arena húmeda, haciendo que la casa se desmoronara por las costuras en lugar de a través de los ladrillos.

Por qué esto es importante
Este nuevo modelo es una mejora significativa porque no se limita a adivinar a dónde va el hidrógeno; lo calcula basándose en la termodinámica real del sistema. Captura con éxito la transición de un tipo de fractura a otro, lo cual es crucial para comprender por qué los materiales fallan en presencia de hidrógeno.

Los autores señalan que, aunque este modelo es un gran paso adelante, actualmente se centra en un mecanismo específico (el debilitamiento del pegamento por el hidrógeno). El trabajo futuro tendrá que añadir otros factores complejos, como la forma en que el metal se dobla y se retuerce (plasticidad) y cómo interactúan otros tipos de defectos con el hidrógeno. Pero por ahora, este modelo proporciona una forma clara, consistente y precisa de ver cómo el hidrógeno convierte un metal fuerte en uno frágil.

Resumen técnico

Planteamiento del Problema
La fragilización por hidrógeno (HE, por sus siglas en inglés) representa un desafío crítico para las tecnologías de hidrógeno, causando una degradación material catastrófica y fracturas inesperadas en los sistemas de energía descarbonizados. Aunque se han propuesto mecanismos como la decohesión mejorada por hidrógeno (HEDE), la plasticidad localizada mejorada por hidrógeno (HELP) y las vacancias inducidas por deformación mejoradas por hidrógeno (HESIV), predecir la HE sigue siendo difícil debido a su naturaleza multifísica que abarca desde la escala atómica hasta la macroscópica. Los modelos de campo de fase (PF) existentes para el agrietamiento asistido por hidrógeno suelen acoplar la mecánica de fractura frágil con el mecanismo HEDE mediante la difusión de hidrógeno. Sin embargo, estos modelos presentan dos limitaciones primarias: no reproducen explícitamente la segregación de hidrógeno (el fenómeno responsable de la reducción de la energía superficial) y dependen de la isoterma de Langmuir-McLean para calcular la cobertura de hidrógeno. Esta última es una relación de equilibrio macroscópica entre fases homogéneas de volumen y superficie, lo que carece de justificación física como ley local para modelar microestructuras heterogéneas que involucran interacciones elásticas.

Metodología
Los autores proponen un modelo de campo de fase termodinámicamente consistente basado en el formalismo de Kim-Kim-Suzuki (KKS) para simular el agrietamiento asistido por hidrógeno en materiales policristalinos. El modelo opera dentro de un marco variacional donde la energía libre total incluye contribuciones de la energía de superficie de la grieta, la energía elástica, el acoplamiento de los bordes de grano (GB) y la energía libre química.

Las características metodológicas clave incluyen:

• Variables de Campo: El sistema se describe mediante un campo de daño ($\eta$), una probabilidad de ocupación de hidrógeno ($c$) y un conjunto de parámetros de orden ($\phi_i$) que representan las orientaciones cristalográficas.
• Formulación KKS: En lugar de imponer una ley local de Langmuir-McLean, el modelo define fases distintas de volumen ($c_b$) y de superficie de grieta ($c_{ck}$). La densidad de energía libre química se construye como una suma ponderada de las contribuciones de volumen y de grieta, restringida por la igualdad de los potenciales químicos ($\mu_b = \mu_{ck}$). Esto permite que las propiedades de la interfaz se desacoplen de los campos de concentración.
• Consistencia Termodinámica: El modelo asume soluciones sólidas ideales para la energía química y deriva las ecuaciones de equilibrio gobernantes (equilibrio mecánico y balance de potencial químico) a partir del funcional de energía libre total.
• Implementación Numérica: Las ecuaciones de equilibrio se resuelven utilizando métodos espectrales de Fourier y un solver de FFT discreto para el equilibrio mecánico.

Contribuciones Clave

1. Modelado Explícito de la Segregación: El modelo contabiliza naturalmente la segregación de hidrógeno en las superficies de grieta y en los bordes de grano sin asumir una isoterma macroscópica como una ley local.
2. Derivación Variacional de la Degradación de la Energía Superficial: La degradación efectiva de la energía de superficie de la grieta inducida por la segregación de hidrógeno se calcula analíticamente dentro del marco termodinámico.
3. Captura de la Transición: El modelo demuestra la capacidad de capturar la transición de un agrietamiento transgranular a uno intergranular asistido por hidrógeno mediante el mecanismo HEDE.

Resultados

• Validación en Monocristales: Las simulaciones numéricas de una grieta preexistente en un monocristal bajo carga de Modo I mostraron que el esfuerzo crítico para el crecimiento de la grieta ($\sigma_c$) sigue el criterio de Griffith. Cuando se activó el mecanismo HEDE, $\sigma_c$ disminuyó debido a una reducción en la energía superficial efectiva ($\Delta \gamma_{ch} \approx -4.29$ J/m$^2$), lo cual es consistente con las predicciones analíticas derivadas de las ecuaciones del modelo.
• Aplicación Policristalina: Las simulaciones en un material policristalino basado en níquel revelaron una transición distintiva en las trayectorias de propagación de grietas. En un entorno libre de hidrógeno, la grieta se propagó de forma transgranular debido a la alta tenacidad de los bordes de grano (GB). Tras la introducción de hidrógeno, el modelo predijo un cambio hacia el agrietamiento intergranular. Esto se atribuyó a una reducción drástica en la relación de energía de fractura de GB frente al volumen (de 0.5 a 0.07), causada por energías de segregación de hidrógeno más altas en las superficies de tipo GB en comparación con las de tipo volumen. El mapa de concentración de hidrógeno confirmó una segregación significativa en los GB y en las superficies de grietas de los GB.

Significancia
El artículo afirma que este trabajo proporciona un método riguroso y genérico para modelar el mecanismo HEDE a escala microestructural. Al utilizar el formalismo KKS, el modelo reproduce la segregación de hidrógeno y la resultante reducción de la energía superficial sin depender de la isoterma de Langmuir-McLean, la cual está limitada termodinámicamente para situaciones heterogéneas. El estudio valida que la isoterma emerge naturalmente solo en límites específicos homogéneos y estáticos, justificando así el nuevo enfoque para microestructuras complejas. El modelo captura con éxito la transición observada experimentalmente de un agrietamiento transgranular a uno intergranular utilizando datos de entrada realistas para el níquel, ofreciendo una herramienta robusta para comprender la falla asistida por hidrógeno en metales policristalinos. Los autores señalan que el trabajo futuro abordará la anisotropía de la energía superficial y la incorporación de otros mecanismos de HE, como el HELP, mediante el modelado de dislocaciones.