[論文レビュー] Some closed form solutions to the Navier Stokes equations

R. M. Kiehn|ArXiv.org|Feb 1, 2001

Rheology and Fluid Dynamics Studies参考文献 1被引用数 7

ひとこと要約

本稿では、ナビエ＝ストークス方程式を積分因子を伴う Pfaffian 系に再定式化することで、コhomological 手法を用いて閉形式解を生成する新しいアプローチを提示する。回転対称性を有する速度場に対称性制約を課すことにより、系は積分可能となる2次元 Pfaffian 形式に簡略化され、間歇的現象、Hopf 分岐、流れの逆転、ヒステリシスを示す正確解が得られる。特に、Langford の鞍節点および透過的 Hopf 系統と関連している。

ABSTRACT

An algorithm for generating a class of closed form solutions to the Navier-Stokes equations is suggested, with examples. Of particular interest are those exact solutions that exhibit intermittency, tertiary Hopf bifurcations, flow reversal, and hysteresis.

研究の動機と目的

• 微分トポロジーとコhomology を用いて、ナビエ＝ストークス方程式に対する系統的アルゴリズムを構築し、正確な閉形式解を生成すること。
• 有効力場 f が Frobenius 積分可能性条件を満たすように速度場成分に課すべき制約条件を特定し、2関数表現を可能にする。
• 可変密度を積分因子として導入することにより、回転フレームや圧縮性流れへのこの手法の拡張を図ること。
• 正確なナビエ＝ストークス解と既知の力学系（特に Langford 系統）を結びつけることにより、三次分岐や間歇的現象を含む複雑な流体挙動をモデル化すること。
• 積分因子と Pfaffian 形式が特定されれば、非圧縮性または粘性領域においても、圧力が四則演算により計算可能であることを示すこと。

提案手法

• ナビエ＝ストークス方程式をコhomological 表記に再定式化する：f_V - f_D = -dP/ρ ここで f_V と f_D は、対流項と粘性項を表す不完全1次微分形式である。
• 時間依存性が方位成分にのみ存在するものと仮定し、円柱座標系における速度場 V を成分 Λ(r,z)、Θ(r,z)、Φ(r,z,t)/r として定義する。
• Pfaffian 形式 W = f·dr を2次元に簡略化するため、f の方位成分が0であるという制約 γ(r,z,t) = 0 を課す。これにより積分可能性が保証される。
• 2次元 Pfaffian 形式 α(r,z)dz + β(r,z)dr = -dP/ρ に積分因子が存在することを用い、これを完全形式にすることで、圧力 P を積分により得られるようにする。
• コリオリ力および遠心力項を含めた回転フレームへの適用により、回転ナビエ＝ストークス方程式と整合性を保つ。
• f の成分から導かれる常微分方程式を、Langford の鞍節点および透過的 Hopf 系統と一致させることで、解を既知の力学系にマッピングする。

実験結果

リサーチクエスチョン

• RQ1 コhomological フレームワークを用いて、ナビエ＝ストークス方程式の閉形式解を系統的に生成することは可能か？
• RQ2 速度場成分にどのような制約を課すと、有効力場 f が Frobenius 積分可能性を満たし、2関数表現が可能になるか？
• RQ3 仕事1形式の Pfaffian 形式を2次元系に簡略化し、積分因子を用いて積分可能とすることで、正確な圧力解を得られるか？
• RQ4 正確なナビエ＝ストークス解は、間歇的現象、流れの逆転、ヒステリシスといった複雑な流体現象をどの程度正確にモデル化できるか？
• RQ5 可変流体密度を積分因子として導入することで、この手法を圧縮性流れに拡張できるか？

主な発見

• 制約条件 γ(r,z,t) = 0 を用いることで、3次元ナビエ＝ストークス系を2次元 Pfaffian 形式に簡略化し、積分可能性を保証する。これにより正確解が得られる。
• ポアズイユ型解（a = b = B = 0）では、圧力は P = b(dr⁴/2 - (r² - 2z²)U - bz⁴)/2 - ν(4d - 2b)z で与えられ、ベクトルラプラシアンが非ゼロのとき粘性依存性を示す。
• 調和的状態（2d - b = 0）では、圧力は ν に依存せず、非散乱状態を示し、粘性流れの極限集合に対応する。
• 透過的 Hopf 分岐解では、P = -1/4cbr⁴ - 1/8A²r² + 1/2A²z² - Az³c + 1/2c²z⁴ - ν(4b - 2c)z となり、圧力が粘性および非粘性成分に分解される。
• ヒステリシス的 Hopf 解では、z = 2d/3A でベクトルラプラシアンの符号が反転し、散乱挙動の転移を示し、非グローバル等体積性のため間歇的挙動を示す。
• すべての解は Langford 系統と関連しており、dz/dt = Λ(r,z), dr/dt = Θ(r,z), dθ/dt = ω と表され、既知の分岐モデルと力学的に同等であることが確認される。