A型ゼオライト吸着筒のマルチスケール熱流体シミュレーション: 3D-RISM-SCF計算の応用#

Advance/PHASE Logo

デシカント空調や調湿建材などの高効率な熱エネルギー・空調システムにおいて、A型ゼオライトに代表される多孔質吸着材が重要な役割を担っています。しかし、ゼオライトが水分を吸着する際に発生する「吸着熱」によって材料温度が上昇し、マクロな吸着性能が著しく低下するという工学的な課題があります。本事例では、第一原理計算ソフトウェアAdvance/PHASEに実装された3D-RISM-SCF法によって得られたミクロな「飽和吸着量」と「吸着エネルギー」を、マクロな1次元熱流体・物質収支モデルへ直接橋渡しするマルチスケール連携手法を紹介します。これにより、実機サイズの吸着筒内部で生じる熱波（Heat wave）と物質波（Mass wave）の分離現象をシミュレーションで検証することが可能となりました。

Keywords: 第一原理計算 (DFT), 3D-RISM-SCF, A型ゼオライト, 水分吸着, マルチスケール連携, 1次元充填層モデル, 破過曲線

解析モデルと理論背景#

1. 3D-RISM-SCF法を用いた第一原理計算：微視的な吸着特性の評価#

解析の第一段階として、Advance/PHASEに実装された3D-RISM-SCF法（統計力学的溶液論と第一原理計算の連成手法）を用い、A型ゼオライト（Na-A型）の空孔内における水分子の吸蔵状態をシミュレーションしました。連続誘電体モデルでは取り扱いが困難な「溶媒分子の明示的な配置と相関」を高精度に計算しています [1]。

解析の結果、ユニットセル内の水分子数が27.5個に達した時点で相対エネルギーが極小値（安定点）を迎えることが確認されました。この結果から、マクロな理論飽和吸着量 $q_\text{max} = 0.291$ kg/kg が導出されます。また、水分子数に対するエネルギー曲線の勾配（N=20〜25付近の平均で約-1.7 kcal/mol）は、「液体の水」を基準状態とする3D-RISM法における吸着エネルギー（約390 kJ/kg）を表しています。空調モデルで扱う「水蒸気（気相）」からの吸着熱を求めるため、これに水の蒸発潜熱（約2442 kJ/kg）を加算し、最終的な微分吸着熱 $\Delta H_\text{ads} \approx 2832$ kJ/kg を算出しました。なお、文献（Gorbachら [2]）で報告されているゼオライト4Aの等量吸着熱（約3050 kJ/kg）と比較しても、第一原理計算に基づく本予測値が妥当な物理スケールを示していることが分かります。これらの値が、後述のマクロ解析における「容量」と「熱源」の主要パラメータとなります。

2. 線形推進力（LDF）モデル：ミクロとマクロの橋渡し#

ゼオライト細孔内への水分子の拡散・吸着速度をマクロシミュレーションに組み込むため、Glueckaufらによって提唱された線形推進力（Linear Driving Force: LDF）モデルを採用しました [3]。このモデルは、複雑な細孔内拡散をシンプルな一次遅れ系で近似する実用的な手法です。

$\frac{\partial q}{\partial t} = k_\text{LDF} (q^* - q)$

ここで、 $q$ は現在の吸着量、 $q^*$ は空気の湿度に依存する平衡吸着量、 $k_\text{LDF}$ は総括物質移動係数です。本解析では、低湿度域で急激に飽和に達するType-Iの吸着等温線の挙動を模し、第一原理計算から得られた最大吸着量 $q_\text{max}$ を上限とする、単純化したLangmuir型等温線を採用しています。また、水分拡散の遅れを表現する総括物質移動係数 $k_\text{LDF}$ を一定として組み込みました。これにより、ミクロな限界容量をマクロな速度論的挙動に反映させています。

3. 1次元熱物質収支モデル（1D Packed Bed Model）#

吸着筒（充填層）全体を通した空気の温湿度変化をシミュレートするため、Ruthvenらの古典的な吸着工学理論 [4] に基づく1次元熱・物質収支方程式を構築しました。空気の流れ（移流）に伴う質量保存則と、吸着熱および壁面からの放熱を考慮したエネルギー保存則を連立して解きます。

