吸着熱の第一原理計算データを用いたMOF水素吸蔵タンクのマルチスケール熱設計解析#

Advance/PHASE Logo

水素社会の実現に向けた課題の一つは、高密度かつ安全な水素貯蔵技術の確立です。金属有機構造体（MOF：Metal Organic Frameworks）は多孔質構造により優れた水素吸蔵能力を持ちますが、吸着時に発生する「吸着熱」が材料自体の温度を上昇させ、吸蔵効率を著しく低下させるという熱管理上の課題があります。本事例では、第一原理計算ソフトウェアAdvance/PHASEによって得られた微視的な吸着エネルギーを、マクロの1次元熱化学反応器モデルへ橋渡しするマルチスケール連携手法を紹介します。これにより、充填中のタンク内部に生じる「ホットスポット」の温度分布と、それに伴う実効的な水素吸蔵容量のロスを定量的に算出することが可能となります。

Keywords: 第一原理計算（DFT）, Metal-Organic Framework (MOF), 水素吸蔵, マルチスケール連携, 1次元熱化学反応器モデル, 熱管理, ホットスポット解析

解析モデルと理論背景#

1. Advance/PHASEを用いた第一原理計算：吸着エネルギーの精密評価#

解析の第一段階として、第一原理計算ソフトウェアAdvance/PHASEを用い、MOFの一種である Ni₂(m-dobdc) のオープンメタルサイトに対する水素分子の吸着状態をシミュレーションしました。水素吸蔵においては、分子が化学結合せず物理吸着するため、弱い相互作用を記述できるvan der Waals補正（DFT-D3）を考慮した計算を行っています。解析の結果、水素分子は解離することなく分子状で安定化し、その吸着エンタルピー変化は $\Delta H = -11.53$ kJ/mol （発熱反応、吸着熱: 11.53 kJ/mol）と算出されました [1]。この値は、後述するマクロモデルにおいて、タンク内で発生する熱量の根拠となる重要なパラメータです。

2. Tri-site Langmuirモデル：ミクロとメゾの橋渡し#

MOF内部の吸着挙動をより現実に即して記述するため、中性子回折実験および理論計算に基づく「Tri-site Langmuirモデル」を採用しました [2]。これは、MOF内の異なる3つの吸着サイト（Site A, B, C）を区別して扱うモデルです。全吸着量 $q_\text{total}$ は、各サイトの独立した吸着量の総和として以下のように表されます。

$q_\text{total} = \sum_{i \in \{A, B, C\}} \frac{q_{\text{sat},i} \cdot b_i(T) \cdot P}{1 + b_i(T) \cdot P}, \quad\quad b_i(T) = b_{0,i} \exp\left(\frac{E_i}{R T}\right)$

ここで、 $q_{\text{sat},i}$ は飽和吸着量、 $b_i$ は吸着平衡定数です。本解析のハイライトは、最も強く水素を引き寄せるSite A（オープンメタルサイト）の吸着熱 $E_A$ として、Advance/PHASEの計算値（ $11.53$ kJ/mol）を直接適用している点にあります。これにより、第一原理計算結果をシステム設計に活かすことが可能になります。

図1. Tri-siteモデルの模式図。左：MOF細孔内におけるSite A（Ni原子直上）、Site B（リンカー壁面）、Site C（細孔中心）の物理的配置。右：各サイトの吸着エンタルピー変化。Site AにAdvance/PHASEの結果を、Site BとSite Cに文献値 [2] を適用しています。

3. 1次元熱化学反応器モデル：マクロスケールへの展開#

タンク内部の巨視的な温度・濃度分布を予測するため、半径方向への1次元熱伝導および吸着速度を考慮したマクロモデルを構築しました。このモデルは、多孔質体充填層（Packed Bed）内の熱物質移動現象を記述する標準的な手法に基づいており、XiaoらによるMOF-5タンクのシミュレーション研究 [3] を参考に、円筒座標系における質量保存則とエネルギー保存則を連立させています。

$\rho C_p \frac{\partial T}{\partial t} = \frac{1}{r} \frac{\partial}{\partial r} \left( r \lambda \frac{\partial T}{\partial r} \right) - \sum_i \left( \frac{dq_i}{dt} \cdot \Delta H_i \right)$

