[論文レビュー] An exact axisymmetric spiral solution of incompressible 3D Euler equations

Liang Sun|arXiv (Cornell University)|Jan 29, 2011

Computational Fluid Dynamics and Aerodynamics参考文献 5被引用数 5

ひとこと要約

本稿では、分離変数法を用いて、3次元非圧縮性Euler方程式に対する一般解としての軸対称らせん解を提示する。これにより、多層ラジアル傘状渦、多放物面型渦といった、有界で物理的に意味のある渦構造が得られる。解はベルヌーイ面におけるらせん流線を捉え、台風、 tornados、メソスケール渦の明示的モデルを提供する。最大風速やらせん角といった観測可能な流れの特徴と一致する3つの調整可能なパラメータを備える。

ABSTRACT

Spiral structure is one of the most common structures in the nature flows. A general steady spiral solution of incompressible inviscid axisymmetric flow was obtained analytically by applying separation of variables to the 3D Euler equations. The solution, depending on 3 parameters, describes the spiral path of the fluid material element on the Bernoulli surface, whereas some new exact solutions were obtained to be bounded within the whole region. The first one is a continued typhoon-like vortex solution, where there are two intrinsic length scales. One is the radius of maximum circular velocity $r_m$, the other is the radius of the vortex kernel $r_k=\sqrt{2}r_m$. The second one is a multi-planar solution, periodically in $z$-coordinate. Within each layer, the solution is a umbrella vortex similar to the first one. The third one is also a multi-planar solution in $z$-coordinate. In each layer, it is a combination of two independent solutions like the Rankine vortex, which is also finite but discontinued for either vertical or horizontal velocity. The fourth one is a multi-paraboloid vortex solution finite for $z$-coordinate but infinite for $r$-coordinate. Besides, some classical simple solutions (Rankine vortex, Batchelor vortex, Hill spherical vortex, etc.) were also obtained. The above explicit solutions can be applied to study the radial structure of the typhoon, tornados and mesoscale eddies. Both the solutions and approaches used here could also be applied to other complex flows by the Navier-Stokes equations.

研究の動機と目的

• 軸対称的、定常的、無粘性な流れの条件下で、3次元非圧縮性Euler方程式の正確な解析的解を導出すること。
• 従来の解が無限大に発散するか、独立していたという限界を克服し、半径方向、方位角方向、軸方向の速度成分の間の結合を強制すること。
• 自然のらせん流れパターンを再現する、物理的に意味があり、有限で有界な渦解を構築すること。
• 分離変数法を用いた体系的な手法を提供し、観測可能な渦の特性（例：最大風速、螺旋角）と一致する調整可能なパラメータを備えた新しい正確な解を生成すること。

提案手法

• 円筒座標系 (r, θ, z) における3次元非圧縮性Euler方程式に分離変数法を適用し、軸対称性と定常流れを仮定する。
• 速度成分 (Vr, Vθ, Vz) を r と z の関数として表現し、Vθ が Ψ/r に比例する形に設定する。ここで Ψ はストリーム関数である。
• 渦度が速度と平行となる一般化されたベルトラミー流条件を導出し、非線形Euler方程式を解ける形に簡略化する。
• 半径関数 R(r) と軸方向関数 H(z) を導入し、ベルヌーイ関数 B = −λR(r)H(z) が流れの構造を定義する。
• 質量保存則と運動量方程式を用いて、R(r) に関する2階常微分方程式を導出し、特定のパラメータ制約下で解けることを示す。
• 内層と外層などの異なる半径領域における解を組み合わせ、界面半径で微分が一致する連続的で区分的解を構築する。

実験結果

リサーチクエスチョン

• RQ1 定常的、無粘性、軸対称な条件下で、3次元非圧縮性Euler方程式に対する一般解としての正確な軸対称らせん解を導出可能か？
• RQ2 有界で有限な渦解を生じる物理的制約とパラメータ依存性は何か？ これは、実際の気象・海洋渦に類似した解を生じる。
• RQ3 分離変数法を結合された非線形Euler方程式に適用することで、ベルヌーイ面上のらせん軌道を持つ明示的解を生成できるか？
• RQ4 新規に得られた渦解の構造的・運動的性質は何か？ これに対してランキン渦やバッチェロール渦といった古典的モデルと比較するとどうなるか？
• RQ5 得られた解を用いて、台風や tornados の半径方向の速度プロファイルと角運動量分布をモデル化・解釈できるか？

主な発見

• 3次元非圧縮性Euler方程式に対する一般解としての正確ならせん解が導出され、流体要素がベルヌーイ面上をらせん経路に沿って移動することを記述する。
• 3つの新しい有界渦解が得られた：2つの固有長尺度（rm と rk = √2 rm）を有する多平面的傘状渦、2つのランキン型渦を組み合わせた多平面的渦、z で有限だが r で無限大の多放物面型渦。
• 解はランキン渦、バッチェロール渦、ヒルの球形渦を特別な場合として含み、既知の古典的解と整合性を確認する。
• 角運動量は m = rVθ = −λR(r)H(z) で与えられ、R(r) の半径依存性により、最大風速が rm で発生するような観測渦構造のモデル化が可能である。
• 3つのパラメータ（rm：最大風速半径、Vθm：最大接線速度、k/λ：らせん角）を調整することで、観測データに適合可能であり、台風や tornados の流れ構造による分類が可能である。
• この手法は非定常な粘性流や、f-平面における地球流体力学的流れ、磁気力学的流れなどの他の複雑系へも拡張可能である。