$\frac{\partial Y}{\partial t} + v \frac{\partial Y}{\partial x} = - \frac{1-\epsilon}{\epsilon} \frac{\rho_s}{\rho_g} \frac{\partial q}{\partial t}$

$\rho_\text{eff} C_{p,\text{eff}} \frac{\partial T}{\partial t} + \epsilon \rho_g C_{p,g} v \frac{\partial T}{\partial x} = (1-\epsilon) \rho_s \Delta H_\text{ads} \frac{\partial q}{\partial t} - \frac{4 h}{D_\text{col}} (T - T_\text{env})$

ここで、 $Y$ は空気の絶対湿度、 $T$ は温度、 $v$ は間隙流速（Interstitial velocity）、 $\epsilon$ は空隙率です。エネルギー収支式の右辺第2項は壁面からの熱伝達（放熱）を表し、 $h$ は総括熱伝達係数、 $D_\text{col}$ はカラム直径、 $T_\text{env}$ は環境温度です（断熱条件の場合は $h=0$ とします）。右辺のソース項に含まれる吸着熱 $\Delta H_\text{ads}$ と吸着速度 $\partial q / \partial t$ に、3D-RISM-SCF計算の結果がダイレクトに反映されます。図1に、本マルチスケール連携の全体構成となる模式図を示します。

1次元充填層モデルとミクロ連携の模式図
図1. マルチスケール連携による1次元充填層吸着モデルの模式図。3D-RISM-SCF計算から得られたミクロパラメータ（最大吸着量と吸着熱）は、マクロなカラム内の熱・物質収支式にソース項として入力されます。内部の波形は、熱波（赤）と物質波（青）の進行速度の差を概念的に表しています。

4. 吸着筒の仕様と運転条件#

デシカント空調のロータや吸着筒を想定し、長さ0.2 mの充填層に対して25℃の湿潤空気を流し込む条件でシミュレーションを実施しました。

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

カテゴリー	パラメータ	設定値と物理的根拠（出典）
ミクロ連携入力 (3D-RISM-SCF)	最大吸着量 ( $q_\text{max}$ )	0.291 kg/kg (安定点27.5個から換算) [1]
	吸着熱 ( $\Delta H_\text{ads}$ )	2832 kJ/kg (エネルギー勾配などより算出) [1]
マクロ物性値	ゼオライトかさ密度 ( $\rho_s$ )	720 kg/m³ (市販A型ゼオライトの代表値) [2]
	空隙率 ( $\epsilon$ )	0.4 (充填層の一般的な空隙率) [2]
	LDF物質移動係数 ( $k_\text{LDF}$ )	0.005 1/s (文献測定値[2]における濃度依存性を考慮した外挿値)
装置・境界条件	充填層長さ ( $L$ )	0.2 m
	空気間隙流速 ( $v$ )	0.5 m/s
	入口空気状態 ( $Y_\text{in}, T_\text{in}$ )	15.0 g/kg, 25.0 °C

解析結果と考察#

1. 破過曲線のシミュレーションと熱波の先行現象（断熱条件）#

図2は、断熱条件で構築した1次元熱物質収支モデルにより計算されたカラム出口における空気の湿度（青線）と温度（赤線）の時間変化（破過曲線：Breakthrough curve）です。解析開始直後、ゼオライトが水分を急速に吸着することで第一原理計算から与えられた大きな「吸着熱」が放出されます。この熱は空気の流れに乗って素早く下流へ運ばれるため、湿度が上昇するよりもはるか手前の段階（約300秒付近）で、出口温度が70℃近くまで急上昇する「熱波（Heat wave）」の到達が確認できます。一方、水分子は吸着層内にトラップされながら進み、A型ゼオライトの強力な吸着力によって細孔内の吸着サイトが飽和するまで捕捉され続けるため、出口の湿度は長時間ゼロに保たれます。結果として、出口湿度が元の15.0 g/kgに戻る「物質波（Mass wave）」の到達（破過）には長時間を要します。この断熱層特有の熱波と物質波の完全な分離現象は、Panらの古典的理論 [5] と整合するものです。

