[論文レビュー] A Robust Compressible APIC/FLIP Particle Grid Method with Conservative Resampling and Adaptive APIC/PIC Blending

Jiansheng Yao, Yingkui Zhao|arXiv (Cornell University)|Mar 4, 2026

Block Copolymer Self-Assembly被引用数 0

ひとこと要約

この論文は、 conserva t ive resampling と滑らかな APIC/PIC ブレンドを追加して、長時間RTIのスパイク頭ボイドを修正しつつ渦動力学を保持することで FLIP–APIC にテンソル人工粘性を拡張する。

ABSTRACT

Modeling inviscid compressible flows with shocks and vortex dominated dynamics remains challenging for particle grid methods due to moving discontinuities, cell crossing noise, and quadrature degradation under strong deformation. Building on a FLIP/APIC framework with vorticity aware tensor artificial viscosity, we identify a long time RTI failure mode: particle depletion at spike heads degrades quadrature and particle grid coupling, producing nonphysical, void-like dents. Standard mitigations (CPDI lite and subcell-jittered seeding) reduce but do not eliminate this artifact. We therefore add two sampling-aware controls: (i) conservative split resampling that replenishes depleted cells while exactly conserving mass, momentum, and internal energy; and (ii) a soft-switch that attenuates only the APIC affine term when local support is insufficient. Tests on the Sod shock tube and single/multi mode RTI show that the method removes spike head voids in long-time RTI while preserving vortex roll up, and matches reference Euler growth metrics

研究の動機と目的

• 長時間RTIおよび衝撃波支配流における粒子-格子ソルバーのサンプリングに起因する故障モードを対処する。
• 局所的な四分割品質を回復し、アフィン APIC 成分を安定化させるサンプリング対応制御を開発する。
• 適切にサンプルされている領域で渦豊かな輸送を保持しつつ、テンソル人工粘性の頑健性を保つ。
• ソド・ショックチューブおよび単一/複数モードRTIベンチマークで有効性を示す。

提案手法

• 不完全な粒子分布下で量子分布性のバイアスを低減するための CPDI-lite 転送を採用する。
• 質量・運動量・内部エネルギーの厳密な保存とともに depleted セルを補充する保守的分割再サンプリングを導入する。
• 局所的な粒子サポートが不十分な場合に APIC のアフィン項を PIC 系列へと徐々に減衰させるセル毎のソフトスイッチを実装する。
• 局所サポートと流れの特徴に基づく APIC アフィン成分の二者指標ゲーティングを用いる。
• FLIP–APIC 転送と渦動性認識テンソル人工粘性を組み合わせてショックを安定化させつつ渦動構造を保存する。
• 過少サンプル領域でのみ APIC アフィン項を漏出させ、不正なエネルギー注入を防ぐ。

Figure 1: Richtmyer–Meshkov instability in the FLIP–APIC + vorticity-aware tensor-AV baseline. Particles are colored by material (red/blue), illustrating post-shock interfsace roll-up at different time.

実験結果

リサーチクエスチョン

• RQ1 長時間RTIにおける spike heads 近傍のサンプリング不足が、圧縮性粒子-格子シミュレーションに非物理的なへこみやボイドを引き起こすか？
• RQ2 保守的リサンプリングと適応的 APIC/PIC ブレンディングはスパイクヘッドボイドを緩和しつつ渦巻きの巻き上げ忠実性を損なわないか？
• RQ3 CPDI-lite 転送とサブセルジッタリングは提案制御とどう相互作用して長時間の安定性を改善するか？
• RQ4 拡張手法はベンチマークの Euler 成長指標とショーカプチャー性能を保つか？

主な発見

• 保守的な分割再サンプリングは depleted セルを補充し、質量・運動量・内部エネルギーを厳密に保存して局所の四分割品質を回復する。
• APIC のアフィン項を PIC 類似の転送へとソフトスイッチすることにより、枯渇領域でのアフィインバリデーションを防ぎつつ、よくサンプルされた領域で APIC の利点を維持する。
• CPDI-lite 転送とサブセルジッタリングは格子の imprinting と四分割バイアスを低減し、スパイクヘッドアーティファクトを緩和する。
• 完全手法は長時間RTI におけるスパイクヘッドのへこみを除去しつつ渦巻きの巻き上げを保持し、参考となる Euler 成長指標に整合する。
• ベースライン RM/RTI テストは、未改良の APIC/FLIP フレームワークと比較して頑健なショック捕捉と改良されたせん断層忠実度を示す。

Figure 2: Single-mode RTI in the FLIP–APIC + tensor-AV baseline: a nonphysical spike-head dent/void-like depression develops in long-time evolution (highlighted region).