ナビエ–ストークス方程式

連続体力学
法則 質量保存の法則 運動量保存の法則 エネルギー保存の法則 クラウジウス–デュエムの不等式 固体力学 固体 · 変形 · 弾性 · 弾性波 · 弾塑性 · 塑性 · フックの法則 · 応力 · ひずみ · 有限変形理論 · レオロジー · 粘弾性 · 超弾性 流体力学 流体 · 流体静力学 流体動力学 · 粘度 · ニュートン流体 非ニュートン流体 表面張力 科学者 ニュートン · ストークス · ナビエ · コーシー · フック · ベルヌーイ
表 話 編 歴

ナビエ–ストークス方程式（ナビエ–ストークスほうていしき、英: Navier–Stokes equations）は、流体の運動を記述する2階非線型偏微分方程式であり、流体力学で用いられる。[1][2]アンリ・ナビエとジョージ・ガブリエル・ストークスによって導かれた[3][4]。日本語の文献だと「NS方程式」とも略される[5]。ナビエ・ストークス方程式は、ニュートン力学における運動の第2法則に相当する。

導出

流体の質量保存の法則と運動量保存の法則を表す連続の方程式

${\displaystyle {\frac {\partial \rho }{\partial t}}+\operatorname {div} (\rho {\boldsymbol {v}})=0}$

${\displaystyle {\frac {\partial (\rho {\boldsymbol {v}})}{\partial t}}+\operatorname {div} (\rho {\boldsymbol {v}}{\boldsymbol {v}})=\operatorname {div} {\boldsymbol {\sigma }}+\rho {\boldsymbol {g}}}$

を用いると、流れの速度場 ${\displaystyle {\boldsymbol {v}}}$ の物質微分は

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}={\frac {1}{\rho }}\operatorname {div} {\boldsymbol {\sigma }}+{\boldsymbol {g}}}$

と導かれる。ここで、 ${\displaystyle \rho }$ は密度場、 ${\displaystyle {\boldsymbol {\sigma }}}$ は応力場、 ${\displaystyle {\boldsymbol {g}}}$ は流体の単位質量あたりに作用する外力場（加速度場）である。

ここで、ニュートン流体を仮定すれば、応力場が

${\displaystyle {\boldsymbol {\sigma }}=\left(-p+\chi \Theta \right)\mathbf {1} +\mu \left({\boldsymbol {e}}-{\frac {2}{3}}\,\Theta \,\mathbf {1} \right)=(-p+(\chi -{\frac {2}{3}}\,\mu )\Theta )\mathbf {1} +\mu {\boldsymbol {e}}}$

で与えられる。ただし、 ${\displaystyle p}$ は圧力（静圧）、 ${\displaystyle \chi }$ は体積粘性率、 ${\displaystyle \mu }$ は剪断粘性率である。 ${\displaystyle {\boldsymbol {e}}}$ は対称化した速度勾配（英語版）で、デカルト座標の下で成分表示をすれば

${\displaystyle e_{ab}={\frac {\partial v_{a}}{\partial x_{b}}}+{\frac {\partial v_{b}}{\partial x_{a}}}}$

で表され、 ${\displaystyle \Theta }$ は速度場の発散

${\displaystyle \Theta =\operatorname {div} {\boldsymbol {v}}={\frac {1}{2}}\operatorname {tr} {\boldsymbol {e}}}$

である。

この形の応力場 ${\displaystyle {\boldsymbol {\sigma }}}$ を用いると、速度場 ${\displaystyle {\boldsymbol {v}}}$ の物質微分が