さらにグラフ後半（2500秒以降）の挙動に注目すると、出口湿度が上昇し始める（破過が始まる）のと同期して、出口温度が徐々に低下しています。これは、ゼオライトが水分を吸いきって飽和状態に近づくにつれて「新たな吸着熱が発生しなくなる」ためです。熱源を失った筒内は、入口から供給され続ける25℃の空気によって冷却されていきます。すなわち、「熱波の到達」は物質波よりも圧倒的に早いものの、「熱波の減衰」は物質の飽和とほぼ完全に同期して起こるという、物理的に妥当な現象が本シミュレーションによって捉えられています。

なお、本シミュレーションで示された約70℃という高温の熱波（温度スパイク）は、文献結果と比較することができます。例えば、Panらの断熱層における水分吸着実験 [5] では、アルミナ吸着材を用いた場合でも入口温度約38℃（100°F）、入口絶対湿度約19.4 g/kg（モル比 $Y = 0.0312$ ）の条件で約106℃（223°F）の熱波到達温度が実測されており、より低湿度の条件でも60〜70℃を超える高温が報告されています。第一原理計算に基づく本計算結果（約70℃への到達）は、物理的に妥当なスケール感であると言えます。

破過曲線と温度変化の初期シミュレーション結果破過曲線と温度変化の長時間シミュレーション結果
図2. 出口における空気絶対湿度（青）と温度（赤）の初期シミュレーション（上）、および長時間シミュレーション（下）の結果。第一原理計算に基づく吸着熱により、出口温度が一時的に約70℃までスパイクするという、断熱層特有の熱波の先行現象が再現されています。下図の後半で湿度の上がり始めと温度の下がり始めが同期しているのは、ゼオライトが飽和して吸着熱の発生が停止したためです。※グラフ上部における線分の交差は、左右の縦軸スケールの違いによる見かけ上のものです。

2. 非断熱条件（自然放熱）における熱波と物質波の挙動#

次に、実際の吸着筒（金属やプラスチック製のカートリッジなど）が室内に置かれ、周囲の空気へ自然に放熱（自然対流）しているケースを想定し、非断熱モデルでのシミュレーションを行いました。総括熱伝達係数 $h = 10 \text{ W/(m}^2\text{K)}$ （自然対流）、カラム直径 $D_\text{col} = 0.01$ m を設定してシミュレーションを実施しました。

自然放熱条件における破過曲線と温度変化
図3. 自然放熱条件（ $h = 10 \text{ W/(m}^2\text{K)}$ ）における出口空気絶対湿度（青）と温度（赤）のシミュレーション結果。断熱条件（図2）と比較して、ピーク温度の低下や波形のブロード化が確認できます。

図3のシミュレーション結果を断熱条件と比較すると、熱流体・吸着工学的に興味深い3つの特徴が現れています。

ピーク温度の現実的な低下： 壁面からの放熱により、断熱条件で約70℃に達していた温度スパイクが約41℃にまで抑え込まれています。これは、「吸着に伴ってカートリッジに触れると少し熱く感じる」という現実のデバイスでよく見られる妥当な温度スケールです。
熱波と物質波の「中間の分離（ブロード化）」： 断熱条件では熱波が物質波より完全に先行して分離していましたが、自然放熱条件では赤い温度のピーク（約5200秒付近）が青い湿度の立ち上がりの中心（約6600秒付近）よりも先行しています。「熱が前方に運ばれる力」と「壁から外へ逃げる力」が拮抗し、熱波が前方にブロードに広がった様子が表現されています。
初期の過渡的な挙動（肩の形成）： グラフ左側の250秒付近に、温度上昇の小さな段差（肩）が見られます。これは計算開始直後にカラム入口付近で爆発的に吸着熱が発生し、それがまだ壁面からの放熱とバランスする前の過渡的な挙動を示しています。

3. 材料設計とシステム設計の連携価値#

本結果は、単に「材料の空孔が何個の水分子を吸う」というミクロな事実だけでなく、「実際の空調システムに組み込んだ際、断熱条件で除湿された空気が70℃の熱風になって排出される」といったマクロなシステム上の挙動や課題を、試作を行う前にシミュレーションによって浮き彫りにできることを意味します。また、自然放熱を考慮したモデルを活用することで、「放熱によって温度上昇が現実的な41℃程度に緩和され、熱波がブロードに広がる」といった、実際のデバイスの設置環境により即した評価も可能となります。

