← 戻る ライブラリ › 微分方程式と力学系 › 偏微分方程式 微分方程式と力学系
ナビエ–ストークス方程式 流体の流れを支配する方程式であり、三次元における解の滑らかさはミレニアム懸賞問題の一つである。
直観 直感:なぜ流体の流れは予測が難しいのか? パドルの周りで渦を巻く水、上昇する煙の渦、羼の上を流れる空気 — これらはすべて一つの方程式系で記述される。ナビエ–ストークス方程式は、速度場 u u u (各点で流体がどれだけ速く、どの方向に動くか)と圧力場 p p p を追跡し、粘性 ν \nu ν (流体がどれだけ「粘っこい」か)を用いて、押し、ねじり、減速させる力の釣り合いを表す。
らぞん面(ヘリコイド) — 空間をねじれながら伸びる渦線を様式化した図。三次元では、渦管が引き伸ばされるとともに回転が速くなることがある(渦の引き伸ばし)。これは三次元の滑らかさ問題を非常に難しくしている極端な振る舞いを引き起こすと考えられている機構である。 大学 方程式 定義: 非圧縮性ナビエ–ストークス方程式
非圧縮性の粘性流体が領域を溡たすとき、速度場 u ( x , t ) u(x,t) u ( x , t ) と圧力 p ( x , t ) p(x,t) p ( x , t ) (一定の密度で規格化)は、運動量のつり合いと非圧縮性の制約を満たす。これらは合わせて、流体要素へのニュートンの第二法則と、流体が移動しても体積が保存されるという要求を表している。
∂ t u + ( u ⋅ ∇ ) u = − ∇ p + ν Δ u \partial_t u+(u\cdot\nabla)u=-\nabla p+\nu\Delta u ∂ t u + ( u ⋅ ∇ ) u = − ∇ p + ν Δ u 左から順に読むと: ∂ t u \partial_t u ∂ t u は固定した点での速度の時間変化、( u ⋅ ∇ ) u (u\cdot\nabla)u ( u ⋅ ∇ ) u は移流項 — 流体が自身の運動量を運びながら動くことを表し、困難の大部分を生む非線形項である。− ∇ p -\nabla p − ∇ p は圧力差によって流体を高圧側から低圧側へ押す力、ν Δ u \nu\Delta u ν Δ u は粘性拡散であり、急激な速度勾配を滑らかにする。
∇ ⋅ u = 0 \nabla\cdot u=0 ∇ ⋅ u = 0 非圧縮性条件 ∇ ⋅ u = 0 \nabla\cdot u=0 ∇ ⋅ u = 0 は、流体がどこかに溯まったり消えたりしないこと — 移動しても体積が保存されることを意味する。結果として得られる流れが滑らか(層流)になるか乱雑(乱流)になるかは、主にレイノルズ数 R e = U L / ν \mathrm{Re}=UL/\nu Re = U L / ν によって決まる。これは流れの代表速度 U U U と代表長さ L L L を粘性 ν \nu ν と比較したものである:R e \mathrm{Re} Re が小さいと粘性が支配的で乱れは減衰し、R e \mathrm{Re} Re が大きいと慣性が支配的となり小さな乱れが乱流へと成長しうる。
レイノルズ数による流れの体制 体制 典型Re 振る舞い 層流 R e ≲ 10 3 \mathrm{Re}\lesssim 10^{3} Re ≲ 1 0 3 滑らかな平行層、層間の混合なし 遷移 10 3 ≲ R e ≲ 10 4 10^{3}\lesssim\mathrm{Re}\lesssim 10^{4} 1 0 3 ≲ Re ≲ 1 0 4 こちらは滑らかな流れの中に断続的に乱れが発生する 乱流 R e ≳ 10 4 \mathrm{Re}\gtrsim 10^{4} Re ≳ 1 0 4 多すかのスケールにわたるカオス的な渦、混合と抗力が急増
大学 主要な定理 有限の運動エネルギーを持つ任意の発散なし初期速度 u 0 u_0 u 0 に対し、全空間上の3次元ナビエ–ストークス方程式は、エネルギー不等式を満たす大域時間弱解 u u u を持つ。
なぜ正しいのか? 非線形項のため方程式を厳密に解くことは難しいが、粘性は常に運動エネルギーを減らし続ける。この単一の定量的事実だけで、極限が一意か滑らかかは分からなくても、近似解を構築して極限をとるのに十分である。
証明 ガレルキン射影によって方程式を近似する:ストークス演算子の固有関数を基底として選び、方程式を最初の n n n 個が張る空間へ射影し、得られる有限次元常微分方程式系を解いて近似解 u n u_n u n を得る。標準的な常微分方程式の理論により、短い時間区間での解の存在が得られる。
ガレルキン方程式に u n u_n u n 自身を作用させる。非線形項と圧力項は消える(非圧縮性の制約を使えば u n u_n u n に直交するため)ため、エネルギー等式 1 2 d d t ∥ u n ∥ L 2 2 + ν ∥ ∇ u n ∥ L 2 2 = 0 \frac{1}{2}\frac{d}{dt}\|u_n\|_{L^2}^{2}+\nu\|\nabla u_n\|_{L^2}^{2}=0 2 1 d t d ∥ u n ∥ L 2 2 + ν ∥∇ u n ∥ L 2 2 = 0 が得られる。時間積分すると、u n u_n u n は n n n に依らず L ∞ ( L 2 ) ∩ L 2 ( H 1 ) L^{\infty}(L^2)\cap L^2(H^1) L ∞ ( L 2 ) ∩ L 2 ( H 1 ) で一様に有界であり、特に短時間解は全時間に延長できる。
この一様的な有界性により、u n u_n u n の部分列で L ∞ ( L 2 ) L^{\infty}(L^2) L ∞ ( L 2 ) で弱・収束し、L 2 ( H 1 ) L^2(H^1) L 2 ( H 1 ) で弱収束するものを取り出せ、ある極限 u u u に収束する。オビゥーン–リオンのコンパクト性補題により、これは有界な時空領域での L 2 L^2 L 2 強収束に強化され、二次非線形項 ( u n ⋅ ∇ ) u n → ( u ⋅ ∇ ) u (u_n\cdot\nabla)u_n\to(u\cdot\nabla)u ( u n ⋅ ∇ ) u n → ( u ⋅ ∇ ) u への極限移行に必要なものである。
ガレルキン方程式の弱形式で極限を取ることで、u u u が超関数の意味でナビエ–ストークス方程式を満たすことが示され、極限下でのノルムの弱下半連続性によりエネルギー不等式が保たれる。この u u u がルレーの大域弱解であり、それが一意かつ滑らかかどうかが未解決のミレニアム問題の本質である。
二次元空間において、有限エネルギーを持つ任意の滑らかな発散なし初期速度に対し、ナビエ–ストークス方程式はすべての時間で滑らかであり続ける一意解を持つ — 未解決の3次元の場合とは対照的に、爆発は決して起こらない。
なぜ正しいのか? 3次元の渦度は ω t + u ⋅ ∇ ω = ( ω ⋅ ∇ ) u + ν Δ ω \omega_t+u\cdot\nabla\omega=(\omega\cdot\nabla)u+\nu\Delta\omega ω t + u ⋅ ∇ ω = ( ω ⋅ ∇ ) u + ν Δ ω に従い、項 ( ω ⋅ ∇ ) u (\omega\cdot\nabla)u ( ω ⋅ ∇ ) u は引き伸ばしによって渦度を増幅させうる。二次元では渦度は運動平面に垂直なスカラー量であり、この引き伸ばし項は恒等的に消え、輸送と拡散だけが残る — 無制限に増幅させるものは何もない。
証明 渦度(2次元ではスカラー)を ω = ∂ x u 2 − ∂ y u 1 \omega=\partial_x u_2-\partial_y u_1 ω = ∂ x u 2 − ∂ y u 1 と書き、運動方程式を ∂ t ω + u ⋅ ∇ ω \partial_t\omega+u\cdot\nabla\omega ∂ t ω + u ⋅ ∇ ω に代入する。u u u が非圧縮性かつ二次元であるため、( u ⋅ ∇ ) u (u\cdot\nabla)u ( u ⋅ ∇ ) u から生じる ω \omega ω を引き伸ばすことになるすべての項が直接計算で消え、ソース項のない輸送拡散方程式 ω t + u ⋅ ∇ ω = ν Δ ω \omega_t+u\cdot\nabla\omega=\nu\Delta\omega ω t + u ⋅ ∇ ω = ν Δ ω が残る。
ソースのない輸送拡散方程式は最大原理に従う:流体粒子の軌道 X ˙ ( t ) = u ( X ( t ) , t ) \dot X(t)=u(X(t),t) X ˙ ( t ) = u ( X ( t ) , t ) に沿って d d t ω ( X ( t ) , t ) = ν Δ ω ( X ( t ) , t ) \frac{d}{dt}\omega(X(t),t)=\nu\Delta\omega(X(t),t) d t d ω ( X ( t ) , t ) = ν Δ ω ( X ( t ) , t ) となり、拡散は最大値を増やしたり最小値を減らしたりできない。したがってすべての t ≥ 0 t\ge0 t ≥ 0 で ∥ ω ( ⋅ , t ) ∥ L ∞ ≤ ∥ ω 0 ∥ L ∞ \|\omega(\cdot,t)\|_{L^\infty}\le\|\omega_0\|_{L^\infty} ∥ ω ( ⋅ , t ) ∥ L ∞ ≤ ∥ ω 0 ∥ L ∞ となり、渦度は決して爆発しない。
渦度が時間に一様に有界であれば、流を関数(ω \omega ω から u u u を復元)に楕円型正則性を適用することで、u u u のすべての空間微分が ∥ ω ∥ L ∞ \|\omega\|_{L^\infty} ∥ ω ∥ L ∞ によって抑えられ、u u u は最初に属していたすべてのソボレフ空間に留まる。
この論議を繰り返すことで — 有界な渦度が速度勾配を有界にし、それが方程式自体を通じて高次微分を制御する — すべての t > 0 t>0 t > 0 で解が滑らかであり続け、有限時間での爆発は起こりえないことが示される。これはまさに、渦の引き伸ばし項 ( ω ⋅ ∇ ) u (\omega\cdot\nabla)u ( ω ⋅ ∇ ) u が存在する三次元では未知である事柄である。
大学 実世界での応用と具体例 ナビエ–ストークス方程式は数値流体力学の中核である:航空機や船舶の船体は、建造前に候補設計周りの気流や水流をシミュレーションして形状が決められる。天気・気候モデルはこれらの方程式のバージョンを大気と海洋全体で積分する。循環器系の研究者は手術計画のために動脈内の血流をシミュレーションし、天体物理学者は圧縮性の変形を用いて降着円盤や恒星内部のガス流をモデル化する。
例: 2枚の平板間の速度分布(平面ポアズイュ流)
一定の圧力勾配によって駆動される定常・非圧縮性の流れが、y = − h y=-h y = − h と y = h y=h y = h にある2つの無限平行平板の間を、両壁でのノースリップ条件で流れる。駆動項(定数)を G = − 1 ρ d p d x > 0 G=-\frac{1}{\rho}\frac{dp}{dx}>0 G = − ρ 1 d x d p > 0 、動粘性係数を ν \nu ν として、速度分布 u ( y ) u(y) u ( y ) を求めよ。
解答 y y y のみに依存する定常流で u = ( u ( y ) , 0 ) u=(u(y),0) u = ( u ( y ) , 0 ) とすると、u u u は流れの方向 x x x に沿って変化しないため、非線形項 ( u ⋅ ∇ ) u (u\cdot\nabla)u ( u ⋅ ∇ ) u は消える。運動方程式は 0 = − 1 ρ d p d x + ν u ′ ′ ( y ) 0=-\frac{1}{\rho}\frac{dp}{dx}+\nu u''(y) 0 = − ρ 1 d x d p + ν u ′′ ( y ) 、すなわち ν u ′ ′ ( y ) = − G \nu u''(y)=-G ν u ′′ ( y ) = − G に帰着する。
y y y について二回積分すると u ( y ) = − G 2 ν y 2 + C 1 y + C 2 u(y)=-\frac{G}{2\nu}y^{2}+C_1y+C_2 u ( y ) = − 2 ν G y 2 + C 1 y + C 2 となり、定数 C 1 , C 2 C_1,C_2 C 1 , C 2 は境界条件で決まる。
ノースリップ条件は u ( h ) = u ( − h ) = 0 u(h)=u(-h)=0 u ( h ) = u ( − h ) = 0 を要求する。両方を代入して引くと C 2 C_2 C 2 が消え C 1 = 0 C_1=0 C 1 = 0 となり(両式は線形項の符号のみが異なるため)、加えると C 2 = G h 2 2 ν C_2=\frac{Gh^{2}}{2\nu} C 2 = 2 ν G h 2 を得る。
代入し直すと、速度分布は u ( y ) = G 2 ν ( h 2 − y 2 ) u(y)=\frac{G}{2\nu}\left(h^{2}-y^{2}\right) u ( y ) = 2 ν G ( h 2 − y 2 ) :壁でゼロ、中心線 y = 0 y=0 y = 0 で最大となる古典的なポアズイュ分布である。
例: レイノルズ数:泳ぐ細菌と高速道路上の自動車
(a) 水中を U = 2 × 10 − 5 m/s U=2\times10^{-5}\text{ m/s} U = 2 × 1 0 − 5 m/s で泳ぐ長さ L = 2 × 10 − 6 m L=2\times10^{-6}\text{ m} L = 2 × 1 0 − 6 m の細菌(ν ≈ 10 − 6 m 2 / s \nu\approx10^{-6}\text{ m}^2/\text{s} ν ≈ 1 0 − 6 m 2 / s )、および (b) 空気中を U = 30 m/s U=30\text{ m/s} U = 30 m/s で走る長さ L = 4 m L=4\text{ m} L = 4 m の自動車(ν ≈ 1.5 × 10 − 5 m 2 / s \nu\approx1.5\times10^{-5}\text{ m}^2/\text{s} ν ≈ 1.5 × 1 0 − 5 m 2 / s )のレイノルズ数 R e = U L / ν \mathrm{Re}=UL/\nu Re = U L / ν を見積もれ。それぞれの値はナビエ–ストークス方程式のどの項が支配的かについて何を語るか?
解答 細菌の場合: R e = ( 2 × 10 − 5 ) ( 2 × 10 − 6 ) 10 − 6 = 4 × 10 − 5 \mathrm{Re}=\dfrac{(2\times10^{-5})(2\times10^{-6})}{10^{-6}}=4\times10^{-5} Re = 1 0 − 6 ( 2 × 1 0 − 5 ) ( 2 × 1 0 − 6 ) = 4 × 1 0 − 5 で、1よりはるかに小さい。
自動車の場合: R e = ( 30 ) ( 4 ) 1.5 × 10 − 5 = 8 × 10 6 \mathrm{Re}=\dfrac{(30)(4)}{1.5\times10^{-5}}=8\times10^{6} Re = 1.5 × 1 0 − 5 ( 30 ) ( 4 ) = 8 × 1 0 6 で、1よりはるかに大きい。
R e ≪ 1 \mathrm{Re}\ll1 Re ≪ 1 と小さいということは、粘性項 ν Δ u \nu\Delta u ν Δ u が非線形項 ( u ⋅ ∇ ) u (u\cdot\nabla)u ( u ⋅ ∇ ) u を圧倒することを意味する:細菌にとって水ははちみつのように粘っこく、慣性は無関係であり、このような低 R e \mathrm{Re} Re では単に前後に振るだけでは正正运動が生まれないため、もみ抹じのように進む必要がある。
R e ≫ 1 \mathrm{Re}\gg1 Re ≫ 1 と大きいということは、自動車表面近くの薄い境界層を除きほぼすべての場所で非線形移流項が粘性を圧倒することを意味する:流れは容易に剥離し、後方の尾流で乱流化する。だからこそ高速走行時に空気抵抗と形状設計が重要になる。
よくある誤り. ミレニアム問題のよくある誤解は、数学者はナビエ–ストークス方程式の解がそもそも「存在する」かどうかさえ分かっていないというものだ。その部分は既に解決済みであり、1934年にルレーが3次元で大域弱解の存在を証明している。未解決の問題はもっと鳴い:滑らかな初期データに対して、その解が一意であり、すべての時間で滑らかさ(爆発がない)を保つかどうかである。存在は既知だが、一意性と滑らかさは未知である。 研究の最前線 2026年時点
クレイ数学研究所のミレニアム懸賞問題 — 3次元ナビエ–ストークス方程式の滑らかで有限エネルギーの解が常に滑らかなままであるか、有限時間で爆発し得るか — は2026年現在未解決である。Caffarelli、Kohn、Nirenberg(1982)は、特異点が形成されるならばその集合は小さい(放物線的ハウスドルフ次元が高々す1)ことを示し、Beale–Kato–Majda基準(1984)は最大渦度の時間積分が発散する場合にのみ爆発が起こりうることを示し、特異点を排除または発見するための具体的な目標を与えている。 2016年、テレンス・タオは、3次元方程式のもっともらしい平均化版が有限時間で爆発し得ることを証明し、本当の方程式の大域正則性の証明には、その粗い大きさだけではなく、正確な非線形項の構造を利用する必要があることを示した。現在の研究は、特異点に近い渦構造の高解像度数値探索と、真の特異点を見つけるか排除することを目指す新しい関数空間手法を組み合わせている。
非圧縮性条件 ∇ ⋅ u = 0 \nabla\cdot u=0 ∇ ⋅ u = 0 は物理的に何を意味するか?
速度がどこでも一定である 流体体積が保存される:どの点でも正正の流入・流出はない 圧力がどこでもゼロである 流れは常に乱流である
ある管内流の代表速度が U = 0.5 m/s U=0.5\text{ m/s} U = 0.5 m/s 、代表長さが L = 0.02 m L=0.02\text{ m} L = 0.02 m 、動粘性係数が ν = 1 × 10 − 6 m 2 / s \nu=1\times10^{-6}\text{ m}^2/\text{s} ν = 1 × 1 0 − 6 m 2 / s である。R e = U L / ν \mathrm{Re}=UL/\nu Re = U L / ν を計算し、流れが層流か乱流かを答えよ。
R e = 100 \mathrm{Re}=100 Re = 100 , 層流R e = 10,000 \mathrm{Re}=10{,}000 Re = 10 , 000 , 乱流R e = 0.01 \mathrm{Re}=0.01 Re = 0.01 , 層流R e = 500,000 \mathrm{Re}=500{,}000 Re = 500 , 000 , 乱流巡航速度での羼上の気流シミュレーションが、はちみつのような遅い流れに比べて計算コストが高いのはなぜか?
高速では圧力場 p p p が存在しないから 高レイノルズ数では非線形項 ( u ⋅ ∇ ) u (u\cdot\nabla)u ( u ⋅ ∇ ) u が支配的となり、すべて解像すべき多すかのスケールの微細な乱流渦を生じるから 高速では粘性 ν \nu ν が負になるから 歩行速度を超えると流体は非圧縮性でなくなるから ナビエ–ストークスのミレニアム問題は正確に何を問うのか?
3次元で弱解がそもそも存在するかどうか 滑らかで有限エネルギーの初期データが与えられたとき、3次元で滑らかな解がすべての時間一意に存在するか、あるいは有限時間で爆発しうるか 方程式がコンピュータなしで手計算で解けるかどうか 十分大きな粘性 ν \nu ν を選べば乱流を排除できるかどうか