${\displaystyle {\begin{aligned}{\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=&-{\frac {1}{\rho }}\operatorname {grad} p+{\frac {\mu }{\rho }}\Delta {\boldsymbol {v}}+{\frac {\chi +{\frac {1}{3}}\,\mu }{\rho }}\operatorname {grad} \Theta +{\frac {\Theta }{\rho }}\operatorname {grad} (\chi +{\frac {1}{3}}\,\mu )\\&+{\frac {1}{\rho }}\operatorname {grad} ({\boldsymbol {v}}\cdot \operatorname {grad} \mu )+{\frac {1}{\rho }}\operatorname {rot} ({\boldsymbol {v}}\times \operatorname {grad} \mu )-{\frac {1}{\rho }}\,{\boldsymbol {v}}\Delta \mu +{\boldsymbol {g}}\\\end{aligned}}}$

で与えられる。この方程式がナビエ–ストークス方程式である[1][2][6][注 1]。この3本[注 2]の連立偏微分方程式を解いて3次元ベクトル ${\displaystyle {\boldsymbol {v}}}$ とスカラー ${\displaystyle p}$ の計4つの未知関数の一般解が常に存在する（もしくは存在しないケースがある）ことを証明せよ、という問題が「ナビエ–ストークス方程式の解の存在と滑らかさ（ミレニアム懸賞問題の1つ）」である。加えて、それらの解が「時間大域的かつ滑らかな解」なのかどうかも、非常に重要な論点となる。

なお、速度場の物質微分の第二項は「対流項」あるいは「移流項」と呼ばれ、ベクトル解析の公式により、

${\displaystyle ({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=\operatorname {grad} \left({\frac {1}{2}}{\boldsymbol {v}}^{2}\right)-{\boldsymbol {v}}\times {\boldsymbol {\omega }}}$

と変形することができる。ここで ω は速度場の回転

${\displaystyle {\boldsymbol {\omega }}=\operatorname {rot} {\boldsymbol {v}}}$

であり、渦度と呼ばれる[7]。

単純化した方程式

ナビエ–ストークス方程式は非線形であり、複雑過ぎるので解を求めることは困難である[1][2][6]。このため、いくつかの仮定をして問題を簡単化することが多い[8]。しかし簡単化された方程式ですら解析的な方法では解が得られないことが普通であり、解の存在性などの定性的な議論[9]を超えて、具体的な解の様子を知るためには、ほとんどの場合に数値的な近似解法が必要になる[注 3][10]。

非圧縮性流れ

非圧縮性流れ[11]では、速度場の発散 Θ がゼロなので、速度場の発散を含む項を落として

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=-{\frac {1}{\rho }}\operatorname {grad} p+{\frac {\mu }{\rho }}\Delta {\boldsymbol {v}}+{\frac {1}{\rho }}\operatorname {grad} ({\boldsymbol {v}}\cdot \operatorname {grad} \mu )+{\frac {1}{\rho }}\operatorname {rot} ({\boldsymbol {v}}\times \operatorname {grad} \mu )-{\frac {1}{\rho }}{\boldsymbol {v}}\,\Delta \mu +{\boldsymbol {g}}}$

となる。

粘性率が一定の流れ

粘性率 μ や χ は温度や圧力の関数であり一定ではないが、粘性率を定数と仮定する場合は、粘性率の勾配を含む項を落として、

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=-{\frac {1}{\rho }}\operatorname {grad} p+{\frac {\mu }{\rho }}\Delta {\boldsymbol {v}}+{\frac {\chi +{\frac {1}{3}}\,\mu }{\rho }}\operatorname {grad} \Theta +{\boldsymbol {g}}}$

となる。また、体積粘性率 χ は非常に小さいので、χ = 0 と仮定すると、

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=-{\frac {1}{\rho }}\operatorname {grad} p+\nu \Delta {\boldsymbol {v}}+{\frac {\nu }{3}}\operatorname {grad} \Theta +{\boldsymbol {g}}}$

となる（ストークスの仮説）。ここで ν = μ/ρ は動粘性率である。

粘性率が一定の非圧縮性流れ

粘性率が一定の非圧縮性流れでは、

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=-{\frac {1}{\rho }}\operatorname {grad} p+\nu \Delta {\boldsymbol {v}}+{\boldsymbol {g}}}$

