Compressible and immiscible fluids with arbitrary density ratio

本論文は、圧縮性かつ非混和性の流体を任意に高い密度比でモデリングするための新しい熱力学に基づく理論を提示するものであり、真の運動量進化を適切に考慮した密度進化方程式を導出することによって、従来のナビエ・ストークス方程式およびオイラー方程式の限界に対処するものである。

要約

二人の非常に異なるパートナーによるダンスを想像してみてください。それは、重くゆっくりと動く象（水）と、軽快に素早く動く羽（空気）のダンスです。流体力学の世界では、これらは「混ざり合わない（immiscible）」、つまり油と水のように混ざり合うことはありませんが、両者が出会う場所には「ぼやけた境界」が存在します。

長い間、科学者たちは、密度の差が極めて大きい（例えば、象は羽よりも1,000倍重い）系における「ダンスのルール（方程式）」を記述することに苦心してきました。従来のルールには大きな欠陥がありました。それは、象の重さを、たとえ象が羽へと変化していくあの「ぼやけた境界」においてさえも、変化しない一定の数値として扱っていたことです。これは、人が幽霊へと変わっていく様子を描写する際に、「彼はずっと固体の人間のままである」と言うようなもので、全く理にかなっていません。

以下に、この論文がその問題をどのように解決しようとしているのか、簡潔な内訳を説明します。

1. 旧来のルールの問題点

流体の動きを計算する伝統的な手法（ナビエ・ストークス方程式およびオイラー方程式を用いる方法）は、**ブジネスク近似（Boussinesq approximation）**と呼ばれるショートカットに依存しています。

• 比喩： あなたがショッピングカートを押しているところを想像してください。カートにレンガが詰まっていれば重いです。空っぽであれば軽いです。旧来のルールは、もしあなたが「レンガが半分入っていて、半分は空気が入っているカート」を押しているとしても、押している最中にカートの重さは決して変わらない、と仮定していました。つまり、重さは固定された平均値であると想定しているのです。
• 欠陥： 現実には、カートが動き、中のレンガが移動（拡散）するにつれて、重さは変化します。旧来のルールは、境界において質量そのものが変化しているために「運動量（質量 × 速度）」が変化するという事実を無視していました。また、空気と水が占める空間の量は完全に予測可能であると仮定していましたが、それらが境界で混ざり合う際、圧力によって占有する空間が実際に揺らいだり変化したりすることを無視していました。

2. 新しいアプローチ：エネルギー最小化

著者は、密度の変化を推測する代わりに、ある根本的な原理から出発しています。それは、**「自然界は常に、可能な限り最小のエネルギーを使おうとする」**という原理です。

• 比喩： 丘を転がり落ちるボールを考えてみください。ボールはどのようなステップを踏むかには関心がありません。ただ、最も低いエネルギー状態（丘の底）に行きたがっているのです。著者は、この「エネルギーの丘」という概念を用いて、水と空気がどのように相互作用するかについての新しいルールを導き出しています。
• 核心となる革新： 著者は**「過剰体積（Excess Volume）」**という概念を導入しています。
  • 水のバケツと空気のバケツがあるとします。これらを一緒に注いだとき、合計の体積は正確に二つのバケツの和になると予想するかもしれません。しかし、微視的なレベルでは、それらが合流する際、分子がより密に、あるいはより疎に集まることがあり、それによって「余分な」あるいは「足りない」空間が生じます。
  • 旧来のルールは、この余分な空間はどこでもゼロであると仮定していました。しかし本論文は、「いいえ、この余分な空間は存在し、場所ごとに変化し、それが密度に影響を与えるのです」と述べています。

3. 新しい「ダンスのルール」（結果）

この「余剰体積」の変化を考慮し、エネルギー最小化を用いることで、著者は主に3つのことを行う新しい方程式のセットを導き出しました。

A. 音の聞こえ方の新しい方法（音速）
本論文は、音速は単なるランダムな数字ではなく、流体が押しつぶされる際にそのエネルギーがどのように変化するかから直接導かれるものであることを示しています。

• 比喩： 音を、群衆の中を伝わる波紋と考えてみてください。その波紋の速度は、人々（分子）がいかに密集しているか、そして彼らがどれほどのエネルギーを持っているかに依存します。新しい公式は、事前に音速を教えられることなく、この速度を自然に計算します。さらに、気体においては、音速はガス分子が跳ね回る平均速度とほぼ等しくなることも示唆しています。

