Thermodynamically consistent phase field model for hydrogen-assisted cracking

Diese Arbeit präsentiert ein thermodynamisch konsistentes Phasenfeldmodell, das wasserstoffunterstützte Rissbildung in polykristallinen Materialien simuliert, indem es die Rissausbreitung mit Wasserstoffsegregation und der Reduktion der Grenzflächenenergie koppelt und dabei den Übergang von transkristallinem zu interkristallinem Versagen unter Mechanismen der wasserstoffverstärkten Dekohäsion erfolgreich erfasst.

Das Wesentliche

Stellen Sie sich vor, Sie haben eine sehr starke Metallstruktur, wie etwa den Rahmen eines Autos oder einer Brücke. Man würde erwarten, dass sie unter Druck standhält, aber manchmal schlüpfen unsichtbare Wasserstoffatome in das Metall und lassen es unerwartet zerbrechen. Dieses Phänomen wird als Wasserstoffversprödung bezeichnet. Es ist, als würde das Metall im Geheimen von innen heraus „vergiftet“, was es spröde macht und dazu neigt, plötzlich zu brechen.

Wissenschaftler haben versucht, Computermodelle zu entwickeln, die genau vorhersagen, wie und wo dieses Metall brechen wird. Das bisherige Problem der Modelle war jedoch ein schwerwiegender Fehler: Sie behandelten das Verhalten des Wasserstoffs wie eine einfache, einheitliche Regel, die überall gilt, obwohl die interne Struktur des Metalls in Wirklichkeit ein komplexes Patchwork-Quilt aus verschiedenen Körnern und Korngrenzen ist.

Das neue „intelligente“ Modell
Die Autoren dieser Arbeit haben eine neue, anspruchsvollere Computersimulation (ein sogenanntes „Phasenfeldmodell“) entwickelt, die wie eine hochauflösende, thermodynamisch konsistente Karte fungiert. So funktioniert es, unter Verwendung einiger Alltagsanalogien:

• Das Metall als Menschenmenge: Stellen Sie sich das Metall als einen überfüllten Raum voller Menschen (Metallatome) vor. Die „Korngrenzen“ sind die unsichtbaren Linien, die verschiedene Gruppen von Menschen voneinander trennen. Der „Riss“ ist eine wachsende Lücke in der Menge.
• Der Wasserstoff als klebriger Gast: Wasserstoffatome sind wie klebrige Gäste, die es lieben, sich in den leeren Zwischenräumen zwischen den Menschen zu verstecken. Sie haben eine besondere Vorlieền: Sie lieben es, sich noch mehr an die Ränder des Risses und an die Linien zwischen den Gruppen (Korngrenzen) zu heften als in der Mitte der Menge.
• Das „Kleber“-Problem: In einem gesunden Metall ist der „Kleber“, der die Risskanten zusammenhält, stark. Wenn sich diese klebrigen Wasserstoffgäste jedoch an den Rändern des Risses sammeln, wirken sie wie ein schmieriges Öl, das den Kleber schwächt. Dies macht es viel einfacher, den Riss zu öffnen.
• Der alte vs. der neue Ansatz:
  • Alte Modelle: Nutzten ein generisches Regelwerk (die Langmuir-McLean-Isotherme), das davon ausging, dass der Wasserstoff gleichmäßig verteilt und überall im perfekten Gleichgewicht ist. Das ist so, als würde man davon ausgehen, dass alle Menschen in dem überfüllten Raum stillstehen und gleichmäßig verteilt sind, was nicht der Fall ist, wenn sich ein Riss bildet.
  • Neues Modell: Verwendet ein flexibles, „variationales“ Framework (basierend auf dem Kim-Kim-Suzuki-Formalismus). Anstatt eine starre Regel aufzuzwingen, lässt es den Wasserstoff natürlich dorthin „wandern“, wo er hin möchte (an die Risskanten und Korngrenzen), basierend auf den lokalen Bedingungen. Es berechnet exakt, wie sehr der „Kleber“ in Echtzeit geschwächt wird, während sich der Wasserstoff ansammelt.

Was sie herausgefunden haben
Das Team hat ihr neues Modell mit zwei Hauptszenarien getestet:

1. Der Einzelriss-Test: Sie simulierten einen Riss in einem einzelnen Stück Metall. Oh生 Wasserstoff wuchs der Riss genau so, wie es die Physik vorhersagt (gemäß dem Griffith-Kriterium). Als sie Wasserstoff hinzufügten, zeigte das Modell, dass der Riss viel leichter wuchs, weil der Wasserstoff die Oberflächenenergie geschwächt hatte. Die Ergebnisse stimmten perfekt mit den theoretischen Vorhersagen überein, was bewies, dass das Modell funktioniert.

