技術要約：非線形波と3D流体のハミルトニアン二方向結合

1. 問題提起

コンピュータグラフィックスにおける大規模な自由表面流のシミュレーションでは、ハイブリッドなアプローチが必要となることが多い。すなわち、高精度で計算コストの高い3Dナビエ・ストークス（NS）ソルバーを局所的な領域（船体や落下する物体の周囲など）に限定し、周囲の開けた海洋には、より計算コストの低い2D表面モデルを用いる手法である。

この2D–3D分解における主要なボトルネックは、結合界面における波のダイナミクスの根本的な不一致にある。既存の効率的な2D波モデルには以下の2つのカテゴリーが存在するが、いずれも強力な非線形3Dソルバーとシームレスに結合させることはできない。

• 線形モデル（例：エアリー理論）: 正しい分散（波長依存の速度）を捉えるが、非線形相互作用（波の急峻化など）を切り捨てる。これを非線形3Dソルバーと結合させた場合、2D側が3D側から生成される振幅依存の分散や波峰の急峻化を表現できないため、界面で目に見える偽の反射やアーティファクトが発生する。
• 非分散モデル（例：浅水波方程式 - SWE）: 非線形な移流は捉えるが、分散関係を崩壊させ、すべての波を波長に関わらず同じ速度で伝播させる。これは、実際の海洋が持つ明確な分散挙動を再現できない。

これらを解決するためのBEM（境界要素法）などの従来の手法は、非線形な精度を取り戻せるものの、膨大な計算コスト（$O(N^3)$）と複雑なメッシュ追跡を必要とするため、大規模なリアルタイムまたはインタラクティブなアプリケーションには実用的ではない。

2. 手法

著者らは、局所的な3D NSソルバーを、正準なザカロフ定式化（Zakharov formulation）に基づく非線形かつ分散的な2Dオイラー型波モデルと結合させる新しいフレームワークを提案している。

2.1 2D波ソルバー：ハミルトニアン・スペクトル・アプローチ

本手法の核となるのは、正準なペアである表面標高（$\eta$）と表面速度ポテンシャル（$\psi$）を、正確な水波ハミルトニアンの下で進化させる2D波ソルバーである。

• ザカロフ定式化: システムは、表面変数のみで表現された運動エネルギーとポテンシャルエネルギーの和であるハミルトニアンに従うハミルトンの方程式によって支配される。
• ディリクレ・ネマン演算子（DNO）: システムを閉じるために、本手法はCraig–Sulem展開によるDNOを採用している。これにより、毎ステップで完全な3Dラプラス方程式を解くことなく、移動する表面上での速度ポテンシャルの法線微分を評価することが可能になる。
• 高次スペクトル（HOS）展開: DNOは、$\eta$の同次べき級数のテイラー展開によって展開される。著者らは項を第2次または第3次（HOS-2またはHOS-3）まで保持し、非線形補正を高速フーリエ変換（FFT）と規則的なグリッド上の点別積へと簡略化している。これにより、$O(N^2 \log N)$の計算複雑性を達成している。
• 制御可能な非線形性（$\varepsilon$）: スカラー摂動パラメータ $\varepsilon \in [0, 1]$ が非線形項を調整する。
  • $\varepsilon = 0$ の場合：線形エアリー理論に帰着する。
  • $\varepsilon = 1$ の場合：完全に非線形な（切り捨てられた）ザカロフ系を再現する。
    これにより、ユーザーは安定性と精度のトレードオフを制御でき、急峻な砕波に近い波が発生するシナリオでは $\varepsilon$ を下げて数値的不安定性を回避できる。
• 時間積分: ソルバーは厳密線形アダムス・バッシュフォース（Adams–Bashforth）スキームを使用する。線形部分（分散）は剛性を避けるためにフーリエ空間で厳密に積分され、非線形強制力は2次のアダムス・バッシュフォース外挿を用いて前進される。

2.2 二方向結合メカニズム

結合は**正準的に一貫（canonically consistent）**するように設計されており、これは情報が物理的な対称性を壊すことなく、2Dの正準状態 $(\eta, \psi)$ と3Dバルク状態の間をスムーズに流れることを意味する。