B. 圧力と速度に関する新しいルール（ベルヌーイの法則）
ベルヌーイの原理として、「流体の速度が上がると、圧力は下がる」という言葉を聞いたことがあるでしょう。

• ひねり： 旧来のルールは、管の中を流れる水には非常にうまく機能しますが、密度の劇的な変化（水が空気に当たる場合など）があるときには破綻します。著者は、**「一般化されたベルヌーイの法則」**を作成しました。
• 比喩： 川が滝へと流れ込む様子を想像してください。旧来のルールは、エネルギーは一定であると言います。しかし新しいルールは、「待ってください。水が霧（空気）に変わる際、水の『詰まり具合』が変化するため、一部のエネルギーが失われたり変換されたりするのです」と指摘します。新しい方程式はこのエネルギーの変化を考慮しており、流体の性質が重いものから軽いものへと変化する場合でも正確に機能します。

C. 密度の「隆起（バンプ）」
これはおそらく、最も視覚的な結果です。

• 旧来の視点： 水と空気の境界を見るとき、旧来のモデルでは、密度は「重い水」から「軽い空気」へと、まるで坂道を下るように滑らかに減少するとされていました。
• 新しい視点： 著者の数学的モデルは、**「隆起（バンプ）」**を予測します。境界を横切る際、密度は空気のレベルまで下がる前に、一度わずかに上昇します。
• 比喩： 人々の群れ（水）が、空の廊下（空気）へ向かうドアを通り抜けようとしている場面を想像してください。ドアを通り抜ける瞬間、彼らは広がる前に、一瞬だけぎゅっと密集するかもしれません。新しい理論はこの「パッキング・バンプ（密集による隆起）」を予測しており、これは高度なコンピュータ・シミュレーション（密度汎関数理論）がすでに示している現象とも一致していますが、従来の単純なモデルでは見落とされていたものです。

まとめ

この論文は、重さが大きく異なる流体（水と空気のような）のための物理法則を記述する、新しい方法を提案しています。

1. 境界において重さが一定であるという仮定を止めました。
2. 分子が占める空間が変化すること（過剰体積）を認めました。
3. 「最小エネルギー」の原理を用いて、音の伝わり方、圧力の変化、そして水と空気の境界において密度が実際に小さな「隆起」を持つ理由を説明する新しいルールを導き出しました。

著者は、この新しいフレームワークが、重さの異なるあらゆる流体の混合に対して機能すると主張しており、雨が降る様子、波が砕ける様子、あるいは泡が上昇する様子といった現象の、より正確なコンピュータ・シミュレーションへの道を開くものです。

技術概要

技術要約：任意の密度比を持つ圧縮性および非混和性流体：拡散界面理論

問題提起
現在の数値流体力学（CFD）モデルは、水と空気（比率 $\approx 1000$）のような高密度比を持つ非混和流体系の正確なシミュレーションにおいて困難に直面している。既存のアプローチは主に2つの手法に依存しているが、これらには重大な限界が存在する。

1. ブシネスク近似（Boussinesq Approximation）： 界面内で密度（$\rho$）が一定であると仮定し、$\rho(du/dt) = -\nabla p + \dots$ という方程式を用いる。この定式化は、真の運動量進化 $d(\rho u)/dt$ を考慮しておらず、密度の変化とそれに伴う運動量転送を無視している。さらに、密度差から自然に生じる重力や浮力を、アドホックな追加なしに生成することはできない。
2. 補間法（Interpolation Method）： 密度を体積含有率の線形関数（$\rho = \rho_w \phi + \rho_a(1-\phi)$）として定義する。この手法は、非圧縮性（$d\rho/dt = 0$）と拡散（$d\phi/dt \neq 0$）の間に根本的な矛盾を生じさせる。また、ドメイン全体で過剰体積が一定（多くの場合ゼロ）であると仮定しているが、これは局所的な過剰体積が局所的な圧力に依存するという物理的事実に反している。その結果、標準的な連続の式や線形補間では、水と空気の界面で見られる非単調な密度プロファイルを捉えることができない。

手法
著者らは、従来のナビエ・ストークス方程式やオイラー方程式の導出とは異なる、エネルギー最小化の熱力学的原理に基づいた新しい理論的枠組みを提案している。その手法は以下の通りである。

