[論文レビュー] The role of microscopic fluctuations in transition prediction

Paolo Luchini|ArXiv.org|Apr 13, 2008

Phase Equilibria and Thermodynamics参考文献 4被引用数 7

ひとこと要約

本稿は、フラクチュアーション・ドレスタン定理に起因する微視的熱揺らぎが、せん断流れにおける乱流転移を誘発する物理的に整合的で十分な初期摂動源であることを示している。この揺らぎは、経験的に観測された $e^N$ 転移基準 $N \approx 9\text{--}10$ ネーパーと一致する。確率的線形安定性理論を用いて、熱ノイズが境界層不安定性を破綻の閾値まで駆動できることを示し、転移予測モデルにおける外部摂動の必要性を排除した。

ABSTRACT

The commonly accepted description of transition to turbulence in shear flows requires the presence of an external source of disturbances that get amplified by an essentially linear mechanism up to the point where breakdown to turbulence occurs. Microscopic fluctuations are shown here to provide just the right amount of initial disturbances to match the predictions of linear stability theory.

研究の動機と目的

• 微視的熱揺らぎが、乱流転移に必要な初期摂動の物理的に根拠のある源として機能できるかどうかを調査すること。
• 外部摂動（例：音響、渇度）が仮定されているが現実的に指定されていない、長年の転移予測のギャップを埋めること。
• 熱ノイズがアンプ入力の根本的下限として機能する場合、経験的に観測された $e^N$ 転移閾値を再現できるかどうかを検証すること。
• 確率的線形システム理論を用いて、熱ノイズとせん断流れにおける乱流発生の間の理論的かつ定量的関係を確立すること。

提案手法

• 熱揺らぎを表す確率的応力テンソルを含むナビエ＝ストークス方程式に、確率的に駆動される線形系の形式的取り扱いを適用する。
• ラプラス変換と固有モード分解を用いて、白色ノイズ入力に対する系の応答を表現し、主な不安定モードに注目する。
• フラクチュエーション・ドレスタン定理を用いて、速度摂動の平均二乗（rms）振幅を導出する。ノイズのパワースペクトルは温度および流体の性質と結びつける。
• レイノルズ数 $Re$、基準長さ $h$、および新しい物性長さ $\lambda = k_B T / (\rho \nu^2)$ を用いた次元なしスケーリングを施し、応答の次元なし形式を得る。
• 準平行流れにおける空間的受動問題を定式化し、距離にわたる増幅度として $e^N$ 法と比較可能にする。
• 空間的受動問題の数値的統合により、熱ノイズのみで $u_{\text{rms}} \approx 1\%$ の自由流速度に達することが確認され、$N \approx 9\text{--}10$ で転移が発生する。

実験結果

リサーチクエスチョン

• RQ1 微視的熱揺らぎのみで、せん断流れにおける乱流転移を誘発する初期摂動を十分に説明できるか？
• RQ2 熱ノイズに対する境界層の応答が、経験的に観測された $e^N$ 転移基準 $N \approx 9\text{--}10$ ネーパーを再現するか？
• RQ3 熱ノイズ強度と境界層内におけるトールマイエン＝シュリヒト波の増幅度との間の定量的関係は何か？
• RQ4 流体の性質から導かれる次元なしパラメータ $\lambda/h$ が、転移閾値にどのように影響するか？
• RQ5 $e^N$ 法の経験的閾値が、線形アンプ系における熱ノイズが定める理論的下限と整合的か？

主な発見

• 熱ノイズによる rms 速度揺らぎは、$N \approx 9\text{--}10$ で $e^N$ 法が転移を予測する位置と同一の場所で、自由流速度の $1\%$ に達する。
• 理論的推定式 $N = 0.5 \log(h/\lambda)$ により、300 K の空気中で境界層厚さ $h = 1$ mm の場合、$N \approx 9.00$ が得られ、経験的閾値と正確に一致する。
• 標準空気（300 K、大気圧）における特徴的長さ $\lambda = 1.508 \times 10^{-11}$ m は、$k_B T / (\rho \nu^2)$ から導かれる物性量である。
• モデルは、境界層を含むいかなるアンプに対しても、入力ノイズの有限で越えられない下限を設定することを示しており、外部摂動の必要性を転移予測から排除する。
• 数値結果により、大気音や渇度が存在しない状況でも、熱揺らぎのみで $N \approx 9\text{--}10$ で転移が発生することが確認された。
• 応答の次元なし形式 $\left<|C^{(1)}|^2\right> \propto \frac{\lambda}{h} Re^{-3} \times \text{空間的積分}$ は、$\lambda/h \sim e^{-2N}$ のとき、$e^N$ 閾値が熱揺らぎから自然に出現することを示している。