• 2D $\to$ 3D（リフティング）: 2Dの状態 $(\eta, \psi)$ を3Dエンベディングボックスの境界へと持ち上げる。Craig–Sulem展開を用いて、平坦な基準面 ($y=0$) におけるバルク速度ポテンシャル $\Phi$ を再構成する。その後、有限水深の調和カーネルを用いてこれを水深層へと拡張し、3DソルバーのFAB（Fluxed Animated Boundary）層に対して水深分解された速度境界条件を提供する。
• 3D $\to$ 2D（フォールディング）: 3Dソルバーは表面標高 $\eta_{3D}$ を露出させる。不連続性を引き起こす直接的な上書きではなく、本手法ではザカロフ系の運動方程式に対して暗黙的な緩和ソースを適用する。このソースは、遷移帯内で2D標高 $\eta$ を $\eta_{3D}$ へと指数関数的に減衰させる。極めて重要な点は、共役変数 $\psi$ が、緩和された $\eta$ とハミルトニアン進化から導出された表面垂直速度を用いて再構成されることである。これにより、2D状態が緩和後も有効な正準ペア $(\eta, \psi)$ として維持され、偽の反射を防ぐことができる。

3. 主な貢献

1. 非線形的かつ分散的な2D波モデル: 非線形な波の急峻化と正しい分散を同時に捉える、グラフィックスのための新しいオイラー型波モデル。BEMの $O(N^3)$ の複雑さやメッシュ操作を回避しつつ、正準なザカロフ・ペアを $O(N^2 \log N)$ のコストで進化させる。
2. 正準的に一貫した二方向結合: 波モデルのハミルトニアン構造を活用した結合戦略。界面交換を正準変数間の写像として扱い、暗黙的緩和を用いることで、界面のアーティファクトを抑制し、エネルギー整合的な情報伝達を実現する。
3. 連続的に調整可能な非線形性: パラメータ $\varepsilon$ の導入により、連続的なファミリーのソルバーを提供し、ユーザーがシーンごとに非線形な忠実度と数値的安定性のバランスを取ることを可能にする。

4. 結果と検証

本手法は、高解像度の3Dオイラー参照モデルおよび従来の結合技術を用いて、様々なベンチマークにおける検証が行われた。

• 精度: 波の伝搬テストにおいて、提案モデルはSWE、BEM、およびエアリーのベースラインと比較して、平均波高誤差を1.7〜5倍減少させた。ストークス波のテストでは、完全非線形モデル ($\varepsilon=1$) は解析解と位相が一致し続けたのに対し、線形モデルは著しくドリフトした。
• パフォーマンス: 2DソルバーはBEMよりも1000倍以上高速に動作する。結合シナリオ（例：単一のボート）では、フルシステムは同サイズのフルドメインGPU NB-FLIPシミュレーションよりも4倍以上高速に動作する。
• 結合の忠実度:
  • 分散のマッチング: 提案手法は、最小限の継ぎ目アーティファクトを持つ対称的な波場を生成した。一方、従来のメソッド（SWE、エアリーベース）では、非対称性や偽の反射が見られた。
  • 複雑なシナリオ: 水上飛行機の着水、潜水艦の浮上、戦艦の荒天航行、および二隻のボートの追跡をシミュレートすることに成功した。すべてのケースにおいて、2D–3D界面には目立ったアーティファクトがほとんどなく、航跡が境界を越えてスムーズに伝播した。
• 安定性: 完全非線形モデル ($\varepsilon=1$) は非常に高精度であるが、極端な強制力（例：戦艦のシーンにおける砕波に近い波）に対しては、安定性を維持するために $\varepsilon$ を下げる（例：0.2）必要があることが著者らによって指摘されている。

5. 意義と主張

本論文は、グラフィックスにおける流体シミュレーションにおける長年の構造的な不一致に対処することを目的としている。非線形かつ分散的な2Dモデルを導入することで、著者らは、複雑な流体・固体相互作用を伴う大規模な海洋シミュレーションに不可欠な、実用的で正確であり、かつ界面アーティファクトがほとんどない2D–3D分解パイプラインを可能にした。

その意義は、結合のハミルトニアン構造にある。2Dソルバーが正準ペア $(\eta, \psi)$ を進化させるため、結合の方向（3Dへのリフティングと3Dからのフォールディング）が、アドホックなデータ転送ではなく、自然な数学的操作（ディリクレ・ネマン写像）となる。これにより、従来の線形または非分散モデルでは視覚的な劣化なしには扱えなかった、強力な非線形相互作用（船の航跡やしぶきなど）を扱うことが可能になる。

著者らは本手法を一つのトレードオフとして位置づけている。すなわち、線形ソルバーと比較して計算コストのわずかな増加と無条件の安定性のわずかな低下を受け入れる代わりに、大規模な海洋シミュレーションに求められる非線形な精度と結合の忠実度を獲得している。また、単一値の高さ場に起因する「オーバーターン（転倒）を伴う砕波」のシミュレーションができないことや、極端なシナリオにおける $\varepsilon$ の手動調整の必要性といった限界も認めているが、本手法が幅広いプロダクションレベルのシナリオにおいて実用的であることを主張している。