[論文レビュー] Identifying the active flow regions that drive linear and nonlinear instabilities

Olivier Marquet, Lutz Lesshafft|arXiv (Cornell University)|Aug 30, 2015

Nonlinear Dynamics and Pattern Formation参考文献 16被引用数 4

ひとこと要約

本論文は、振動流における線形および非線形不安定性を駆動する特定の流れ領域を同定するための新しいフレームワーク「エンドゲニティ（endogeneity）」を導入する。時間微分された直接流れ状態と随伴モードの局所的内積を計算することで、各空間的点がグローバル固有周波数および増幅率に与える寄与度を定量化し、ギンツブルク＝ランダウ方程式や円柱のうねり流れといった代表的な流れにおける不安定性メカニズムの詳細な物理的解釈を可能にする。

ABSTRACT

A new framework for the analysis of unstable oscillator flows is explored. In linear settings, temporally growing perturbations in a non-parallel flow represent unstable eigenmodes of the linear flow operator. In nonlinear settings, self-sustained periodic oscillations of finite amplitude are commonly described as nonlinear global modes. In both cases the flow dynamics may be qualified as being endogenous, as opposed to the exogenous behaviour of amplifier flows driven by external forcing. This paper introduces the endogeneity concept, a specific definition of the sensitivity of the global frequency and growth rate with respect to variations of the flow operator. The endogeneity, defined both in linear and nonlinear settings, characterizes the contribution of localized flow regions to the global eigendynamics. It is calculated in a simple manner as the local point-wise inner product between the time derivative of the direct flow state and an adjoint mode. This study demonstrates for two canonical examples, the Ginzburg-Landau equation and the wake of a circular cylinder, how an analysis based on the endogeneity may be used for a physical discussion of the mechanisms that drive a global instability. The results are shown to be consistent with earlier 'wavemaker' definitions found in the literature, but the present formalism enables a more detailed discussion: a clear distinction is made between oscillation frequency and growth rate, and individual contributions from the various terms of the flow operator can be isolated and separately discussed. [Truncated]

研究の動機と目的

• 非平行流れにおけるグローバル不安定性の空間的起源を物理的に解釈可能なフレームワークとして開発すること。
• 線形および非線形不安定性メカニズムにおける周波数と増幅率への寄与を区別すること。
• 平均流れや摂動成分といった流れ作用素内の個々の項の役割を分離・定量化する感度ベースの手法を提供すること。
• 時間平均された平均流れの線形安定性解析が、円柱のうねり流れなど特定の流れにおいて非線形飽和振動数を正確に予測できる理由を明らかにすること。
• 随伴モードと離散固有値問題を用いて、局所的寄与を安定かつ座標不変に計算する手法を確立すること。

提案手法

• エンドゲニティを、直接流れ状態の時間微分と随伴モードの局所的点ごとの内積として定義し、不安定性寄与の空間的マップを生成する。
• 離散随伴固有値問題を用いてエンドゲニティを定式化し、線形化作用素の転置共役を用いて随伴固有ベクトルを解く。
• 非一様メッシュや座標系（例えば円柱座標）に対応するため、重み行列 Q を用いた一般化内積を導入する。
• 内積の定義に依存しないエンドゲニティの不変性を示し、結果が計量行列 Q に依存せず、随伴固有ベクトルにのみ依存することを確認する。
• 線形ギンツブルグ＝ランダウ方程式と非線形円柱のうねり流れの2つの代表的系にこの手法を適用する。
• 非線形流れ作用素を時間平均成分と摂動成分に分解し、準線形ダイナミクスと高調波相互作用を別々に分析する。

実験結果

リサーチクエスチョン

• RQ1 流れ場内のどの局所的領域が、周波数および増幅率の観点からグローバル不安定性を主に駆動しているか？
• RQ2 流れ作用素内の異なる項（例：平均流れ対摂動）が不安定性に与える寄与をどのように分離・定量化できるか？
• RQ3 なぜ時間平均された平均流れの線形安定性解析が、円柱のうねり流れにおいて非線形飽和振動数を正確に予測できるのか？
• RQ4 エンドゲニティフレームワークは、より詳細かつ物理的に解釈可能な形で、古典的な「波源（wavemaker）」概念をどの程度回復または精緻化できるか？
• RQ5 エンドゲニティ測定値は、異なる座標系やメッシュ離散化に対しても安定的かつ不変に計算可能か？

主な発見

• エンドゲニティフレームワークは、不安定性を駆動する空間的領域を効果的に同定しており、古典的な波源概念と整合的であるが、物理的明確性が向上している。
• この手法は、固有値の実部（増幅率）と虚部（周波数）への寄与を明確に分離でき、増幅と振動メカニズムの別々の物理的解釈を可能にする。
• 円柱のうねり流れにおいては、非線形ダイナミクスが主に平均流れ上の準線形相互作用によって支配されており、高調波相互作用は二次的な役割を果たしていることが明らかになった。
• 時間平均された平均流れの線形安定性解析が非線形飽和周波数を正確に予測できるのは、主要な不安定性源が平均流れの局所的不安定性特性に起因しているためである。
• エンドゲニティは、内積定義の変更に対しても不変であり、特に円柱座標系においても有効であることが確認され、実用的応用における堅牢性が裏付けられた。
• 数値ソルバーへの実装を簡素化するために、計量行列 Q を必要とせず、随伴固有ベクトルから直接エンドゲニティを計算可能である。