[論文レビュー] Numerical analysis of the thermal relaxation of the dense gas between two parallel plates: the free energy monotonicity for the Enskog equation

Shigeru Takata, Soma Sakata+2 more|arXiv (Cornell University)|Mar 25, 2026

Gas Dynamics and Kinetic Theory被引用数 0

ひとこと要約

論文は Enskog 方程式を元の Enskog 因子（OEE）と修正因子（EESM）と比較し、並平板間の高密度ガスを対象に、EESM では自由エネルギーが単調減少することを示す一方、OEE はそうならない。

ABSTRACT

The thermal relaxation problem between two parallel plates with the same temperature is investigated, aiming to study the behavior of the free energy of the dense gas described by the Enskog equation. Two types of Enskog equation have been used: one is the Enskog equation with the original Enskog factor, while the other is that with a modified Enskog factor proposed recently in Takata & Takahashi, Phys. Rev. E 111, 065108 (2025). The evaluated free energy is a natural extension of the thermodynamic free energy to the non-equilibrium state. It is observed that this free energy monotonically decreases in time for the modified factor version, while it is not necessarily the case for the original version. Differences are also observed in other quantities in their time evolutions, most typically in the density profile.

研究の動機と目的

• dense-gas シミュレーションにおける熱力学的一貫性を動機づけるため、Enskog 型運動モデルに基づく自由エネルギーの時間発展を検討する。
• 系が熱浴と接触している場合に、拡張エントロピから定義された自由エネルギーが単調減少するかを調べる。
• 原始の Enskog 方程式（OEE）と微小改良を加えた Enskog 方程式（EESM）を自由エネルギー挙動と密度プロファイルの観点から比較する。

提案手法

• Enskog 衝突積分の数値解法として有限差分時間進行法と高速フーリエスペクトル法を組み合わせたハイブリッド数値スキームを用いる。
• 二つの Enskog 変種を実装する：原始 Enskog 因子（OEE）と、Carnahan–Starling エンタルピーを用いた修正 Enskog 因子（EESM）。
• 等温平板間の拡散反射境界条件を持つ一次元の熱緩和問題を解く。
• 拡張エントロピーに基づく自由エネルギー汎関数を定義し、それの時間的単調性を評価する。
• 速度空間を周期的領域で離散化し、時間進行には Adams–Bashforth 外挿法を用い、数値質量を正規化する。

Figure 1: Problem setting. Dashed circle indicates a molecule in contact with the plate surface.

実験結果

リサーチクエスチョン

• RQ1 熱浴と接触している場合、修正 Enskog 因子（EESM）版で自由エネルギーは単調に減少するか？
• RQ2 同じ設定で OEE は自由エネルギーの単調減少を保証するか、進行中の密度プロファイルはどのように変わるか？
• RQ3 過渡緩和中に OEE と EESM の密度・温度プロファイルはどのように異なるか？
• RQ4 分子サイズと初期密度の変動が自由エネルギーおよび他の巨視的場の時間発展に与える影響は何か？

主な発見

• 修正因子版（EESM）では自由エネルギーが時間的に単調に減少する。
• 原始 Enskog 方程式（OEE）では自由エネルギーが必ずしも単調ではない。
• 特に初期密度変動が大きい場合、EESM と OEE の過渡的な密度プロファイルの違いが現れる。
• 最終的な定常状態の密度は非一様である一方、温度は両モデルで均一となり、プレート温度に等しくなる。
• 自由エネルギーの理想気体成分は単調に変化せず、EESM における単調性には運動学的成分以外が重要であることを示す。
• 初期緩和の初期段階では EESM と OEE の密度プロファイルに顕著な過渡的差異が見られるが、最終状態は類似している。

Figure 2: The Enskog factor ${\mathsf{g}}(0,\Delta X_{1})$ around the midpoint between the plates in the initial state in the case of $\eta_{0}=0.25$ , $\sigma/L=0.1$ ( ${\mathrm{Kn}}=0.0227$ ), and $\lambda=0.1$ for $w=0$ , $0.3$ , $0.5$ , and $0.7$ . Solid lines: EESM. Dashed (red) lines: OEE.