A generalized vertical coordinate transformation based on SPH(2) for efficient free surface flow simulations

技術要約：効率的な自由表面流シミュレーションのための一般化垂直座標変換

問題提起
平滑化粒子法（SPH）や移動粒子半隠式法（MPS）などの粒子法は、大きな変形や移動する不連続性を伴う自由表面流のシミュレーションにおいて効果的である。しかし、これらは大きな計算コストの課題に直面している。従来の粒子法では、精度を維持するために全領域にわたって一様な粒子間隔を必要とするため、（津波シミュレーションにおける深海のような）高解像度が重要ではない領域においても大量の粒子を必要とする。一方で、格子法はグリッドサイズを調整してコストを削減できるが、従来の粒子法にはこの柔軟性が欠けている。さらに、複雑な底面境界（海底や河床など）の扱いは、複雑な粒子の配置や境界処理を必要とすることが多く、それが精度や安定性の低下を招くことがある。著者らは、粒子法のラグランジュ的な利点を損なうことなく、非一様な粒子解像度（特に垂直方向）を可能にし、境界処理を簡素化することで、計算効率を向上させる手法の必要性を指摘している。

手法
本論文では、SPH法に適用される**一般化垂直座標変換（VCT）**に基づくフレームワークを提案している。核心となる概念は、物理空間（複雑な境界と可変水深を持つ空間）を、垂直座標が変換された正規化された「投影空間」へと投影することである。これにより、投影空間において楕円体粒子や簡略化された境界表現を使用することが可能になる。

実装は、以下の3つの主要な技術的コンポーネントに基づいている：

1. 座標変換： 物理的な垂直座標 $z$ を、スケーリング因子 $\alpha(x,y)$ と高さ補正項 $\beta(x,y)$ を用いて投影座標 $b_z$ にマッピングする変換が定義される。この変換には、変換された空間における支配方程式を正確に解くために、高次の微分（交差微分を含む）の計算が必要となる。
2. 高精度離散化（SPH(2)）： 座標変換によって導入される複雑な微分項を扱うため、著者らは**SPH(2)**と呼ばれる、2次精度の高いSPHモデルを採用している。このモデルは、座標変換に伴う精度と安定性を維持するために不可欠な、高精度な2次微分（交差微分を含む）を計算する。
3. 体積保存： $\alpha$ が変化する場合、座標変換によってヤコビアン行列式（ひいては粒子の代表体積）が変化するため、新しい体積保存補正技術が導入されている。これには、数値的不安定性や体積損失を防ぐために、ヤコビアンを平滑化し、圧力ポアソン方程式の安定化項を通じて粒子の体積を調整するプロセスが含まれる。

提案されている3つの手法
著者らはVCTの概念を以下の3つの特定の手法へと一般化している：

1. BF-SPH (Bottom Boundary-Fitted SPH)： 有限差分法の境界適合座標系に着想を得ている。これは、複雑な底面境界を投影空間における平坦な表面へと変換する（$\beta = -h(x,y)$、$\alpha=1$）。これにより、底面境界条件の適用が簡素化される。
2. BFE-SPH (Bottom Boundary-Fitted Ellipsoidal SPH)： BF-SPHと楕円体粒子モデルを組み合わせたものである。垂直方向のスケーリング因子 $\alpha > 1$ を用いて、物理空間において粒子を垂直方向に伸長させることで、底面適合の利点を維持しつつ、総粒子数を削減する。
3. $\sigma$-SPH： 海洋学（例：プリンストン海洋モデル）で一般的な$\sigma$座標系を利用する。ここでは、垂直スケーリング因子 $\alpha$ が水深に基づいて変化する（$\alpha = \bar{H} / (H - h(x,y))$）。これにより、段階的な解像度の遷移が可能となり、浅海域には高解像度を、遠方の深海域には粗い解像度を提供することで、総粒子数を大幅に削減できる。

主な結果
これらの手法は、静水圧、ダム崩壊問題、および斜面上の動的な流れの2次元および3次元シミュレーションを通じて検証された：

• 精度： VCTを用いたSPH(2)の使用は、標準的なISPHと比較して精度を大幅に向上させた。静水圧テストにおいて、SPH(2)の圧力における$L_2$誤差は、標準的なISPHよりも2桁小さかった。これらの手法は、複雑な底面形状やアスペクト比の変化がある場合でも、静水圧分布を正確に再現した。
• 効率性： $\sigma$-SPH法が最も高い効率を示した。静水圧問題において、$\sigma$-SPHは標準的なBF-SPH法と比較して、同等の精度を維持しながら、粒子数を約1/4に、計算時間を1/6に削減した。
• 安定性と体積保存： 提案された体積保存技術（動的問題における条件C）は、長時間のシミュレーション（最大100秒）においても、総体積誤差を$\pm 1\%$以内に維持することに成功した。一方、平滑化を行わない手法や、不適切な平滑化パラメータを用いた手法では、体積の減少や発散が見られた。
• 3次元能力： 3次元ダム崩壊シミュレーションにより、$\sigma$-SPH法が、低密度誤差と速度ダイバージェンスを維持しながら、波の遡上、反射、渦形成を含む複雑な3次元の流れ現象を堅牢に扱えることが確認された。

意義と主張
本論文は、この研究が粒子法に$\sigma$座標系を適用した最初の試みであることを主張している。著者らは、これらの手法の成功が、座標変換によって要求される2次微分を正確に計算するための**SPH(2)**の使用に大きく依存していることを強調している。

主な意義は、以下の点にある：

• 粒子法に対してVCTを一般化し、ラグランジュ的な粒子法の柔軟性と、座標変換を用いた格子法の効率性の間の溝を埋めたこと。
• 自由表面や境界条件の精度を損なうことなく、可変の垂直解像度を許容することで、海洋学的および沿岸部のシミュレーション（例：津波）における計算コストを劇的に削減できること。
• 複雑な底面地形を計算領域内の平坦な面にマッピングすることで、境界処理を簡素化できること。

著者らは、本研究は主にこれらの手法の3次元における実現可能性を確認するものであるが、$\sigma$-SPH法は非バイアスな手法に対して5倍の計算効率向上を実現しており、沿岸および海洋工学の実用的なアプリケーションにとって有望なツールであると結論付けている。また、垂直解像度を制御できる能力は、将来的に近壁乱流の解像において利点をもたらす可能性があるが、これにはさらなる調査が必要であるとも述べている。