もしこの出口温度が高すぎる場合、設計者はマクロな観点から「システムに冷却機構や放熱フィンを追加する」か、あるいはミクロな観点に戻り「3D-RISM-SCF計算を活用して、骨格金属を置換することで吸着熱（エネルギー曲線の勾配）がより穏やかな新しいゼオライトを探索する」といった具体的な対策を打つことが可能になります。

本解析アプローチの特徴と適用範囲#

吸着等温線の取り扱い：
本解析では第一原理計算から得られた最大吸着量 $q_\text{max}$ を直接組み込むため、Langmuir型等温線を単純化したモデルを採用しています。ただし、Ruthven [4] も指摘するように、ゼオライトのような微細孔（マイクロポア）を持つ吸着材は、細孔が水分子で完全に満たされる明確な飽和限界を持つため、低湿度域で急激に立ち上がる非線形な「Type-I」の吸着等温線を示します。より高度なシミュレーションを行うには、第一原理計算や実験データからGorbachら [2] のような詳細な等温線モデルを導出し、それをLDFモデルの平衡吸着量 ( $q^*$ ) の計算に組み込むことが望ましいです。
※本解析のモデルに実験データに基づく非線形な「Type-I」の吸着等温線 [2] を取り入れてシミュレーションした結果は、今回の単純化したLangmuir型等温線モデルでの結果に近いことが確認されています。
物理モデルの簡略化：
計算コストと実用性を重視し、空間1次元モデルを採用しています。円筒ロータの回転や、ハニカム構造壁面からの詳細な放熱、3次元的な偏流効果などを厳密に評価する必要がある場合は、より汎用的な3D熱流体ソルバーの使用が求められます。
材料スクリーニングへの応用：
本手法は、合成や実験による測定が面倒な「吸着熱」や「飽和吸着量」といった重要パラメータを、第一原理計算によって代替（予測）できる点にあります。密度や比熱などのマクロな基本物性に類似材料の仮定値を置くことで、未知の新材料であっても合成前にデバイス性能のポテンシャル（熱波の影響度など）を概算できます。ここで得られる概算値は、マテリアルズ・インフォマティクス（MI）と組み合わせたハイスループットな材料設計の有力な評価指標となります。

より複雑・詳細な熱流体シミュレーションへ
本事例では、物理現象の理解やマルチスケールの概念実証（Proof of Concept）として、Pythonソルバーによる1次元簡易モデルを適用しました。実際の複雑な製品形状における3次元的な気流解析や詳細設計においては、Advance/FrontFlow/red などの流体解析ソフトウェアとの連成解析が推奨されます。

まとめ#

本事例では、Advance/PHASEの3D-RISM-SCF法による第一原理計算から得られた「最大吸着量」と「吸着熱」を、マクロな1次元熱流体・物質収支方程式へ組み込むマルチスケール連携を実施しました。その結果、断熱条件における熱波と物質波の完全な分離や、非断熱条件における自然放熱の挙動など、ゼオライト吸着筒におけるマクロな出口温湿度の変化を、シミュレーションによって算出可能であることを実証しました。このアプローチは、原子・分子レベルの材料設計が、最終製品としての空調効率や省エネ性能にどう貢献するのかを定量的に評価する「デジタルツイン」構築の第一歩となります。

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

お問い合わせ

参考文献#

3D-RISM-SCFを用いたA型ゼオライトへの第一原理計算の応用
A. Gorbach, M. Stegmaier, and G. Eigenberger, "Measurement and Modeling of Water Vapor Adsorption on Zeolite 4A—Equilibria and Kinetics", Adsorption 10, 29 (2004).
E. Glueckauf, "Theory of chromatography. Part 10.—Formulae for diffusion into spheres and their application to chromatography", Trans. Faraday Soc. 51, 1540 (1955).
D. M. Ruthven, Principles of Adsorption and Adsorption Processes, John Wiley and Sons (1984).
C. Y. Pan and D. Basmadjian, "An analysis of adiabatic sorption of single solutes in fixed beds: pure thermal wave formation and its practical implications", Chem. Eng. Sci. 25, 1653 (1970).