となる。ここで ν = μ/ρ は動粘性率である。各項はそれぞれ、

• 左辺 - 第1項：時間[微分]項、第2項：移流項（対流項）
• 右辺 - 第1項：圧力項、第2項：粘性項（拡散項）、第3項：外力項

と呼ばれる。外力項は、状況によって、重力をはじめ浮力・表面張力・電磁気力などが該当する。

ストークス流れ（クリープ流れ）

粘性率が一定の非圧縮性流れのうち、流体の速度が遅かったりスケールが小さいなど、レイノルズ数が小さい流れを特にストークス流れあるいはクリープ流れという。ストークス流れでは、非線型である対流項 ${\displaystyle ({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}}$ が無視できて、

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}=-{\frac {1}{\rho }}\operatorname {grad} p+\nu \Delta {\boldsymbol {v}}+{\boldsymbol {g}}}$

となる。この式はストークス方程式（Stokes equations）と呼ばれている[12][13]。

オイラー方程式

粘性のない（χ = μ = 0）流れでは

${\displaystyle {\frac {D{\boldsymbol {v}}}{Dt}}={\frac {\partial {\boldsymbol {v}}}{\partial t}}+({\boldsymbol {v}}\cdot \nabla ){\boldsymbol {v}}=-{\frac {1}{\rho }}\operatorname {grad} p+{\boldsymbol {g}}}$

となる。この式はオイラー方程式と呼ばれている。[14][15]

ポテンシャル流れ

渦度（速度場の回転）がない流れ

${\displaystyle {\boldsymbol {\omega }}=\operatorname {rot} {\boldsymbol {v}}=\mathbf {0} }$

の場合には、ベクトル解析の定理により

${\displaystyle {\boldsymbol {v}}=\operatorname {grad} \Phi }$

となる速度ポテンシャル Φ が存在する。

近似

ブシネスク近似

熱輸送を伴う流れにおいて、温度による密度変化が大きくないとして扱う近似法をブシネスク近似という。[16]

境界層近似

流れが主流方向を持ち（逆流、再循環および剥離がない）、幾何的な変形が緩やかなときに行う近似法を境界層近似という。

数値シミュレーション

もし一般解が求まれば、流体の挙動を完全に知る事ができることになるが、未だに一般解は発見されていない。また、解の存在可能性についても明らかとはなっておらず、物理学と数学の両方にまたがる重要な課題の一つとなっている[1][2]（詳細は、ミレニアム懸賞問題、ナビエ–ストークス方程式の解の存在と滑らかさを参照）。従って、極めて特殊な制約条件の問題を除いて数値解析によって近似的に解を求める[17][18]。

流体の数値シミュレーション（数値流体力学、CFD）では、このナビエ–ストークス方程式と連続の式、その他必要に応じてエネルギー保存の法則（熱対流）やマクスウェルの方程式（磁気流体力学）、状態方程式などを連立して、数値的に解くことで流体の挙動を予測する。[10][19][20][21]

移流と拡散両方に関係している現象であるので、クーラン数、拡散数の両方を満たすようにシミュレーションを行う必要がある。

性質

乱流

乱流は流体の多くの流れで見られる時間依存のカオス的な振る舞いである。[22][23][24][25]全体としての流体の慣性にそれがしたがうことが一般に信じられている。それゆえ慣性の効果が小さな流れは層流となる傾向がある。[26]移流と粘性の強さの比率はレイノルズ数と呼ばれる無次元量であり、レイノルズ数がある閾値を越えると、微小なかく乱が移流項の非線型性により拡大していき、流れ場は非定常な乱流となる[27]。