• 体積と密度の再定式化： 過剰体積をゼロと仮定する代わりに、本モデルでは界面内に空間的に変化する過剰体積（$v_e$）を導入する。この体積は流体成分間で再分配され、その結果、部分密度（$\rho'_w, \rho'_a$）は一定ではなく、局所的な圧力と混合に依存するものとなる。これには、体積含有率（$\phi$）と密度（$\rho$）の2つの結合された進化方程式が必要となる。
• 散逸・保存原理： 著者らは、熱力学第一法則（保存）と第二法則（散逸）を結びつける定理を確立した。全システムエネルギー（$E$）を、自由エネルギー（$F$）、圧力エネルギー（$P$）、運動エネルギー（$K$）、および壁の自由エネルギー（$W$）の総和として定義することで、エネルギー散逸を最小化することにより運動方程式を導出している。
• 修正ギブス・デュエム関係式： 既存モデルにおける空間勾配と関数微分との間の数学的不整合を解決するために、修正された等温ギブス・デュエム関係式が導出された。これは、速度依存のエントロピーと慣性仕事（inertial work）を明示的に考慮しており、自由エネルギーと速度エントロピーからの寄与を含む厳密な圧力の定義をもたらす。
• 運動エネルギーの定義： 運動エネルギーは、従来の $\frac{1}{2}\rho u^2$ ではなく、$K = \int \rho u \cdot u \, d\Omega$ （線形運動量密度 $\zeta = \rho u$）として定式化される。この区別により、運動量バランス（$u$ と $\rho$ の両方を変化させる）と速度散逸（$u$ のみを変化させる）の分離が可能になる。

主な貢献
この導出は、一般化された密度進化方程式と、3つの主要な理論的進展をもたらす。

1. 一般化されたベルヌーイの法則： 標準的なベルヌーイの原理（$p + \frac{1}{2}\rho u^2 = \text{const}$）は、特殊なケースであることが示されている。一般化された形式は、速度エントロピーの変化と密度の進化を考慮しており、$p + \rho u^2 = \text{constant}$ （特定の定常状態の非粘性条件下において）となり、密度差によって界面を横切る際にエネルギーの跳躍が生じることを認めている。
2. 普遍的な音速定式化： 自由エネルギーの密度に関する微分に基づく、音速の一般式が導出された：$c = \sqrt{d\Psi/d\rho}$。理想気体の場合、これは分子の平均熱速度に近似される。
3. 一般化された状態方程式（EOS）： 本理論は、密度進化方程式が、あらかじめ規定されたEOSを必要とせずに、流体のEOS（例：理想気体における $p \sim \rho \ln \rho$）を自然に生成することを実証している。これは標準的な拡散界面法とは対照的である。

結果と検証
論文では、以下の3つのケーススタディを通じて理論を検証している。

• 音波伝播： 導出された音速式は摂動解析と一致しており、空気中の音波の分散関係を正しく予測している。
• 理想気体のEOS： 定常状態において、モデルは理想気体の状態方程式を再現しており、密度進化方程式が熱力学的状態関係を本質的に捉えていることを確認している。
• 水・空気界面のダイナミクス： 管内を移動する平坦な水・空気界面の解析解により、本モデルが非単調な密度プロファイルを予測することが示された。具体的には、密度はバルク水の値から界面の中心での最大ピークに向かって増加し、その後、バルクガスの値へと減少する。この結果は、密度汎関数理論（DFT）の知見 [11, 12] と一致しており、線形補間法が予測する単調減少とは矛盾するものである。また、モデルは曲面界面におけるヤング・ラプラス方程式も正しく再現している。

意義
本論文は、圧縮性および非圧縮性の両方の領域に適用可能な、任意の高い密度比を持つ非混和流体のための一般的なフレームワークを提供すると主張している。密度進化を、過剰体積および粘性散逸と結合した動的な変数として扱うことで、本モデルはブシネスク近似と線形補間法の矛盾を解決している。著者らは、粘性散逸が本質的に密度を変化させることを示唆しており、散逸流における厳密な非圧縮性の仮定に疑問を投げかけている。著者らは、この理論が、密度比が極端な流体－流体界面における運動量およびエネルギー転送の物理的に一貫した記述を提供し、CFDに新たな窓を開くものであると断じている。