オイラー方程式 vs. ナビエ–ストークス方程式:違いは何か?
オイラー方程式は粘性を無視する。ナビエ–ストークス方程式はそれを含む。この一つの違いが、物理学、数学、そして百万ドルの問題を根本的に変える。
簡潔な答え
オイラー方程式は、内部摩擦のない — 粘性のない — 流体の運動を記述します。ナビエ–ストークス方程式は、同じ流体を粘性付きで記述します。
数学的に言えば、オイラー方程式とナビエ–ストークス方程式の違いは、一つの項 $\nu \Delta u$、すなわち粘性拡散項です。これを取り除けばナビエ–ストークスはオイラーになります。これを残せば、方程式は物理学と解析の両方を根本的に変える平滑化メカニズムを獲得します。
この一つの項が、煙が散逸する理由、境界層が表面に沿って形成される理由、そしてナビエ–ストークス・ミレニアム問題が対応するオイラーの問題とは異なる性格を持つ理由です。
$\mathbb{R}^3$ 上の非圧縮オイラー方程式は
$$\partial_t u + (u \cdot \nabla)u = -\nabla p, \qquad \nabla \cdot u = 0.$$
非圧縮ナビエ–ストークス方程式は
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0,$$
ここで $\nu > 0$ は動粘性係数です。形式的に、ナビエ–ストークスで $\nu = 0$ とすればオイラーが回復されます。しかしこの形式的な代入は、$\nu \to 0$ の極限が特異的であるという事実を隠しています:粘性項 $\nu \Delta u$ は系における最高次の空間微分であり、これを落とすとPDEの型と解が存在する関数空間が変わります。
二つの方程式の比較
比較が明確になるように、両方の方程式を標準的な非圧縮形式で示します:
オイラー:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p$$
ナビエ–ストークス:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u$$
両方の方程式は同じ左辺を共有しています:速度の時間変化率と非線形自己輸送項 $(u \cdot \nabla)u$。両方とも $\nabla \cdot u = 0$ により非圧縮性を強制します。唯一の構造的な違いは、ナビエ–ストークスの右辺にある粘性項 $\nu \Delta u$ です。
パラメータ $\nu$ は動粘性係数であり、流体の物理定数です。蜂蜜は大きな $\nu$ を持ちます。空気は小さな $\nu$ を持ちます。オイラー方程式は $\nu = 0$ に対応します:完全に摩擦のない理想化であり、境界から離れた高レイノルズ数流れの近似にはなりますが、通常の流体で厳密に実現されることはありません。
両系は双線形形式 $(u \cdot \nabla)u$ と、発散なし拘束条件により暗黙的に決定される圧力を共有しています。運動量方程式の発散をとり $\nabla \cdot u = 0$ を用いると、圧力のポアソン方程式
$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j),$$
が得られ、これは両方の系で同一です。圧力はいずれの場合も $u$ の非局所的汎関数です。
粘性項 $\nu \Delta u$ は各速度成分に作用する2次線形楕円型作用素です。これは放物型正則化であり、その存在によりナビエ–ストークスは半線形放物型系となり、一方で非圧縮オイラーは非局所的な圧力結合を持つ一次非線形輸送方程式となります。PDEの型のこの違いが、正則性理論におけるその後のほぼすべての違いを駆動しています。
粘性が物理的に何をするか
粘性は、流体の隣接層間の内部摩擦です。速く動く層が遅く動く層の隣にあるとき、粘性は両者の間で運動量を伝達し、速度差を滑らかにします。
これにより、ナビエ–ストークスをオイラーから分離する3つの主要な物理的帰結が生まれます:
- 散逸。運動エネルギーが熱に変換されます。コーヒーカップをかき混ぜて止めると、コーヒーはやがて静止します。オイラー方程式はこれを予測できません。理想化されたオイラーモデルでは、運動を熱に散逸させる粘性メカニズムがありません。
- 境界層。実際の流体は表面に付着します(すべりなし条件)。これにより壁面近くに速度が急変する薄い層 — 境界層 — が形成され、航空機の抗力、管内の摩擦、高速での乱流を生み出します。オイラー流はすべりなしではなくすべり条件を壁面で満たすため、粘性壁面抗力や実際の流体で見られる境界層物理学の多くを見逃します。
- 小スケールの平滑化。粘性は最も急峻な速度勾配を抑制します。それがなければ、流れが無限に微細な構造を発展させることを妨げるものはありません。この平滑化こそが、ナビエ–ストークスの正則性の問題をオイラーの場合と異なるものにしているのです。
$\mathbb{R}^3$(または周期領域)上のナビエ–ストークスのエネルギー恒等式は
$$\frac{d}{dt}\frac{1}{2}\|u\|_{L^2}^2 = -\nu \|\nabla u\|_{L^2}^2,$$
であり、運動エネルギーは単調に散逸します。オイラー($\nu = 0$)では、$u$ の $L^2$ ノルムは形式的に保存されます。
境界条件の観点では、すべりなし条件($u|_{\partial \Omega} = 0$)を持つナビエ–ストークスは適切に設定された放物型境界値問題ですが、オイラーは不浸透条件($u \cdot n = 0$)のみを要求します。$\nu \to 0$ での境界条件の不一致がプラントル境界層を生成し、これはプラントル(1904)以来研究されている特異摂動現象です。
物理的には、粘性は高周波フィルターとして機能します:フーリエモードを $\nu |k|^2$ の速度で減衰させ、小スケールを優先的に抑制します。このスペクトル的な減衰が、乱流におけるコルモゴロフ散逸スケール $\eta \sim (\nu^3 / \varepsilon)^{1/4}$ の背後にあるメカニズムです。完全なスケーリングの描像については、レイノルズ数と乱流をご覧ください。
オイラーは粘性ゼロのナビエ–ストークスにすぎないのか?
形式的にはそうです。ナビエ–ストークス方程式で $\nu = 0$ とすればオイラー方程式が得られます。しかしそこで立ち止まると誤解を招きます。
$\nu \to 0$ の極限は滑らかではなく特異的です。粘性は方程式における最高次の微分を担っています。それを取り除くことは小さな変化ではなく — PDEの根本的な性格を変えるのです。境界層は優雅に薄くなるわけではありません。乱流になりうるのです。ナビエ–ストークスのもとで滑らかだった解が、オイラーのもとでは全く異なる振る舞いを示す可能性があります。
アナロジー:力学系から摩擦を取り除くことは、単にものが滑りやすくなるだけではありません。系を質的に異なるものにしうるのです — 安定だった軌道がカオス的になり、摩擦によって強制されていた拘束が完全に解放されます。ここでも同じことが成り立ちます。
したがって、オイラーとナビエ–ストークスは数学的DNAを共有していますが、粘性ゼロ極限は流体力学における最も深く微妙な問題の一つであり、単純な代入ではありません。
非粘性極限 $\nu \to 0$ は特異摂動です:$\nu \Delta u$ が系における最高次の空間微分を担っており、$\nu = 0$ とするとPDEの次数が下がります。境界を持つ領域では、この極限はプラントルの境界層展開の妥当性に結びついており、これは劇的に破綻しうるものです(Grenier 2000、Gérard-Varet & Dormy 2010)。
$\mathbb{R}^3$ または $\mathbb{T}^3$(境界なし)では、状況はよりきれいですが、なお自明ではありません。オイラー解 $u^E$ が $[0,T]$ 上で滑らかであれば、$\nu \to 0$ でナビエ–ストークス解 $u^\nu$ は $L^2$ で $u^E$ に収束します(Kato 1972)。しかし、オイラー解がグローバルに滑らかであり続けるかどうか — 滑らかなデータに対して有限時間爆発が起きるかどうか — は、3次元ではそれ自体が未解決です。
この極限は乱流理論とも相互作用します。コルモゴロフの描像は散逸スケールを定義するために $\nu > 0$ を必要としますが、異常散逸 — $\nu \to 0$ でのエネルギー散逸の持続 — は中心的な予想です。オンサーガーの予想(現在は定理:Isett 2018、鋭さはBuckmaster–De Lellis–Székelyhidi–Vicol 2019による)は、粘性が存在しない場合にオイラー解がいつエネルギーを散逸しうるかを特徴づけています。
オイラー方程式がナビエ–ストークスの代わりに使われるのはいつか?
オイラー方程式は、流れの中の他の力と比べて粘性が無視できるほど小さいときに使われます。これは想像以上に頻繁に起こります:
- 表面から離れた高速空気力学。翼の表面から遠く離れた場所では、気流はほぼ非粘性的に振る舞います。エンジニアはバルク流にオイラーソルバーを使い、壁面近くに境界層モデルを組み込みます。
- 天体物理学的流れ。星間ガス、恒星内部、降着円盤は非常に大きなスケールで動作するため、分子粘性は無関係です(ただし乱流の実効粘性はそうではないかもしれません)。
- 圧縮性気体力学。衝撃波、爆轟、超音速流は、しばしば圧縮性オイラー方程式でモデル化されます。支配的な物理学は圧力と慣性であり、摩擦ではありません。
- 理論的解析。オイラーは、ナビエ–ストークスへの踏み台として、あるいはそれ自体が興味深いPDE系として研究されることがあります — 幾何学、渦力学、乱流の構造と深い関連があります。
しかし、摩擦、抗力、境界での振る舞いが重要な場合 — 管内流、表面近くの車両空気力学、血流、人間スケールの気象 — ではナビエ–ストークスが適切なモデルです。
オイラー方程式は、高レイノルズ数 $\mathrm{Re} = UL/\nu \gg 1$ での主要次の振る舞いを支配し、粘性効果は薄い境界層と内部せん断層に限定されます。境界から離れた高レイノルズ数バルク流では、オイラーがしばしば主要次の近似を与え、粘性補正は薄い境界層やせん断層に集中します。
圧縮性オイラー方程式 — 有限の伝播速度を持つ双曲型系 — は、衝撃波形成やリーマン問題を含む気体力学の標準モデルです。これらは上で議論した非圧縮オイラー方程式とは異なります:圧縮性オイラーは真に双曲型ですが、非圧縮オイラーは非局所的な圧力結合と無限の伝播速度を持ちます。
数学的解析において、オイラーは粘性消失問題の極限対象として、また独自の正則性理論、保存量(ヘリシティ、微分同相群上のオイラー–アーノルド枠組みを通じたカシミール関数)、幾何学的力学との関連を持つ豊かなPDE系として機能します。
この違いが正則性にとって何を意味するか
正則性の問題は、オイラー対ナビエ–ストークスのギャップが最も重要な場所です。
ナビエ–ストークス・ミレニアム問題は問います:3次元で滑らかで振る舞いの良い流れから出発した場合、ナビエ–ストークスの解はすべての時刻で滑らかであり続けるか? 誰も知りません。
同じ問いがオイラーの3次元においても未解決です。しかし、二つの問題は根本的に異なる性格を持っています:
- ナビエ–ストークスには粘性が味方として働いています — 常に平滑化し、常にエネルギーを散逸し、常に最も急峻な勾配を抑制しています。問題は、この平滑化が非線形項が特異点を作り出すことを防ぐのに十分に強いかどうかです。
- オイラーには平滑化が全くありません。非線形項は勾配を何の制止もなく増幅できます。滑らかな初期データから3次元オイラーが有限時間特異点を発展させるかどうかは未解決です。
2次元では、滑らかな初期データに対して両方の方程式がグローバルに適切に設定されています。2次元の話は解決済みです。謎が存在するのは3次元であり、両方の方程式に対して、しかし異なる形で存在しています。
正則性の状況:
2次元:滑らかな解の大域存在と一意性は、両方の系で知られています。滑らかなデータを持つ2次元オイラーについて、Wolibner(1933)がヘルダー空間での大域存在を証明し、Yudovich(1963)が有界渦度データに対する一意性を確立しました。2次元ナビエ–ストークスについては、ラディジェンスカヤの不等式と渦度の最大値原理から大域正則性が従います。
3次元ナビエ–ストークス:ルレイ(1934)は $L^2$ における弱解の大域存在を証明しましたが、一意性と正則性は未解決のままです。カファレリ–コーン–ニーレンバーグの定理(1982)は、特異点集合が1次元放物型ハウスドルフ測度ゼロであることを示しています — 爆発が起きるとしても、極めて疎です。粘性項は重要なアプリオリ評価 $\int_0^T \|\nabla u\|_{L^2}^2 \, dt \leq C(u_0)$ を与えますが、この $H^1$ 制御は3次元スケーリングに対して劣臨界的であり、ブートストラップ論法を閉じるには不十分です。超臨界性のギャップについては、なぜ難しいのかをご覧ください。
3次元オイラー:滑らかなデータに対する大域理論は存在しません。ソボレフ空間 $H^s$、$s > 5/2$ における局所適切性は古典的です(Kato 1972、Kato–Ponce 1988)。ビール–加藤–マイダの判定条件(1984)は、爆発検出を $\int_0^T \|\omega\|_{L^\infty} \, dt$ に帰着させます:$\int_0^T \|\omega\|_{L^\infty} \, dt < \infty$ ならば、かつそのときに限り、解は $[0,T]$ 上で滑らかです。爆発には渦度の上限が時間で非可積分になるほど急速に増大することが必要です。Elgindi(2021、Annals of Mathematics)は $C^{1,\alpha}$ データに対する有限時間特異点形成を証明しました — 画期的な成果ですが、滑らか($C^\infty$)な閾値には達していません。滑らかなオイラー解が3次元で爆発するかどうかは依然未解決です。
粘性、乱流、そしてカスケード
乱流は、オイラー対ナビエ–ストークスの比較が最も物理的に鮮明になる場所です。
乱流の流れでは、エネルギーは大スケール(管の大きさ、翼、嵐のスケール)で流入し、どんどん小さな渦へとカスケードしていきます。これがエネルギーカスケードです。カスケードの底で、粘性がついに運動エネルギーを熱に変換し、プロセスを停止させます。
オイラー方程式はカスケードの上部と中間 — 大スケールと慣性領域のダイナミクス — をモデル化できます。しかし、粘性が作用する底部はモデル化できません。粘性がなければ、自然な最小スケールは存在せず、エネルギーの行き場がありません。
これが、乱流のモデリングがほぼ常にオイラーではなくナビエ–ストークスを伴う理由です。レイノルズ数 $\mathrm{Re} = UL/\nu$ は、カスケードの幅を制御します:高い $\mathrm{Re}$ は、エネルギー注入と粘性散逸の間に多くのディケードのスケールがあることを意味します。実際の乱流は、非粘性カスケードと粘性カットオフの相互作用の中に存在します。
コルモゴロフの1941年の理論は、大スケール強制も粘性散逸も支配的でない慣性領域において、エネルギースペクトル $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ を予測します。散逸スケール $\eta = (\nu^3/\varepsilon)^{1/4}$ がこの範囲の下限を設定します。
オイラー方程式は $\nu = 0$ の極限、すなわち慣性領域がすべてのスケールに及ぶ場合を記述します。この極限でエネルギー散逸が持続するかどうか — 異常散逸 — がオンサーガーの予想です。剛性側($u \in C^{0,\alpha}$、$\alpha > 1/3$ ならば散逸なし)はConstantin–E–Titi(1994)により証明されました。柔軟性側($C^{1/3}$ 以下の散逸的オイラー解の存在)は、De Lellis–Székelyhidiの凸積分プログラムに基づいてIsett(2018)により完成されました。
ナビエ–ストークスでは、カスケードの描像がエネルギー収支に組み込まれています:$\varepsilon = \nu \langle |\nabla u|^2 \rangle$。重要な未解決問題は、この統計理論が数学的に根拠を持つのに十分な時間、ナビエ–ストークス解が滑らかであり続けるかどうかです。存在と滑らかさの問題は、部分的に、粘性の散逸メカニズムがすべてのスケールですべての時刻にわたってカスケードを正則化するのに十分に強いかどうかを問うています。
まとめ:一つの項、二つの異なる世界
オイラー方程式とナビエ–ストークス方程式の違いは一つの項 $\nu \Delta u$ です。しかしその項がすべてを変えます。
| オイラー | ナビエ–ストークス | |
|---|---|---|
| 粘性 | なし($\nu = 0$) | あり($\nu > 0$) |
| エネルギー | 保存(形式的に) | 散逸 |
| 境界層 | 非粘性モデルの範囲外 | あり(すべりなし条件) |
| PDEの型 | 一次非線形 + 非局所的圧力 | 二次放物型 + 非局所的圧力 |
| 2次元正則性 | 解決済み(常に滑らか) | 解決済み(常に滑らか) |
| 3次元正則性 | 未解決 | 未解決(ミレニアム問題) |
オイラーは単純化されたナビエ–ストークスではありません。構造の大部分を共有するが根本的に異なる系です。この比較は実用的な力も持っています:間違ったモデルの選択 — 粘性が重要な場面でオイラーを使うこと、あるいは粘性が不要な場面でナビエ–ストークスを使うこと — は、有用なシミュレーションと無用なシミュレーションの分かれ目になりえます。完全な方程式についてはナビエ–ストークス方程式とは何か?を、障害についてはなぜ難しいのかを、賞金についてはミレニアム問題を、非圧縮流 vs. 圧縮流については非圧縮 vs. 圧縮ナビエ–ストークスをご覧ください。
粘性項 $\nu \Delta u$ は、一次非線形系(オイラー)を半線形放物型系(ナビエ–ストークス)に変換する放物型正則化です。エネルギー散逸、高次のアプリオリ評価、そしてナビエ–ストークス理論の多くの基盤となる解析的半群構造を提供します(Leray 1934、Fujita–Kato 1964)。
しかしこの正則化は、3次元で大域正則性の論法を閉じるのに十分ではありません。クレイ・ミレニアム問題は、ナビエ–ストークスにおける放物型平滑化が、すべてのスケールにわたるすべての時刻で非線形エネルギー転送を制御できるかどうかを正確に問うています。オイラーに対する並行する問い — 平滑化が全く存在しないことが、滑らかなデータからの有限時間爆発につながるかどうか — も同様に未解決であり、同様に根本的です。
両方の問題は、流体の数学的理論の中心にあります。両者の比較は、粘性が何を買い、何を買わないかを明確にし、どちらの系に対しても3次元正則性の問題が解析学における最も困難な未解決問題の一つである理由を示しています。