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

提供: testwiki
2025年3月15日 (土) 02:18時点におけるimported>Wetchによる版 (→粘性率が一定の非圧縮性流れ)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
ナビゲーションに移動 検索に移動

テンプレート:連続体力学

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

導出

流体の質量保存の法則と運動量保存の法則を表す連続の方程式 テンプレート:Indent テンプレート:Indent を用いると、流れの速度場 𝒗 の物質微分は テンプレート:Indent と導かれる。ここで、 ρ は密度場、 σ は応力場、 𝒈 は流体の単位質量あたりに作用する外力場（加速度場）である。

ここで、ニュートン流体を仮定すれば、応力場が テンプレート:Indent で与えられる。ただし、 p は圧力（静圧）、 χ は体積粘性率、 μ は剪断粘性率である。 𝒆 は対称化したテンプレート:仮リンクで、デカルト座標の下で成分表示をすれば テンプレート:Indent で表され、 Θ は速度場の発散 テンプレート:Indent である。

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

なお、速度場の物質微分の第二項は「対流項」あるいは「移流項」と呼ばれ、ベクトル解析の公式により、 テンプレート:Indent と変形することができる。ここで テンプレート:Mvar は速度場の回転 テンプレート:Indent であり、渦度と呼ばれる[7]。

単純化した方程式

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

非圧縮性流れ

非圧縮性流れ[11]では、速度場の発散 テンプレート:Mvar がゼロなので、速度場の発散を含む項を落として

D𝒗Dt=∂𝒗∂t+(𝒗⋅∇)𝒗=−1ρgrad⁡p+μρΔ𝒗+1ρgrad⁡(𝒗⋅grad⁡μ)+1ρrot⁡(𝒗×grad⁡μ)−1ρ𝒗Δμ+𝒈

となる。

粘性率が一定の流れ

粘性率 テンプレート:Mvar や テンプレート:Mvar は温度や圧力の関数であり一定ではないが、粘性率を定数と仮定する場合は、粘性率の勾配を含む項を落として、

D𝒗Dt=∂𝒗∂t+(𝒗⋅∇)𝒗=−1ρgrad⁡p+μρΔ𝒗+χ+13μρgrad⁡Θ+𝒈

となる。また、体積粘性率 テンプレート:Mvar は非常に小さいので、テンプレート:Math と仮定すると、

D𝒗Dt=∂𝒗∂t+(𝒗⋅∇)𝒗=−1ρgrad⁡p+νΔ𝒗+ν3grad⁡Θ+𝒈

となる（ストークスの仮説）。ここで テンプレート:Math は動粘性率である。

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

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

D𝒗Dt=∂𝒗∂t+(𝒗⋅∇)𝒗=−1ρgrad⁡p+νΔ𝒗+𝒈

となる。ここで テンプレート:Math は動粘性率である。各項はそれぞれ、

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

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

ストークス流れ（クリープ流れ）
粘性率が一定の非圧縮性流れのうち、流体の速度が遅かったりスケールが小さいなど、レイノルズ数が小さい流れを特にストークス流れあるいはクリープ流れという。ストークス流れでは、非線型である対流項 (𝒗⋅∇)𝒗 が無視できて、
D𝒗Dt=∂𝒗∂t=−1ρgrad⁡p+νΔ𝒗+𝒈
となる。この式はストークス方程式（Stokes equations）と呼ばれ[12][13]、線形方程式のため基本解が知られている[14]。例としてクエット流れやハーゲン・ポアズイユ流れがある。

オイラー方程式

粘性のない（テンプレート:Math）流れでは

D𝒗Dt=∂𝒗∂t+(𝒗⋅∇)𝒗=−1ρgrad⁡p+𝒈

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

ポテンシャル流れ

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

ω=rot⁡𝒗=𝟎

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

𝒗=grad⁡Φ

となる速度ポテンシャル テンプレート:Mvar が存在する。

近似

ブシネスク近似
熱輸送を伴う流れにおいて、温度による密度変化が大きくないとして扱う近似法をブシネスク近似という。[17]
境界層近似
流れが主流方向を持ち（逆流、再循環および剥離がない）、幾何的な変形が緩やかなときに行う近似法を境界層近似という。

数値シミュレーション

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

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

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