2. Der polykristalline Test (Die große Entdeckung): Sie simulierten ein Metall, das aus vielen winzigen Kristallen (Körnern) mit Grenzen zwischen ihnen besteht.

   • Ohne Wasserstoff: Der Riss bevorzugte es, direkt durch die Körner zu brechen (transgranularer Rissbruch). Es war wie ein Abrissbagger, der durch die Wände eines Hauses bricht, weil die Wände schwächer sind als der Mörtel dazwischen.
   • Mit Wasserstoff: Der Wasserstoff sammelte sich massiv an den Grenzen zwischen den Körnern und schwächte den „Mörtel“ deutlich stärker als die „Wände“. Plötzlich änderte der Riss seinen Pfad. Anstatt durch die Körner zu brechen, begann er, entlang der Grenzen zu mäandern (intergranularer Rissbruch). Es war, als hätte der Wasserstoff den Mörtel in nassen Sand verwandelt, was dazu führte, dass das Haus entlang der Nähte auseinanderfiel, anstatt durch die Ziegel zu brechen.

Warum dies wichtig ist
Dieses neue Modell ist ein bedeutendes Upgrade, da es nicht nur rät, wohin der Wasserstoff geht, sondern es basierend auf der tatsächlichen Thermodynamik des Systems berechnet. Es erfasst erfolgreich den Übergang von einer Art des Rissens zu einer anderen, was entscheidend für das Verständnis dessen ist, warum Materialien in Gegenwart von Wasserstoff versagen.

Die Autoren merken an, dass dieses Modell zwar ein großer Schritt nach vorn ist, sich aber derzeit auf einen spezifischen Mechanismus konzentriert (die Schwächung des Klebers durch Wasserstoff). Zukünftige Arbeiten werden weitere komplexe Faktoren hinzufügen müssen, wie etwa die Verformung und Verdrehung des Metalls (Plastizität) und wie andere Arten von Defekten mit dem Wasserstoff interagieren. Aber für den Moment bietet dieses Modell eine klare, konsistente und genaue Möglichkeit zu sehen, wie Wasserstoff ein starkes Metall in ein zerbrechliches verwandelt.

Technische Zusammenfassung

Problemstellung
Wasserstoffversprödung (Hydrogen Embrittlement, HE) stellt eine kritische Herausforderung für Wasserstofftechnologien dar, da sie zu katastrophalem Materialabbau und unerwartetem Bruch in dekarbonisierten Energiesystemen führt. Obwohl Mechanismen wie die wasserstoffunterstützte Dekohäsion (HEDE), die wasserstoffunterstützte lokalisierte Plastizität (HELP) und die wasserstoffunterstützte durch Dehnung induzierte Vakanzenbildung (HESIV) vorgeschlagen wurden, bleibt die Vorhersage von HE aufgrund ihrer Multiphysik-Natur, die von der atomaren bis zur makroskopischen Skala reicht, schwierig. Bestehende Phasenfeld-Modelle (PF) für wasserstoffunterstützte Rissbildung koppeln typischerweise spröde Bruchmechanik mit dem HEDE-Mechanismus über die Wasserstoffdiffusion. Diese Modelle weisen jedoch zwei primäre Einschränkungen auf: Sie bilden die Wasserstoffsegregation (das Phänomen, das für die reduzierte Oberflächenenergie verantwortlich ist) nicht explizit ab und stützen sich auf die Langmuir-McLean-Isotherme zur Berechnung der Wasserstoffbelegung. Letztere ist eine makroskopische Gleichgewichtsbeziehung zwischen homogenen Bulk- und Oberflächenphasen, was sie als lokales Gesetz für die Modellierung heterogener Mikrostrukturen unter Berücksichtigung elastischer Wechselwirkungen physikalisch ungerechtfertigt macht.

Methodik
Die Autoren schlagen ein thermodynamisch konsistentes Phasenfeldmodell basierend auf dem Kim-Kim-Suzuki-Formalismus (KKS) vor, um die wasserstoffunterstützte Rissbildung in polykristallinen Materialien zu simulieren. Das Modell operiert innerhalb eines Variationsrahmens, bei dem die gesamte freie Energie Beiträge aus der Rissflächenenergie, der elastischen Energie, der Korngrenzenkopplung (GB-Kopplung) und der chemischen freien Energie umfasst.

Zentrale methodische Merkmale sind:

