[論文レビュー] Stratified inclined duct: direct numerical simulations

Lu Zhu, Amir Atoufi+6 more|arXiv (Cornell University)|Jan 24, 2023

Meteorological Phenomena and Simulations参考文献 31被引用数 2

ひとこと要約

本研究は、層密度が異なる傾斜管（SID）の三次元直接数値シミュレーション（DNS）を初めて実施した。リザーバーと非周期的境界条件を正確にモデル化することで、長時間にわたり高分解能で実験と一致する流れの状態（層流、波動、間欠的乱流、完全乱流）を再現し、ヒドラリック制御を駆動する低圧領域を特定した。その結果、混合効率係数Γ ≈ 0.1が得られ、古典的値0.2より低かった。

ABSTRACT

The stratified inclined duct (SID) experiment consists of a zero-net-volume exchange flow in a long tilted rectangular duct, which allows the study of realistic stratified shear flows with sustained internal forcing. We present the first three-dimensional direct numerical simulations (DNS) of SID to explore the transitions between increasingly turbulent flow regimes first described by Meyer \& Linden ( extit{J. Fluid Mech.} extbf{753}, 242-253, 2014). We develop a numerical set-up that faithfully reproduces the experiments and sustains the flow for arbitrarily long times at minimal computational cost. We recover the four qualitative flow regimes found experimentally in the same regions of parameter space: laminar flow, waves, intermittent turbulence, and fully-developed turbulence. We find good qualitative and quantitative agreement between DNS and experiments and highlight the added value of DNS to complement experimental diagnostics and increase our understanding of the transition to turbulence, both temporally (laminar/turbulent cycles) and parametrically (as the tilt angle of the duct and the Reynolds number are increased). These results demonstrate that numerical studies of SID -- and deeper integration between simulations and experiments -- have the potential to lead to a better understanding of stratified turbulence in environmental flows.

研究の動機と目的

• 層密度が異なる傾斜管（SID）実験の物理的境界条件を忠実に再現できる計算効率の高いDNSフレームワークの開発。
• 非周期的流れ方向の力およびリザーバー効果を含む、現実的な境界条件のもとで長期間にわたる三次元層密度乱流のシミュレーション。
• SIDにおける層流から乱流への遷移、特にヒドラリック制御と圧力勾配の役割の解明。
• 実験では得られない高分解能の体積平均データ（乱流運動エネルギー、圧力場、エネルギーフラックス）の提供。

提案手法

• 局所的力による簡略化されたリザーバー・モデルを用いた新規数値設定を開発し、計算コストを最小限に抑えながら長期間にわたる流れの維持を実現。
• 流れ方向に非周期的境界条件を導入し、物理的リザーバーと傾斜に起因する圧力勾配を再現。
• Prandtl数Pr = 7の条件下で、非圧縮性ナビエ–ストークス方程式と移流拡散方程式を解くためにXcompact3dスペクトル要素コードを採用。
• 乱流運動エネルギー（TKE）の分析に移動平均法を適用し、空間的・時間的間欠性のパターンを同定。
• 簡略化された力のアプローチの正確性と頑健性を確認するため、完全に解像されたリザーバーのシミュレーションとDNSを比較検証。
• 流れの安定性と混合のメカニズムを解析するため、勾配リチャードソン数、乱流エネルギーフラックス、圧力場といった主要な診断量を計算。

実験結果

リサーチクエスチョン

• RQ1 直接数値シミュレーション（DNS）は、SID実験で観測された4つの異なる流れ状態（層流、界面波、間欠的乱流、完全乱流）をどのように忠実に再現できるか？
• RQ2 流れ方向の圧力勾配とリザーバーに起因する力は、傾斜管における乱流の維持とヒドラリック制御の実現にどのような役割を果たすか？
• RQ3 乱流運動エネルギー、圧力場、エネルギーフラックスは管内をどのように変化するか。それらは混合メカニズムに何を明らかにするか？
• RQ4 実験的手法では到達できないスケールのダイナミクス（例えば、コherent構造や間欠性）をDNSはどの程度解像できるか？
• RQ5 DNSは混合効率、特にフラックス係数Γをどの程度正確に評価できるか。古典的推定値と比較してどう異なるか？

主な発見

• DNSは、同一のθ–Reパラメータ空間内で、実験で観測されたすべての4つの流れ状態（層流、界面波、間欠的乱流、完全乱流）を再現できた。
• 乱流状態では勾配リチャードソン数がRi_g ≈ 0.10–0.15に安定化し、Prandtl数の違いがあるにもかかわらず実験観測と一致した。
• 非層流状態のすべての流れで、管内に大規模な低圧領域が形成され、管の前半では有利な圧力勾配、後半では逆圧力勾配が生じ、ヒドラリック制御の新たなメカニズムを示唆した。
• 乱流運動エネルギーは空間的・時間的間欠性を示し、流れの速度に同期した再層流化と乱流バーストが管に沿って伝播した。
• 実験から導出されたエネルギーフラックスモデル（5.6）と経験的関係式（5.10）がDNSで裏付けられ、高乱流状態では混合効率係数がΓ ≈ 0.1であることが確認された。
• DNSの結果から、この系では古典的混合効率値0.2が過大評価されていることが判明し、海洋混合モデルにおける仮定の再評価が求められることが示された。