Multispecies Bhatnagar-Gross-Krook models and the Onsager reciprocal relations

本論文は、多成分系に対する既存の速度非依存型ブハナガー・グロス・クロックモデルの大半がオンサーガーの相反関係を満たさないことを示しており、それによって特定の流体に一致するように輸送特性を較正することが複雑化している。

要約

2 種類の異なるビー玉（赤と青としましょう）を箱の中で揺らしたときにどのように混ざり合うかを予測しようとしていると想像してください。科学者たちは、これらのビー玉の振る舞いに関する「規則集」を持っており、それをボルツマン方程式と呼びます。これは驚くほど正確ですが、あまりにも複雑で、これを解くことは、ハリケーンが吹き荒れる間に海岸の砂粒を一粒ずつ数えようとするようなものです。

これを容易にするため、科学者たちはこの規則集の簡略化されたバージョンである BGK モデル（バトナガル、グロス、クロークにちなんで命名）を作成しました。BGK モデルは、すべての重たい数学計算を行うことなく、ビー玉の複雑な振る舞いを近似する「カンニングペーパー」や「近道」と考えてください。これは、翼を流れる空気から恒星内のプラズマに至るまで、あらゆるものをシミュレーションするために数十年にわたって使用されてきました。

問題：壊れたコンパス

E. S. ベニロフによって書かれたこの論文は、異なる気体の混合物（酸素と窒素、または水蒸気と空気など）を扱う際、この「カンニングペーパー」の最も一般的なバージョンに重大な欠陥があることを指摘しています。

著者は、これらの人気のある BGK モデルが、オンサーガーの相反関係として知られる物理学の根本法則に違反していることを発見しました。

オンサーガーの関係についての簡単な比喩：
2 つの町の間を交通が往来する一方通行ではない道路を想像してください。

• 規則 A: 町 A に丘を作ると、町 B での車の走行速度に影響を与えます。
• 規則 B: 町 B に丘を作ると、町 A での車の走行速度に影響を与えます。

オンサーガーの関係は、これら 2 つの影響が互いに完全に「調和」していなければならないと述べています。町 A の丘が町 B の交通を 10% 遅くさせるなら、町 B の丘は数学的に関連付けられた量だけ町 A の交通を遅くしなければなりません。これは対称性の規則であり、宇宙はこれらの相互作用が釣り合うことを要求します。

論文が見つけたこと

ベニロフは、この規則に対して標準的な BGK「カンニングペーパー」をテストしました。その結果、以下のことがわかりました。

1. モデルは非対称である: BGK モデルでは、町 A の「丘」（温度変化）は、町 B の「交通」（質量流れ）にゼロの影響しか与えません。しかし、町 B の「丘」（密度変化）は、町 A の「交通」（熱流）に影響を及ぼします。
2. 不一致: 方程式の一方がゼロで他方がゼロでないため、対称性が破れています。このモデルは、片側に永久に傾いた秤のようです。
3. 結果: モデルがこの根本的な規則を破るため、それを「較正」することは不可能です。較正とは、明確な信号を得るためにラジオを調整するようなものです。特定の流体の実世界データと一致させるために BGK モデルのノブ（パラメータ）を調整しようとしても、できません。モデルは根本的に破損しており、どのように微調整しても決して完全に正確になることはありません。

「水蒸気」の例外（そしてなぜそれが救いにならないのか）

「まあ、これは奇妙な気体だけにしか関係ないのかもしれない。水蒸気と空気のような一般的な混合物はどうだろう？」と思うかもしれません。

論文もこれを検証しています。温度が質量流れに与える影響が微小であっても（水蒸気と空気の場合、実際そうである）、モデルは依然として失敗します。この特定のケースでモデルを機能させるには、ダイヤルを無限大に回さなければならず、それは実質的にすべての運動を停止させることでモデル全体を壊してしまいます。したがって、このモデルは複雑な混合物だけでなく、単純な混合物に対しても失敗します。

優れたモデルは存在するか？

論文は、規則に従ういくつかのより複雑な BGK モデルが存在することに言及していますが、それらには独自の問題（エントロピーが常に増加することを保証する「H 定理」など、他の物理法則を破るなど）があります。

著者は結論として、現在、完璧な BGK モデルは存在しないと述べています。完璧なモデルには以下が必要でしょう。

• 質量、運動量、エネルギーを保存すること。
• 熱力学の法則（エントロピー）に従うこと。
• 同一の粒子を公平に扱うこと。
• 温度と濃度を正に保つこと。
• 科学者があらゆる実在の流体と一致するように調整できるようにすること。
• そして、オンサーガーの相反関係（対称性の規則）に従うこと。

現在、私たちが持っているすべてのモデルは、これらのテストの少なくとも 1 つに失敗しています。

結論

この論文は、これらのモデルを使用する科学者への警告です。気体混合物をシミュレートするために標準的な BGK モデルを使用している場合、物理法則と根本的に「調律外れ」のツールを使用していることになります。それは何が起こっているかの大まかなアイデアを与えるかもしれませんが、自然の核心的な対称性規則に違反しているため、実在の流体に対する正確で較正された結果を提供するとは信頼できません。著者は、将来、これらの問題のすべてを修正する「完璧な」モデルが誰かによって構築されることを願っています。

技術概要

技術的概要：多種類 Bhatnagar–Gross–Krook モデルとオンサーガーの相反定理