• Feldvariablen: Das System wird durch ein Schadensfeld ($\eta$), eine Wasserstoffbelegungswahrscheinlichkeit ($c$) und einen Satz von Ordnungsparametern ($\phi_i$) beschrieben, die kristallographische Orientierungen repräsentieren.
• KKS-Formalismus: Anstatt ein lokales Langmuir-McLean-Gesetz aufzuzwingen, definiert das Modell distinkte Bulk- ($c_b$) und Rissflächenphasen ($c_{ck}$). Die chemische Energiedichte wird als gewichtete Summe von Bulk- und Rissbeiträgen konstruiert, wobei die Gleichheit der chemischen Potentiale ($\mu_b = \mu_{ck}$) als Nebenbedingung dient. Dies ermöglicht die Entkopplung der Grenzflächeneigenschaften von den Konzentrationsfeldern.
• Thermodynamische Konsistenz: Das Modell nimmt ideale feste Lösungen für die chemische Energie an und leitet die zugrunde liegenden Gleichgewichtsgleichungen (mechanisches Gleichgewicht und chemisches Potentialgleichgewicht) aus dem gesamten freien Energiefunktional ab.
• Numerische Implementierung: Die Gleichgewichtsgleichungen werden mittels Fourier-Spektralmethoden und einem diskreten FFT-Solver für das mechanische Gleichgewicht gelöst.

Wesentliche Beiträge

1. Explizite Modellierung der Segregation: Das Modell berücksichtigt die natürliche Wasserstoffsegregation an Rissflächen und Korngrenzen, ohne eine makroskopische Isotherme als lokales Gesetz vorauszusetzen.
2. Variationale Ableitung der Oberflächenenergie-Degradation: Die effektive Degradation der Rissflächenenergie, induziert durch Wasserstoffsegregation, wird analytisch innerhalb des thermodynamischen Rahmens berechnet.
3. Erfassung des Übergangs: Das Modell demonstriert die Fähigkeit, den Übergang von transkristalliner zu wasserstoffunterstützter interkristalliner Rissbildung über den HEDE-Mechanismus zu erfassen.

Ergebnisse

• Validierung in Einkristallen: Numerische Simulationen eines präexistierenden Risses in einem Einkristall unter Mode-I-Belastung zeigten, dass die kritische Spannung für das Risswachstum ($\sigma_c$) dem Griffith-Kriterium folgt. Als der HEDE-Mechanismus aktiviert wurde, sank $\sigma_c$ aufgrund einer Reduktion der effektiven Oberflächenenergie ($\Delta \gamma_{ch} \approx -4,29$ J/m$^2$), was mit den analytischen Vorhersagen aus den Gleichungen des Modells übereinstimmt.
• Polykristalline Anwendung: Simulationen in einem Nickel-basierten polykristallinen Material offenbarten einen deutlichen Übergang in den Rissausbreitungspfaden. In einer wasserstofffreien Umgebung propagierte der Riss transkristallin aufgrund hoher Korngrenzenzähigkeit. Nach Einführung von Wasserstoff sagte das Modell einen Wechsel zu interkristalliner Rissbildung voraus. Dies wurde einer drastischen Reduktion des effektiven Verhältnisses der Korngrenzen-zu-Bulk-Bruchenergie (von 0,5 auf 0,07) zugeschrieben, verursacht durch höhere Wasserstoffsegregationsenergien an Korngrenzentypen im Vergleich zu Bulk-Typen. Die Wasserstoffkonzentrationskarte bestätigte eine signifikante Segregation an Korngrenzen und gebrochenen Korngrenzenflächen.

Bedeutung
Das Paper behauptet, dass diese Arbeit eine rigorose und generische Methode zur Modellierung des HEDE-Mechanismus auf der Mikrostrukturebene bietet. Durch die Verwendung des KKS-Formalismus reproduziert das Modell die Wasserstoffsegregation und die resultierende Reduktion der Oberflächenenergie, ohne sich auf die thermodynamisch begrenzte Langmuir-McLean-Isotherme für heterogene Situationen zu verlassen. Die Studie validiert, dass die Isotherme natürlich nur in spezifischen homogenen, statischen Grenzfällen hervortritt, wodurch der neue Ansatz für komplexe Mikrostrukturen gerechtfertigt wird. Das Modell erfasst erfolgreich den experimentell beobachteten Übergang von transkristalliner zu interkristalliner Rissbildung unter Verwendung realistischer Eingangsdaten für Nickel und bietet somit ein robustes Werkzeug zum Verständnis des wasserstoffunterstützten Versagens in polykristallinen Metallen. Die Autoren merken an, dass zukünftige Arbeiten die Oberflächenenergieanisotropie sowie die Einbeziehung anderer HE-Mechanismen, wie etwa HELP durch Versetzungsmodellierung, adressieren werden.