さらに本解析では、より現実的な熱設計を評価するため、タンク壁面においてRobin境界条件（対流熱伝達）を適用しました。これは「壁温一定」という理想条件ではなく、外部を流れる冷媒（液体窒素等）とタンク壁との間の熱の移動抵抗を考慮したものです。

$-\lambda \left. \frac{\partial T}{\partial r} \right|_{r=R} = h (T_\text{surface} - T_\text{coolant})$

ここで、 $h$ は熱伝達係数、 $T_\text{coolant}$ は冷媒温度です。これにより、冷却性能が不足した場合の影響も評価可能となります。なお、本解析ではKumarらによる1次元簡略化モデルのアプローチ [4] を参考に、現象の本質（半径方向の熱移動律速）に着目して計算コストを抑えた1次元円筒座標モデルを採用しています。また、吸着速度式には計算安定性と物理精度のバランスに優れた線形推進力（LDF）モデルを用いました。図2にこのモデルの離散化概念を示します。

1D Reactorモデルの模式図
図2. タンク半径方向の離散化モデル。中心から壁面（冷却境界）までの温度勾配と、各ノード内での独立した吸着・発熱を同時に計算します。タンク内部の熱流動現象において、MOFの熱伝導率が非常に低く、圧力勾配による流動タイムスケールは熱拡散に比べて無視できるほど小さいため、空間均一圧力近似に基づく1次元非定常熱物質収支モデル（1D Reactor Model）を採用しました。壁面には外部冷媒との熱伝達（Robin境界）を設定しました。

4. タンクの仕様と運転条件#

実用的な水素ストレージシステムを想定し、以下の幾何学的寸法および運転条件を設定しました。これらの条件は、熱伝導の遅れによるホットスポット現象を顕著に観察するためのベンチマークとして設定されています。熱伝導率については、MOF-5単結晶の測定値 0.32 W/(m·K) [5] と、粉末成形体の測定値 0.088 W/(m·K) [3] を参照し、粒子間接触や充填状態を考慮した実効的な設計値として 0.15 W/(m·K) を採用しました。

表1. 解析に使用した主要パラメータ設定

カテゴリー	パラメータ	設定値と物理的根拠（出典）
吸着エネルギー	Site A 吸着熱 ( $E_A$ )	11.53 kJ/mol (Advance/PHASE計算値 [1])
吸着容量	Site A 飽和容量 ( $q_{\text{sat},A}$ )	3.78 mol/kg (文献 [2] のRietveld解析値)
	総飽和容量 ( $q_\text{total}$ )	11.97 mol/kg (文献 [2] の実験フィッティング値)
MOF物性	充填密度 ( $\rho_\text{mof}$ )	760 kg/m³ (結晶密度1.17 g/cm³ [2] に充填率を考慮)
	熱伝導率 ( $\lambda$ )	0.15 W/(m·K) (MOFペレットの一般的な代表値)
幾何形状	タンク半径 ( $R_\text{tank}$ )	10.0 cm (直径20cmの円筒型試験タンクを想定)
	解析ドメイン	半径方向 1次元 (中心 $r=0$ から壁面 $r=R$ )
冷却・境界条件	冷媒温度 ( $T_\text{coolant}$ )	77 K (液体窒素温度)
	熱伝達係数 ( $h$ )	200 W/(m²·K) (強制対流冷却を想定)
	充填圧力 / 時間	1.0 MPa (10 bar) / 300 秒

解析結果と考察#

1. 77K充填時のホットスポット形成と冷却限界#

シミュレーションの結果、タンク外部を77Kの冷媒で冷却しているにもかかわらず、中心部の温度は最大で140K近くまで上昇することが明らかになりました（図3）。これは、MOFの熱伝導率が極めて低い（ $\lambda=0.15$ ）ため、内部で発生した吸着熱が壁面まで到達するのに時間がかかり、熱が内部に閉じ込められてしまう「ホットスポット現象」です。この傾向は、Chakrabortyらによるクライオ吸着水素貯蔵システムの熱管理シミュレーション [6] において、低熱伝導率の吸着材では中心部の温度上昇が顕著になるという報告ともよく一致しています。Robin境界条件を適用したことで、壁面温度自体も冷媒温度（77K）より若干上昇しており、現実的な冷却の限界が見て取れます。