問題提起
Bhatnagar–Gross–Krook（BGK）モデルは、特に気体混合物に対して広く用いられるボルツマン運動論方程式の現象論的近似である。元の単一種類モデルおよび混合物への拡張（例：Gross および Krook [2]）は、その定性的な妥当性と適応性から評価されているが、オンサーガーの相反定理への適合性については厳密に検証されていない。これらの関係は輸送係数に根本的な熱力学的制約を課し、質量流における温度勾配の係数（熱拡散またはソレー効果）と、熱流における密度勾配の係数（デュフォア効果）を相互に結びつける。第一原理から導出されたモデルはこれらの関係を自動的に満たすが、現象論的モデルは満たさない可能性がある。これらの関係への不適合は、モデルの物理的妥当性に対する疑念を招くだけでなく、決定的な点として、特定の実際の流体の輸送特性に一致させるためのモデルパラメータの較正を不可能にする。

手法
本論文は、単原子気体の二成分混合物に対する最も一般的な多種類 BGK モデルを定式化する。この定式化には、種に依存する速度と温度を持つ混合マクスウェル分布（$M_{ij}$）への緩和を定義する衝突頻度（$\nu_{ij}$）およびパラメータ（$\alpha_i, \beta_i, \gamma_i$）が含まれる。著者は、これらのパラメータに特定の制約を課すことで、モデルが質量、運動量、およびエネルギーを保存することを保証する。

輸送特性を分析するために、本論文は標準的な流体力学的近似ではなく、拡散近似を採用する。これには以下の手順が含まれる：

1. スケーリング：遅く、弱い拡散流に適した微小パラメータ $\epsilon$ を導入し、空間座標（$\nabla \to \epsilon \nabla$）および時間（$\partial_t \to \epsilon^2 \partial_t$）を伸長する。
2. 展開：分布関数、巨視的速度、温度、および密度を $\epsilon$ のべき級数として展開する。
3. 導出：これらの展開を再スケーリングされた BGK 方程式に代入し、分布関数の一次補正（$f^{(1)}_i$）を解く。
4. 流束計算：$f^{(1)}_i$ を積分し、密度勾配および温度勾配を用いた質量流束（$J_i$）および熱流束（$Q$）の明示的な式を導出する。
5. 比較：導出された BGK 流束を、オンサーガー関係を満たす標準的な流体力学的式（Enskog–Chapman 理論）と比較する。

主要な貢献と結果
本論文の中心的な結果は、衝突頻度およびモデル係数が分子速度に依存しない、最も一般的な多種類 BGK モデルのバージョンが、オンサーガーの相反定理を満たさないことを実証した点である。

具体的には、分析は根本的な矛盾を明らかにする：

• オンサーガーの関係によれば、温度勾配を質量流束に結びつける係数（熱拡散率、$B$）は、密度勾配を熱流束に結びつける係数（デュフォア効果、$C$）が非ゼロである場合、非ゼロでなければならず、かつ $C = pB$ によって結びつけられる。
• BGK モデルの導出は、特定の非物理的な極限を除き、熱拡散項が厳密にゼロ（$B=0$）となる質量流束をもたらす。
• しかし、得られる熱流束には、密度勾配に比例する項（デュフォア効果）が含まれており、これは非ゼロであり、かつ種の質量差に比例する。
• 結果として、BGK モデルはゼロのソレー効果と並んで非ゼロのデュフォア効果を予測し、オンサーガー関係 $C = pB$ に違反する。

さらに、本論文は、この不一致がモデルパラメータ（相対速度緩和を支配するパラメータ $\beta$ など）の較正によって解決できないことを示している。一般的な混合物において、標準的な輸送方程式とオンサーガー関係を同時に満たす $\beta$ の値は存在しない。熱拡散率が無視できる流体（$B \approx 0$）の特定のケースであっても、BGK モデルは、拡散流束全体をゼロに強制しない限り、正しい熱流束の挙動を再現できない。

意義と主張
本論文は、この不適合が多種類 BGK モデルの実用的応用に対して重大な問題をもたらすと主張している。モデルから導出された輸送係数が必要な熱力学的相反性を遵守しないため、モデルを「較正」すること、すなわちモデルの輸送特性を特定の実際の流体（例えば水蒸気と空気など）のそれらに一致させるパラメータ値を選択することは不可能である。

著者は、モデルがオンサーガー関係を厳密には満たさないものの、以下の 2 つの特定の極限において漸近的に満たすことを指摘している：

1. 質量比 $m_1/m_2 \to 0$ の場合（電子質量がイオンに比べて無視できる電離プラズマに関連する）。
2. 質量比 $m_1/m_2 \to 1$ の場合（熱流束式における第一項が消滅する）。

本論文は最後に、「完璧な」多種類運動論モデルの要件を列挙して結論付けている。これらには、保存則、H 定理、無差別性の原理、温度・濃度の正値性、任意のプラントル数・シュミット数を表現する能力、およびオンサーガー関係への適合が含まれる。著者は、現在までに、これらすべての要件を同時に満たすことが示された BGK 型モデルは存在しないと述べている。これらの知見は、現在の BGK モデルは有用な近似ではあるものの、オンサーガー関係を満たさない能力が、その物理的忠実度および一般流体力学応用における精密な較正への有用性を制限していることを示唆している。