Thermodynamically consistent phase field model for hydrogen-assisted cracking

本論文は、亀裂の進展を水素偏析および界面エネルギーの減少と結合させることにより、多結晶材料における水素脆化割れをシミュレートする熱力学的に整合したフェーズフィールドモデルを提示しており、水素誘起脱結合メカニズム下での粒内破壊から粒界破壊への遷移を捉えることに成功している。

要約

非常に強固な金属構造、例えば車のフレームや橋のようなものを想像してみてください。それらは圧力に耐えられるはずだと期待されますが、時として、目に見えない水素原子が金属の内部に忍び込み、予期せぬ破壊を引き起こすことがあります。この現象は「水素脆化（すいそぜいか）」と呼ばれます。これは、金属が内側から密かに「毒されている」ようなもので、金属を脆くし、折れやすくしてしまうのです。

科学者たちは、金属がまさにいつ、どこで壊れるのかを予測するための高度なコンピュータモデルを構築しようと試みてきました。しかし、従来のモデルには大きな欠陥がありました。それは、水素の振る舞いを、単純で一様なルールとして扱っていたことです。実際には、金属の内部構造は異なる結晶粒や粒界（結晶の境界）が入り混じった複雑なパッチワークのキルトのようなものなのですが、従来のモデルではそこが考慮されていませんでした。

新しい「スマート」なモデル
この論文の著者たちは、より洗練された新しいコンピュータ・シミュレーション（「フェーズフィールドモデル」と呼ばれる）を作成しました。これは、熱力学的に整合性の取れた、高精細な地図のように機能します。ここでは、日常的な例えを用いてその仕組みを説明します。

• 群衆としての金属： 金属を、人々（金属原子）で満たされた混雑した部屋だと想像してください。「結晶粒界」は、異なるグループを隔てる目に見えない線のことです。「亀裂（クラック）」は、その群衆の中に広がる隙間です。
• 粘着質なゲストとしての水素： 水素原子は、人々の間の空いたスペースに隠れ込むのが大好きな、粘着質なゲストのようなものです。彼らには特別な好みがあります。彼らは、群衆の真ん中にいるよりも、亀裂の端やグループ間の境界線に張り付くことをより好みます。
• 「接着剤」の問題： 健康な金属では、亀裂の端を繋ぎ止めている「接着剤」は強力です。しかし、この粘着質な水素ゲストが亀裂の端に集まると、彼らは滑りやすいオイルのように振る舞い、接着剤を弱めてしまいます。これにより、亀裂が開きやすくなるのです。
• 旧来のアプローチ vs 新しいアプローチ：
  • 旧モデル： 一般的なルールブック（ラングミュア・マクリーン吸着等温式）を使用しており、水素が至る所で均一に分布し、完璧なバランス状態にあると仮定していました。これは、亀裂が形成されている最中であるにもかかわらず、混雑した部屋の全員が静止して均等に間隔を空けて立っていると仮定するようなものです。
  • 新モデル： 柔軟な「変分」フレームワーク（キム・キム・鈴木形式に基づく）を使用しています。硬直したルールを押し付けるのではなく、局所的な条件に基づいて、水素が自然に「移動」したい場所（亀裂の端や粒界）へと向かうように設計されています。このモデルは、水素が集まるにつれて、どれだけリアルタイムで「接着剤」が弱まるかを正確に計算します。

彼らが発見したこと
チームは、新しいモデルを用いて主に2つのシナリオでテストを行いました。

1. 単一亀裂テスト： 単一の金属片における亀裂をシミュレートしました。水素がない場合、亀裂は物理学の予測通り（グリフィスの基準に従って）成長しました。水素を加えると、モデルは水素が表面エネルギーを低下させたために、亀裂がより容易に成長することを示しました。結果は理論的な予測と完璧に一致しており、モデルが機能することを証明しました。

2. 多結晶テスト（大きな発見）： 多くの小さな結晶（結晶粒）と、その間の境界を持つ金属をシミュレートしました。

   • 水素がない場合： 亀裂は結晶粒を真っ直ぐ突き破ることを好みました（粒内破壊）。これは、壁よりもモルタルの方が弱い場合、鉄球が家の壁を粉砕していくようなものです。
   • 水素がある場合： 水素が結晶粒の境界に大量に集まり、「モルタル」を周囲の「壁」よりも大幅に弱めました。突然、亀裂の経路が変わりました。亀裂は結晶粒を突き破るのではなく、境界に沿って蛇行するように進み始めました（粒界破壊）。これは、まるで水素がモルタルを濡れた砂に変えてしまい、レンガそのものを壊すのではなく、継ぎ目に沿って家が崩壊するようにさせたのです。

なぜこれが重要なのか
この新しいモデルは、単に水素がどこへ行くかを推測するのではなく、システムの実際の熱力学に基づいてそれを計算するため、大幅なアップグレードと言えます。このモデルは、亀裂の種類の遷移（変化）を捉えることに成功しており、これは材料が水素の存在下でなぜ破壊されるのかを理解する上で極めて重要です。

