Leveraging unstructured grids for direct numerical simulations of wall turbulence

本論文は、壁乱流の直接数値シミュレーションにおける局所コルモゴロフスケールに合わせてグリッドサイズをスケーリングする非構造化グリッド生成フレームワークである{\eta}-グリッドを導入し、従来の直交グリッドと同等の精度を達成しつつ、特に高レイノルズ数および複雑なリブレット幾何学において計算コストを大幅に削減することを示す。

要約

自動車の上を流れる空気や船の船体を通過する水流をシミュレーションすることを想像してみてください。これをコンピュータ上で正確に行うために、科学者たちは**直接数値シミュレーション（DNS）**と呼ばれる手法を用います。DNS を考える際、流体（空気または水）を数百万個の微小で不可視な立方体（グリッド）に分割する、巨大な 3 次元デジタル顕微鏡を作成しているとイメージしてください。コンピュータは、それぞれの立方体がどのように移動し、隣接する立方体と相互作用するかを計算します。

問題は、表面（船の側面など）に近い流体の流れが、極めてカオス的で詳細に富んでいることです。明確な像を得るためには、表面のすぐ近くにこれらの微小な立方体を大量に配置する必要があります。しかし、表面から離れるにつれてカオスは滑らかになり、それほど微小な立方体は不要になります。

従来の方法：「剛性のレンガの壁」

従来、科学者たちは直交格子を用いていました。これは、同一で剛性のレンガを積み上げて壁を作るようなものです。

• 問題点： 表面近くの微細な詳細を捉えるためには、壁の下部にあるレンガを非常に小さくする必要があります。しかし、これらのレンガは剛性であり、一直線に連結されているため、詳細が重要でない上部に至るまで、同じ微小なレンガを使用せざるを得ません。
• 結果： 結果として、不必要なものが大半を占める、数十億個の微小なレンガで構成された壁ができてしまいます。これにより、コンピュータシミュレーションは信じられないほど遅く、高価なものになります。まるで、潮の満ち引きを測るために、浜辺の砂粒一つ一つを数えようとするようなものです。

新しい解決策：「賢く伸縮する網」

この論文では、**$\eta$-グリッド（イータ・グリッド）**と呼ばれる新しい手法が紹介されています。剛性のレンガの代わりに、賢く伸縮する漁網を想像してください。

• 仕組み： 著者たちは、網の穴のサイズが、必要な詳細度に応じて自動的に変化するシステムを設計しました。
  • 表面付近（「内層」）： 流体の微小でカオス的な渦を捉えるため、網の穴は非常に小さく、きつくなっています。
  • より遠方（「外層」）： 流体が穏やかになるにつれて、網は自動的に伸び、穴ははるかに大きくなります。
• 秘密の要素： これらの穴のサイズは、コルモゴロフスケール（$\eta$と表記）と呼ばれるものに基づいています。これは、任意の高さにおいて流体中に存在し得る「最小の渦」と考えてください。新しいグリッドは単に、「この特定の高さにおける最小の渦を捉えるのに十分な大きさの穴にし、それ以上大きくしない」というルールに従います。

これが画期的な理由

著者たちは、この「賢い網」を 2 種類の異なるコンピュータコード（一つはスペクトル要素法に類似、もう一つは有限体積法に類似）でテストし、従来の「剛性のレンガ」方式と比較しました。

1. 精度： 結果はほぼ同一でした。「賢い網」は、摩擦や速度などの主要な測定値において、1% 未満の差異しか見せず、「剛性のレンガ」と同様に物理現象を捉えました。
2. 莫大な節約： ここに魔法が宿ります。
   • 滑らかな表面（平坦な壁など）の場合、新しいグリッドは高速域において必要な「レンガ」（グリッド点）の数を約**90%**削減しました。
   • 粗い表面（抵抗低減のために設計された「リブレット」と呼ばれる微小な溝のある壁など）の場合、節約効果はさらに劇的で、グリッド点数が97% 減少しました。

「溝のある壁」の比喩

