[論文レビュー] Variation on the Kolmogorov Forcing: Asymptotic Dissipation Rate Driven by Harmonic Forcing

Bertrand Rollin, Yves Dubief+1 more|arXiv (Cornell University)|Mar 11, 2009

Fluid Dynamics and Turbulent Flows参考文献 22被引用数 15

ひとこと要約

本研究は、3次元周期的乱流における漸近的エネルギー散逸率に及ぼす調和的体力の形状の影響を調査する。直接数値シミュレーション（DNS）と変分最適化を用いて、次元なし散逸係数βが力の形状に強く依存することを示した。無限大のレイノルズ数における理論的上限と定性的には一致するが、定量的予測は限定的である。

ABSTRACT

The relation between the shape of the force driving a turbulent flow and the upper bound on the dimensionless dissipation factor $β$ is presented. We are interested in non-trivial (more than two wave numbers) forcing functions in a three dimensional domain periodic in all directions. A comparative analysis between results given by the optimization problem and the results of Direct Numerical Simulations is performed. We report that the bound on the dissipation factor in the case of infinite Reynolds numbers have the same qualitative behavior as for the dissipation factor at finite Reynolds number. As predicted by the analysis, the dissipation factor depends strongly on the force shape. However, the optimization problem does not predict accurately the quantitative behavior. We complete our study by analyzing the mean flow profile in relation to the Stokes flow profile and the optimal multiplier profile shape for different force-shapes. We observe that in our 3D-periodic domain, the mean velocity profile and the Stokes flow profile reproduce all the characteristic features of the force-shape. The optimal multiplier proves to be linked to the intensity of the wave numbers of the forcing function.

研究の動機と目的

• 3次元周期的乱流における、非自明で複数の波数を有する力関数が次元なし散逸係数βの上界に与える影響を検討すること。
• 中程度のレイノルズ数における直接数値シミュレーション（DNS）の結果と、変分最適化から導かれた解析的上限を比較すること。
• 平均速度プロファイル、ストークス流れプロファイル、および異なる力の形状における最適乗数プロファイルとの関係を調査すること。
• ドーリングら（2003）の最適化フレームワークが有限レイノルズ数流れにおける定量的散逸挙動を正確に予測できるかどうかを評価すること。

提案手法

• 調和的力関数 f(y) = sin(ky) + A·sin(3ky) を用いた、3次元周期領域における非圧縮性ナビエ-ストークス系の定式化。
• ドーリングとコスタンチンの手法を用いて、次元なし散逸係数 β = εl/U³ の上界を変分最適化により導出。
• 周期的境界条件を満たすナビエ-ストークス方程式の直接数値シミュレーション（DNS）を実施し、βおよび平均速度プロファイルを計算。
• DNSの結果と解析的上限を比較し、力の振幅および波数成分の組み合わせに伴うβおよび速度プロファイルの依存性に注目。
• 最適乗数プロファイル ψₘ を分析し、力関数および平均速度場との関係を理解すること。
• 正規化された力、平均速度、最適乗数プロファイルを用いて、異なる力の振幅および波数における構造的相関を可視化すること。

実験結果

リサーチクエスチョン

• RQ1 複数の高調波を有する体力の形状は、3次元周期的乱流における次元なし散逸係数βの上界にどのように影響するか？
• RQ2 無限大のレイノルズ数で導かれたβの解析的上限は、有限レイノルズ数におけるDNSにおいて、定性的には成立するか、定量的にも成立するか？
• RQ3 3次元周期領域における乱流の平均速度プロファイルは、力関数の構造的特徴（勾配の変化、局所的極値など）をどのように反映するか？
• RQ4 最適乗数プロファイル ψₘ と力関数、または平均速度プロファイルとの関係は何か？
• RQ5 力関数に高次波数成分を追加することで、βは予測可能な非線形増加を示すか？

主な発見

• 次元なし散逸係数βは力関数の形状に強く依存しており、高次波数成分はその振幅が十分に大きくない限り、散逸に寄与しにくい。
• 変分最適化により導かれたβの解析的上界は、無限大のレイノルズ数における力の形状依存性を定性的に捉えており、中程度のレイノルズ数におけるDNSでもその傾向が定性的に維持されている。
• 乱流における平均速度プロファイルは、力関数の特徴（勾配の変化、局所的極値など）をすべて再現しており、特に複数の波数を含む力の場合顕著である。
• 最適乗数プロファイル ψₘ は主に力関数内の支配的波数に敏感であり、二次的成分の振幅が小さい場合にはそれらの影響をほとんど受けない。
• 二次波数成分の振幅が増加すると、ψₘ は組み合わせた波数構造を反映するよう変化し、その振幅はゆっくり減少し、二次成分が支配的になると安定化する。
• ψₘ と平均速度プロファイルの関係は明確でない。ψₘ は平均速度プロファイルに存在する二次的特徴を反映しないことがあり、乗数と流れ応答との間で分離が生じていることを示唆している。