著者らは、このモデルは大きな一歩ではあるものの、現在は一つの特定のメカニズム（水素による接着剤の弱体化）に焦点を当てていると述べています。将来的には、金属が曲がったりねじれたりする性質（塑性）や、他の種類の欠陥と水素がどのように相互作用するかといった、より複雑な要因を追加する必要があります。しかし、現時点において、このモデルは、水素がいかにして強固な金属を脆いものへと変えてしまうのかを、明確かつ一貫した正確な方法で可視化しています。

技術概要

問題提起
水素脆化（HE）は、水素技術における重大な課題であり、脱炭素化されたエネルギーシステムにおいて壊滅的な材料劣化や予期せぬ破壊を引き起こす。水素誘起のデコヒージョン（HEDE）、水素誘起局所塑性（HELP）、および水素誘起歪誘起空孔（HESIV）といったメカニズムが提案されているが、原子スケールからマクロスケールにわたるマルチフィジックスな性質のため、水素脆化の予測は依然として困難である。既存のフェーズフィールド（PF）モデルは、水素拡散を介して脆性破壊力学とHEDEメカニズムを結合させている。しかし、これらのモデルには主に2つの限界がある。第一に、表面エネルギーの低下の原因となる現象である水素偏析を明示的に再現できていないこと。第二に、水素被覆率を算出するためにラングミュア・マクリーン（Langmuir-McLean）吸着等温式に依存していることである。後者は、均質なバルク相と表面相の間のマクロな平衡関係であり、弾性相互作用を含む不均質な微細構造をモデリングするための局所的な法則としては、物理的に正当化できない。

手法
著者らは、多結晶材料における水素誘戯破壊をシミュレートするために、Kim-Kim-Suzuki（KKS）形式に基づく熱力学的に整合したフェーズフィールドモデルを提案している。このモデルは、クラック表面エネルギー、弾性エネルギー、結晶粒界（GB）結合、および化学自由エネルギーの寄与を含む全自由エネルギーを用いた変分フレームワーク内で動作する。

主な手法の特徴は以下の通りである：

• 場変数： 系は、損傷場（$\eta$）、水素占有確率（$c$）、および結晶方位を表す一連の秩序パラメータ（$\phi_i$）によって記述される。
• KKS定式化： 局所的なラングミュア・マクリーン則を課す代わりに、本モデルは異なるバルク相（$c_b$）とクラック表面相（$c_{ck}$）を定義する。化学自由エネルギー密度は、バルクとクラックの寄与の加重和として構成され、化学ポテンシャルの等価性（$\mu_b = \mu_{ck}$）によって制約される。これにより、界面特性を濃度場から切り離すことが可能となる。
• 熱力学的整合性： 本モデルは化学エネルギーに対して理想溶液を仮定し、全自由エネルギー汎関数から支配的な平衡方程式（力学的平衡および化学ポテンシャルのバランス）を導出する。
• 数値実装： 平衡方程式は、フーリエ・スペクトル法および力学的平衡のための離散FFTソルバーを用いて解かれる。

主要な貢献

1. 偏析の明示的モデリング： 本モデルは、マクロな等温式を局所的な法則として仮定することなく、クラック表面や結晶粒界における水素偏析を自然に考慮する。
2. 表面エネルギー低下の変分導出： 水素偏析によって誘起されるクラック表面エネルギーの実効的な低下は、熱力学的枠組みの中で解析的に計算される。
3. 遷移の捕捉： 本モデルは、HEDEメカニズムを介した粒内破壊から水素誘起粒界破壊への遷移を捉える能力を示している。

結果

• 単結晶における検証： モードI荷重下における単結晶内の既存クラックの数値シミュレーションでは、クラック進展の臨界応力（$\sigma_c$）がグリフィスの基準に従うことが示された。HEDEメカニチズムが活性化されると、実効表面エネルギーの減少（$\Delta \gamma_{ch} \approx -4.29$ J/m$^2$）により$\sigma_c$が低下し、これはモデルの方程式から導出された解析的予測と一致した。
• 多結晶への適用： ニッケル基多結晶材料におけるシミュレーションにより、クラック進展経路の明確な遷移が明らかになった。水素のない環境では、高い粒界靭性によりクラックは粒内進展した。水素を導入すると、モデルは粒界破壊への移行を予測した。これは、バルク型の表面と比較してGB型の表面における水素偏析エネルギーが高いために、実効的な粒界対バルクの破壊エネルギー比が劇的に減少した（0.5から0.07へ）ことに起因する。水素濃度マップは、GBおよび破断したGB表面における顕著な偏析を確認した。

意義
本論文は、微細構造スケールにおけるHEDEメカニズムをモデリングするための厳密かつ汎用的な手法を提供すると主張している。KKS形式を利用することで、本モデルは、不均質な状況において熱力学的に制限のあるラングミュア・マクリーン吸着等温式に頼ることなく、水素偏析とそれに伴う表面エネルギーの低下を再現している。本研究は、等温式が特定の均質かつ静的な極限においてのみ自然に現れることを実証しており、複雑な微細構造に対する新しいアプローチの正当性を裏付けている。モデルは、ニッケルの現実的な入力データを用いて、実験的に観察される粒内から粒界への破壊遷移を捉えることに成功した。著者らは、今後の課題として、表面エネルギーの異方性と、転位モデリングを通じたHELPなどの他のHEメカニズムの組み込みを挙げている。