リブレットの部分を理解するために、微小な平行な溝で覆われた壁（ゴルフボールの質感やサメの皮膚のようなもの）を想像してください。

• 従来の方法： これをシミュレーションする場合、剛性のレンガ方式は、溝がグリッドを全体にわたって微細化することを強いるため、至る所でレンガを微小に保たなければなりませんでした。まるで、布地から遠く離れた部分を含め、セーターの糸一本一本を数えようとするようなものです。
• 新しい方法： 「賢い網」は、これらの微小な溝から数インチ離れると、流れが再び滑らかになることを知っています。そのため、溝の直上の穴を即座に伸ばし、もはや重要ではない微細な詳細を無視します。

結論

著者たちは、流体シミュレーションのための賢いズームレンズとして機能するフレームワークを構築しました。これは、必要な場所（壁の近く）に計算資源を集中させ、不要な場所では緩やかにするのです。

• 滑らかな壁の場合： シミュレーションが大型化するにつれて、必要な計算努力の増加が非常に緩やかになります。
• 粗い壁の場合： さらに優れたスケーリングを実現し、以前は処理できなかった複雑な抵抗低減表面のシミュレーションを、コンピュータ上で可能にします。

要するに、彼らは同じ高品質な作業を、計算能力のほんの一部で達成する方法を見出し、スーパーコンピュータで 1 ヶ月かかるかもしれないタスクを、数日で完了可能なものへと変えました。

技術概要

技術サマリー：壁面乱流の直接数値シミュレーションにおける非構造化格子の活用

問題提起
壁面乱流の直接数値シミュレーション（DNS）は、主に全領域にわたってコルモゴロフ長さスケール（$\eta$）を解像する必要があるため、計算コストが非常に高い。滑らかな面および非滑らかな面（浸没境界法などを用いた場合）の両方に広く用いられる従来の直交格子ソルバーは、通常、全領域にわたって流方向（$\Delta x^+$）および横方向（$\Delta z^+$）の格子解像度を固定、あるいは緩やかに変化させるように強制する。このアプローチは、コルモゴロフスケールが成長する外層領域において過剰な解像をもたらし、計算上のボトルネックを生じさせる。これは、リブルートのような複雑な表面を扱うシミュレーションにおいて特に顕著である。なぜなら、リブルートのような表面では、微小な溝を解像するために表面近傍で極めて微細な格子が必要となるが、直交格子の枠組みでは、この微細な解像度が境界層厚全体に不必要に維持されるからである。

手法
著者らは、$\eta$-grid と呼ばれる新しい非構造化格子生成フレームワークを提案する。その中核となる原理は、壁面法線方向（$\Delta y^+$）および横方向（$\Delta z^+$）の格子サイズを局所的なコルモゴロフスケール（$\eta^+$）に比例させることで、乱流スケールがより大きい外層領域において格子を粗化することを可能にするものである。

1. 格子定式化:

   • 内層: 壁に隣接する薄い層（約 50 粘性単位）は、従来の DNS と同様の粘性スケーリング格子サイズ（$0.3 \lesssim \Delta y^+ \lesssim 4$、$\Delta z^+ \simeq 5$）を使用する。非滑らかな表面（例：リブルート）の場合、格子は最小の表面波長（$\ell^+$）を$\ell^+/30 \lesssim \Delta y^+, \Delta z^+ \lesssim 4$となるように解像する。
   • 外層: 内層より上では、格子サイズは局所的なコルモゴロフスケールに比例して増加する（$\Delta y^+ \simeq \Delta z^+ \simeq 2\eta^+$）。
   • スケーリング則: 著者らは、乱流チャネル流およびゼロ圧力勾配（ZPG）乱流境界層（TBL）の対数領域およびウェイク領域における$\eta^+$の半経験的フィッティングを改良した。彼らは、$y^+ \approx \delta^+/2$で遷移する$\eta^+$のための区分的べき乗則フィッティングを提案する。
   • リブルート適応: リブルート幾何学に対して、このフレームワークは二次流およびケルビン - ヘルムホルツローラーを解像するために正方形セルで満たされた「リブルート亜層（$\delta_r^+$）」を導入し、その上で$\eta$に基づくスケーリングへ滑らかに遷移させる。