一方、右辺の粘性率を含む項（粘性項）は乱流の変動を抑制する効果を持つ。あまり深く理解されていないにもかかわらず、ナビエ‐ストークス方程式が乱流の性質を正確に記述することが信じられている[28][29]。計算に対して計算時間が有意味に解き得るようになるちょうどよい計算メッシュによる解のようなこの要求条件の安定した解または直接数値シミュレーションの、乱流に関するナビエ‐ストークス方程式の数値解は、極度に困難である。[30][31][32][33]難易度はその乱流に含まれている混合長さの尺度の違いに強く依存する。適当に変換するのに役立たない、層流を解くものを用いて乱流の流れを解く試みは、非定常解で典型的な結果を残す。

これに反して、レイノルズ平均ナビエ-ストークス方程式（RANS）のような、乱流モデルを補った時間平均方程式は、乱流をモデル化するときに実用的な数値流体力学（CFD）の応用で用いられる。追加の方程式を加えてRANSを導く、Spalart-Allmaras乱流モデル（英語版）、[34]k‐ω乱流モデル（英語版）、[35]k‐ε乱流モデル（英語版）を含む幾つかのモデルは、Large eddyシミュレーション（英語版）（LES）[36][37][38][39]がこれらの方程式を数値的に解くように用いるようにもできる。RANSよりも計算時間と計算機メモリーの面で、これらのアプローチは電子計算機で行うには大変コストがかかる。しかし、それは陽的に大きな乱流の尺度を分解するのでより良い結果を生み出すのである。

脚注

1. 前述の通り、ナビエ・ストークス方程式は流体がニュートン流体であることを前提としているため、非ニュートン流体に対しては成立しない。
2. 見かけは1本だが、x成分, y成分, z成分に分割すると3本となる。
3. 単純化された方程式を上手く選べば、数値計算の負荷を小さくできるため、依然これらの近似方程式は重要である（Ferziger, Perić, 2003）。

参考文献

