[論文レビュー] Time dependent simulation of the Driven Lid Cavity at High Reynolds Number

Nuno Vieira Cardoso, Pedro Bicudo|ArXiv.org|Sep 18, 2008

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

ひとこと要約

本研究では、1024×1024グリッドとSSPRK(5,4)スキームを用いた2次元駆動渦巣流れの時間依存シミュレーションを、Reynolds数30,000まで実施した。GPUで加速された倍精度計算により、結果は定常解に収束し、先行研究で報告された非定常解が粗いグリッドに起因する数値的アーティファクトであり、物理的挙動ではないことを解明した。

ABSTRACT

In this work, numerical solutions of the two dimensional time dependent incompressible flow, in a driven cavity at high Reynolds number Re, are presented. At high Re, there is a controversy. Some studies predicted that the flow is steady, others found time dependent non-steady flow, either periodic or aperiodic. In this study, the driven lid cavity is successfully solved using a very fine grid mesh, for Re up to $30 000$. We discretize the Vorticity-Stream formulation of the Navier-Stokes equation with the SSPRK(5,4) scheme in a $1024 imes1024$ grid. Using this very fine grid, the results obtained converge to a stationary solution. Detailed results for Re between $5 000$ and $30 000$ are presented. The driven lid cavity problem is solved with a NVIDIA GPU using the CUDA programming environment with double precision.

研究の動機と目的

• 高Reynolds数における駆動渦巣流れが非定常的になるかどうかという長年の論争を解決すること。
• 非常に細かいグリッドを用いた時間依存シミュレーションが、定常解に収束するか、それとも周期的・非周期的挙動を示すかを特定すること。
• 文献で報告された非定常解が、空間分解能が不十分なために生じる数値的アーティファクトであるかどうかを調査すること。
• 十分に細かいグリッドが真の物理的解を捉えるために不可欠であり、粗いメッシュによる誤った振動を回避できることを検証すること。
• GPUで加速された倍精度計算を用いた高分解能直接数値シミュレーションの有効性を示すこと。

提案手法

• 連続の自動満たしを実現するため、流れ関数-渦度形式を用いた2次元非圧縮性ナビエ＝ストークス方程式の数値的解法。
• 時間精度と安定性を確保するため、時間依存渦度輸送方程式をSSPRK(5,4)ルンゲ＝クッタスキームで離散化。
• 特にコーナー領域における小さな渦の形成に注目し、高い空間分解能を確保するため、1024×1024の構造的直交グリッドを採用。
• 倍精度算術を用いて、丸め誤差を最小限に抑え、解の忠実性を向上させるために、NVIDIA GPU上でCUDAを実装。
• 境界条件として、上部のリードレートを除き、すべての壁面にスリップなし条件を適用。
• 残差ノルムと渦度場の時間発展をモニタリングし、解が定常状態に達するまで時間積分を継続。

実験結果

リサーチクエスチョン

• RQ1 Reynolds数30,000に達する高Reynolds数における駆動渦巣流れは、時間依存的挙動を示すのか、それとも定常状態に収束するのか。
• RQ2 先行の直接数値シミュレーション（DNS）で報告された非定常解は、物理的不安定性に起因するものか、それとも粗いグリッド解像度に起因する数値的アーティファクトなのか。
• RQ3 非常に細かいグリッド（1024×1024）を用いることで、主な渦構造、特に小さなコーナー渦（BL3, TL2など）をすべて解像でき、解が安定化するか。
• RQ4 GPU上で倍精度算術を用いることで、低精度手法に比べてシミュレーションの精度と収束性が向上するか。
• RQ5 定常ナビエ＝ストークス方程式における定常解の存在が、物理的現実性を示すのに十分か、それとも時間依存シミュレーションが必要か。

主な発見

• Reynolds数5,000から30,000の範囲で、1024×1024グリッド上での時間依存シミュレーションは定常解に収束し、真の物理的解が定常的であることが示された。
• Re = 5,000では、主渦と一部のコーナー渦（BL1, BR1, TL1, TR1）が明確に可視化されるが、BL3とTL2の渦は解像度不足により欠落している。
• Re = 15,000に達すると、BL3とTL2を含むすべての主要渦が完全に解像され、解は安定している。
• Re = 20,000, 25,000, 30,000でも解は定常のまま維持され、高Reynolds数でも定常解に収束することが確認された。
• 本研究では、先行研究で報告された非定常または周期的流れは、粗いグリッドに起因する数値的振動に起因すると結論づけた。
• 結果は、Erturkの主張である「非常に細かいグリッドが、小スケール渦を解像し、正しい定常解を得るために不可欠である」という主張を支持するものである。