2. 実装:

   • このフレームワークは、マルチブロック非構造化格子を生成するためにオープンソースパッケージGmshを利用する。
   • 格子は、SOD2D（スペクトル要素法）およびOpenFOAM（有限体積法）という 2 つの異なるソルバーを用いてテストされる。
   • シミュレーションは、滑らかな壁および様々なリブルート幾何学（三角形の溝）における乱流チャネル流および ZPG TBL を、摩擦レイノルズ数$\delta_0^+ = 1000$までカバーする。

主要な貢献

• フレームワーク開発: 滑らかな面および複雑な粗面両方に適用可能であり、格子間隔を局所的なコルモゴロフスケールに動的に適応させる、堅牢なマルチブロック非構造化格子生成戦略。
• 改良されたスケーリング則: TBL およびチャネル流のウェイク領域におけるコルモゴロフスケールに対する更新された半経験的フィッティングにより、外層における格子サイズの精度が向上。
• ソルバー非依存性: SEM および FVM という異なる数値離散化スキーム全体で、$\eta$-grid フレームワークが一貫した高忠実度の結果をもたらすことの証明。
• リブルート固有の定式化: 亜層の物理（二次流、KH ローラー）を解像しつつ、外層での急速な粗化を可能にするリブルートに特化した格子戦略。これにより、直交格子/IBM アプローチの限界を克服する。

結果

• 精度: $\eta$-grid は、従来の直交格子および参照となる DNS/実験データと比較して 1% 未満の差で結果を生成する。指標には、摩擦係数（$C_f$）、平均速度プロファイル（$U^+$）、乱流応力、およびスペクトログラムが含まれる。
• 格子効率:
  • 滑らかな壁: 格子点数（$N_\eta$）は$\propto \delta_0^{+2.5}$としてスケーリングされ、双曲線正接伸縮を用いた直交格子の$\propto \delta_0^{+3.0}$と比較される。$\delta_0^+ = 6000$において、$\eta$-grid は標準的な直交格子の格子点数の約 10% しか必要としない（$N_\eta/N_{Tanh} \simeq 0.1$）。
  • リブルート: 効率の向上はさらに顕著である。抗力低減型の三角形リブルートの場合、$N_\eta$は$\propto \delta_0^{+2.0}$としてスケーリングされるのに対し、直交格子は$\propto \delta_0^{+3.0}$としてスケーリングされる。$\delta_0^+ = 6000$において、$\eta$-grid は直交格子の格子点数の約 3% しか必要としない（$N_\eta/N_{Tanh} \simeq 0.03$）。
• リブルート物理: このフレームワークは、仮想原点、二次流、およびケルビン - ヘルムホルツローラーを含むリブルート上の複雑な流れ機構を成功裡に捉え、中程度の格子パラメータ（$C_y=2.0, C_z=2.5, y_{in}^+=50$）を用いて格子収束を達成した。
• TBL 検証: リブルート上の ZPG TBL の DNS は、ウェイク領域における固有の実験-DNS 不一致にもかかわらず、抗力低減率および速度プロファイルに関して実験データ（Baron & Quadrio; Choi & Orchard）と良好な一致を示した。

重要性
本論文は、$\eta$-grid フレームワークが「莫大な格子節約」を提供し、複雑な表面における高レイノルズ数壁面乱流の DNS を計算上実行可能にすると主張している。固定された解像度を維持するのではなく、格子解像度を局所的な物理スケール（コルモゴロフスケール）に合わせることで、この手法は精度を犠牲にすることなく計算コスト（格子点数および時間ステップの制約）を劇的に削減する。著者らは、このアプローチにより、以前は直交格子ベースの手法の prohibitively 高いコストによって阻害されていた、実験条件に一致するレイノルズ数における複雑な表面の TBL の DNS を実行することが経済的に可能になると強調している。この研究は、構造化された直交メッシュの制約を超えて、高忠実度乱流シミュレーションに対する非構造化格子の有効性を検証するものである。