性質

乱流

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

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

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

脚注

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

参考文献

テンプレート:Reflist

関連項目

外部リンク

テンプレート:Authority control

  1. ↑ 1.0 1.1 1.2 1.3 Constantin, P., & Foias, C. (1988). Navier-stokes equations. University of Chicago Press.
  2. ↑ 2.0 2.1 2.2 2.3 小薗英雄. (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. ↑ 6.0 6.1 藤田宏. (1962). Navier-Stokes 方程式の数学的プロフイル. 日本物理学会誌, 17(4), 260-264.
  7. ↑ テンプレート:Kotobank
  8. ↑ テンプレート:Cite
  9. ↑ 柴田良弘、久保隆徹：「非線形偏微分方程式」、朝倉書店、ISBN 978-4-254-11771-4 (2012年)
  10. ↑ 10.0 10.1 テンプレート:Cite
  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. ↑ テンプレート:Cite
  15. ↑ テンプレート:Cite journal
  16. ↑ テンプレート:Cite journal
  17. ↑ Zeytounian, R. K. (2003). Joseph Boussinesq and his approximation: a contemporary view. Comptes Rendus Mecanique, 331(8), 575-586.
  18. ↑ Temam, R. (2001). Navier-Stokes equations: theory and numerical analysis (Vol. 343). American Mathematical Society.
  19. ↑ Girault, V., & Raviart, P. A. (2012). Finite element methods for Navier-Stokes equations: theory and algorithms (Vol. 5). Springer Science & Business Media.
  20. ↑ Anderson, John D. (1995). Computational Fluid Dynamics: The Basics With Applications. Science/Engineering/Math. McGraw-Hill Science. ISBN 978-0-07-001685-9.
  21. ↑ Chung, T. J. (2010). Computational fluid dynamics. Cambridge University Press.
  22. ↑ Wesseling, P. (2009). Principles of computational fluid dynamics. Springer Science & Business Media.
  23. ↑ テンプレート:Kotobank
  24. ↑ H. Tennekes、J. L. Lumley、藤原仁志、荒川忠一訳『乱流入門』東海大学出版会、1998年。ISBN 978-4-486-01440-9。
  25. ↑ Lesieur, M. (2012). Turbulence in fluids (Vol. 40). Springer Science & Business Media.
  26. ↑ Davidson, P. A. (2015). Turbulence: an introduction for scientists and engineers. Oxford University Press.
  27. ↑ テンプレート:Kotobank
  28. ↑ テンプレート:Kotobank
  29. ↑ Foias, C., Manley, O., Rosa, R., & Temam, R. (2001). Navier-Stokes equations and turbulence (Vol. 83). Cambridge University Press.
  30. ↑ テンプレート:Cite book
  31. ↑ 大宮司久明, 三宅裕, & 吉澤徴. (1998). 乱流の数値流体力学. 東京大学出版会.
  32. ↑ 梶島, & 岳夫. (2014). 乱流の数値シミュレーション. 養賢堂.
  33. ↑ Wilcox, D. C. (1998). Turbulence modeling for CFD (Vol. 2, pp. 103-217). La Canada, CA: DCW industries.
  34. ↑ Chen, C. J. (1997). Fundamentals of turbulence modelling. CRC Press.
  35. ↑ Spalart, P. R. and Allmaras, S. R., 1992, "A One-Equation Turbulence Model for Aerodynamic Flows" AIAA Paper 92-0439
  36. ↑ テンプレート:Citation
  37. ↑ Piomelli, U. (1999). Large-eddy simulation: achievements and challenges. Progress in Aerospace Sciences, 35(4), 335-362.
  38. ↑ Mason, P. J. (1994). Large‐eddy simulation: A critical review of the technique. Quarterly Journal of the Royal Meteorological Society, 120(515), 1-26.
  39. ↑ Zhiyin, Y. (2015). Large-eddy simulation: Past, present and the future. Chinese journal of Aeronautics, 28(1), 11-24.
  40. ↑ Sagaut, P. (2006). Large eddy simulation for incompressible flows: an introduction. Springer Science & Business Media.
「https://ja.wiki.beta.math.wmflabs.org/w/index.php?title=ナビエ–ストークス方程式&oldid=836」から取得