半径方向温度分布
図3. タンク内半径方向の温度分布推移（Robin境界条件）。壁面（r=10cm）においても77Kとの温度乖離が見られ、冷却の遅れが確認できます。

2. 熱による「見えないロス」の定量化#

図4は、理想的な等温条件（77K一定）と比較した実際の水素吸蔵量の推移です。グラフ内のグレーの領域は、熱管理が不十分なために失われた吸蔵容量（Thermal Loss）を表しています。最終的な充填量は理想値の約67%に留まっており、高性能な材料を使用しても、熱設計が追いつかなければそのポテンシャルを30%以上無駄にすることが定量的に明らかになりました。

77K充填時の吸着量推移比較
図4. 理想状態（青）と発熱考慮時（赤）の吸着量推移比較。温度上昇により充填効率が著しく低下しています。

本解析アプローチの適用範囲と限界について#

吸着等温線の取り扱い：
本解析では、吸着等温線そのものをモデル内で都度シミュレーション（GCMC法など）しているわけではありません。事前にTri-site Langmuirモデルという「関数」を定義し、その主要パラメータ（Site Aの吸着熱）に第一原理計算の結果を適用しています。したがって、解析精度は適用した吸着モデル（Langmuir式）の妥当性と、入力パラメータの精度に依存します。
物理モデルの簡略化：
計算コストと実用性のバランスを考慮し、吸着速度には線形推進力（LDF）モデルを採用し、タンク形状は1次元円筒モデルへと縮退させています。MOF結晶内の詳細な拡散挙動や、タンク端部での3次元的な熱流動効果を厳密に評価する必要がある場合は、より高次な3Dモデルへの拡張が必要です。
材料探索への適用：
本手法は、新しいMOF材料の「吸着熱」や「細孔容量」が第一原理計算等で予測できれば、その材料をタンクに用いた際のシステム性能（温度上昇や充填効率）を即座に予測可能です。材料スクリーニングの初期段階において有効なツールとなります。

より複雑・詳細な熱流体シミュレーションへ
本事例では、物理現象の理解やプロトタイピング（Concept Proof）には、Pythonソルバーによる1次元熱化学反応器モデル（簡易熱流動解析）を行いました。これは実際のタンク設計において、より複雑・詳細なシミュレーションを行う必要があり、Advance/FrontFlow/redのような流体解析ソフトウェアの使用が推奨されます。

まとめ#

本事例では、Advance/PHASEによる第一原理計算（ミクロ）から得られた「吸着熱」という材料固有のパラメータを、1次元熱化学反応器モデル（マクロ）へと連携させることで、水素貯蔵システムの熱課題を定量化しました。特に、「材料は高性能（吸着熱が強い）であるほど、システム化の際には発熱がボトルネックとなり性能をフルに発揮できない」という、材料開発とデバイス設計の間のジレンマが、シミュレーションによる可視化を通じて浮き彫りになりました。この結果から、熱伝導率向上のためのフィン設計や充填プロトコルの最適化など、材料性能を最大限に活かすためのデバイス設計の重要性が示唆されます。

本解析の詳細や、研究への適用可能性に関するご相談はこちら

お問い合わせ

参考文献#

MOFへの水素吸蔵シミュレーション
M. T. Kapelewski et al., "M2(m-dobdc) (M = Mg, Mn, Fe, Co, Ni) Metal-Organic Frameworks Exhibiting Increased Charge Density and Enhanced H2 Binding at the Open Metal Sites", J. Am. Chem. Soc. 136, 12119 (2014).
J. Xiao, M. Hu, P. Bénard, and R. Chahine, "Simulation of hydrogen storage tank packed with metal-organic framework", Int. J. Hydrog. Energy 38, 13000 (2013).
V. S. Kumar and S. Kumar, "Generalized model development for a cryo-adsorber and 1-D results for the isobaric refueling period", Int. J. Hydrog. Energy 35, 3598 (2010).
B.L. Huang, Z. Ni, A. Millward, A.J.H. McGaughey, C. Uher, M. Kaviany, and O. Yaghi, "Thermal conductivity of a metal-organic framework (MOF-5): Part II. Measurement", Int. J. Heat Mass Transf. 50, 405 (2007).
A. Chakraborty and S. Kumar, "Thermal management and desorption modeling of a cryo-adsorbent hydrogen storage system", Int. J. Hydrog. Energy 38, 3973 (2013).