1. ^ Constantin, P., & Foias, C. (1988). Navier-stokes equations. University of Chicago Press.
2. ^ 小薗英雄. (2002). Navier-Stokes 方程式. 数学, 54(2), 178-202.
3. C. L. M. H. Navier, "Mémoire sur les lois du mouvement des fluides," Mémoires Acad. Roy. Sci. Inst. France, 6, pp.389-440 (1823)
4. G. G. Stokes, "On the Theories of the Internal Friction of Fluids in Motion, and of the Equilibrium and Motion of Elastic Solids," Trans. Camb. Phil. Soc., 8, pp.287-319(1845）original paper
5. 児玉良明. (1996). CFD 入門 (その 1)− NS 方程式の様々な形 とモデル方程式 一. 日本造船学会誌, (805).
6. ^ 藤田宏. (1962). Navier-Stokes 方程式の数学的プロフイル. 日本物理学会誌, 17(4), 260-264.
7. 『渦度』 - コトバンク
8. 寺沢寛一 編『自然科学者のための数学概論 応用編』岩波書店、1960年、640頁。ISBN 4-00-005481-3。 
9. 柴田良弘、久保隆徹：「非線形偏微分方程式」、朝倉書店、ISBN 978-4-254-11771-4 (2012年)
10. ^ Joel H. Ferziger; Milovan Perić 著、小林敏雄、谷口伸行、坪倉誠 訳『コンピュータによる流体力学』シュプリンガー・フェアラーク東京、2003年、12–15頁。ISBN 4-431-70842-1。 
11. Panton, R. L. (2013). Incompressible flow. John Wiley & Sons.
12. Pironneau, O. (1973). On optimum profiles in Stokes flow. Journal of Fluid Mechanics, 59(1), 117-128.
13. Pozrikidis, C. (2001). Interfacial dynamics for Stokes flow. Journal of Computational Physics, 169(2), 250-301.
14. Christodoulou, Demetrios (October 2007). “The Euler Equations of Compressible Fluid Flow”. Bulletin of the American Mathematical Society 44 (4): 581–602. doi:10.1090/S0273-0979-07-01181-0. http://www.ams.org/bull/2007-44-04/S0273-0979-07-01181-0/S0273-0979-07-01181-0.pdf. 
15. Euler, Leonhard (1757). “Principes généraux du mouvement des fluides”. Mémoires de l'académie des sciences de Berlin 11: 274–315. https://scholarlycommons.pacific.edu/euler-works/226. 
16. Zeytounian, R. K. (2003). Joseph Boussinesq and his approximation: a contemporary view. Comptes Rendus Mecanique, 331(8), 575-586.
17. Temam, R. (2001). Navier-Stokes equations: theory and numerical analysis (Vol. 343). American Mathematical Society.
18. Girault, V., & Raviart, P. A. (2012). Finite element methods for Navier-Stokes equations: theory and algorithms (Vol. 5). Springer Science & Business Media.
19. Anderson, John D. (1995). Computational Fluid Dynamics: The Basics With Applications. Science/Engineering/Math. McGraw-Hill Science. ISBN 978-0-07-001685-9.
20. Chung, T. J. (2010). Computational fluid dynamics. Cambridge University Press.
21. Wesseling, P. (2009). Principles of computational fluid dynamics. Springer Science & Business Media.
22. 『乱流』 - コトバンク
23. H. Tennekes、J. L. Lumley、藤原仁志、荒川忠一訳『乱流入門』東海大学出版会、1998年。ISBN 978-4-486-01440-9。
24. Lesieur, M. (2012). Turbulence in fluids (Vol. 40). Springer Science & Business Media.
25. Davidson, P. A. (2015). Turbulence: an introduction for scientists and engineers. Oxford University Press.
26. 『層流』 - コトバンク
27. 『レイノルズ数』 - コトバンク
28. Foias, C., Manley, O., Rosa, R., & Temam, R. (2001). Navier-Stokes equations and turbulence (Vol. 83). Cambridge University Press.
29. R. G. Lerner; G. L. Trigg (1991). Encyclopaedia of Physics (2nd ed.). VHC publishers. ISBN 0-89573-752-3 
30. 大宮司久明, 三宅裕, & 吉澤徴. (1998). 乱流の数値流体力学. 東京大学出版会.
31. 梶島, & 岳夫. (2014). 乱流の数値シミュレーション. 養賢堂.
32. Wilcox, D. C. (1998). Turbulence modeling for CFD (Vol. 2, pp. 103-217). La Canada, CA: DCW industries.
33. Chen, C. J. (1997). Fundamentals of turbulence modelling. CRC Press.
34. Spalart, P. R. and Allmaras, S. R., 1992, "A One-Equation Turbulence Model for Aerodynamic Flows" AIAA Paper 92-0439
35. Wilcox, D. C. (2008), Formulation of the k–ω Turbulence Model Revisited, 46, AIAA Journal, pp. 2823–2838, Bibcode: 2008AIAAJ..46.2823W, doi:10.2514/1.36541 
36. Piomelli, U. (1999). Large-eddy simulation: achievements and challenges. Progress in Aerospace Sciences, 35(4), 335-362.
37. Mason, P. J. (1994). Large‐eddy simulation: A critical review of the technique. Quarterly Journal of the Royal Meteorological Society, 120(515), 1-26.
38. Zhiyin, Y. (2015). Large-eddy simulation: Past, present and the future. Chinese journal of Aeronautics, 28(1), 11-24.
39. Sagaut, P. (2006). Large eddy simulation for incompressible flows: an introduction. Springer Science & Business Media.

関連項目

• ナビエ–ストークス方程式の解の存在と滑らかさ
• ミレニアム懸賞問題
• 流体力学
• バーガース方程式
• 移流拡散方程式
• gifted/ギフテッド（2017年のアメリカ映画。ナビエ–ストークス方程式を題材として扱っている）
• 偏微分方程式
• 微分方程式

外部リンク

• 京都大学 Navier-Stokes方程式の数理とその応用プロジェクト