文献：GeoDualSPHysics: a high-performance SPH solver for large deformation modelling of geomaterials with two-way coupling to multi-body systems

GeoDualSPHysics: a high-performance SPH solver for large deformation modelling of geomaterials with two-way coupling to multi-body systems

AI要約

1. 背景（Background）

•  従来の課題: 地盤工学や防災工学における地滑り、土石流、斜面崩壊などの現象は、土砂が液体のように振る舞う「極大変形（大流動）」を伴います。従来の有限要素法（FEM）などの格子基盤の手法では、メッシュが過度に歪むため、このような大変形を安定して計算することが困難でした。
• 粒子法（SPH）の台頭と限界: メッシュを用いない「粒子法（SPH）」は大変形問題に適していますが、計算コストが非常に高く、数百万〜数千万粒子に及ぶ広大な実地形の大規模シミュレーションを現実的な時間で処理することが困難でした。
•  開発の動機: 既存の非常に強力なオープンソースGPU加速型流体ソルバー「DualSPHysics」の並列計算フレームワークを基盤とし、地盤固有の力学挙動（弾塑性）や構造物との複雑な相互作用を、数億粒子スケールで超高速に計算できる実用的な地盤用ソルバーの確立が必要とされていました。

2. 手法（Methodology）

• •連続体力学のラグランジュ定式化:地盤材料を連続体として捉え、質量保存則および運動量保存則をSPHの粒子離散化形式で解きます。カーネル関数には、粒子の異常な凝集（クランピング不安定性）を防ぎ安定性に優れる「5次Wendland（C2）カーネル」が採用されています。
• 地盤の構成モデル（Constitutive Models）:材料の応力状態を評価するため、弾塑性モデル（例：Drucker-Prager Yield Criterion など）を組み込んでいます。毎タイムステップで、まず弾性試行応力を計算し、降伏条件を超えた場合には「リターンマッピングアルゴリズム（Return Mapping Algorithm）」を用いて、応力を降伏曲面上に引き戻すプラスチックひずみ計算を行います。
•  ノイズ除去応力処理技術（Noise-free Stress Treatment）:従来のSPHで課題となっていた「応力場の数値的なノイズや振動」を抑制・安定化させ、高精度に応力分布を解像する特別な平滑化・拡散技術を導入しています。
• 改良型境界条件（Modified Dynamic Boundary Conditions / mDBC）の拡張:固体壁や地盤境界での粒子配置の乱れを補正するため、境界粒子の外側に仮想的なゴーストノードを投影し、周囲の材料粒子から変数を1次あるいは0次の補間（Shepard補間）で定義する手法を地盤用に拡張・適用しています。
• マルチボディシステム（MBS）との双方向結合:物理演算エンジン Project Chrono と双方向（Two-way）でカップリングしています。SPH側で構造物（強体）粒子にかかる地盤・流体からの力を集計してChronoに渡し、Chrono側で構造物の並進・回転運動や衝突、拘束条件を計算して再びSPH側へ戻すというプロセスを、時間積分ステップごとに高同期でループさせます。
•  2相流（水–土骨格）フレームワークの統合:間隙水圧（pw）と土骨格の変位（u）を連成して解く水・土砂の2相流モデルが組み込まれており、雨水の浸透による有効応力の低下や地盤内の圧密現象（Terzaghiの圧密など）を機械的・水理学的境界条件のもとで同時に扱えるよう拡張されています。

3. 結果（Results）

• 圧倒的な計算パフォーマンス（GPU加速）:CUDA並列化スキームの最適化により、単一のGPU（Single GPU）であっても最大で1億粒子オーダーという極めて高解像度かつ大規模な地盤変形計算を実行できることが実証されました。
• ベンチマークによる検証:1次元圧密試験（Terzaghi理論解との比較）や、球体サンプルの水理負荷試験（Cryer効果の再現）において、計算された間隙水圧の推移や地盤の沈下量が理論値・参照解と高い精度で一致することが確認されました。
• 実規模の崩壊シミュレーション:5メートルおよび8メートル規模のデリケートな粘土質斜面（Sensitive-clay slopes）の崩壊実験では、雨水の浸透に伴って斜面が段階的に後退崩壊（Retrogressive failure）を起こし、最終的に大流動へと至る複雑なシーケンスを破綻なく再現することに成功しました。

4. 考察（Discussion）

• 本手法の優位性:GeoDualSPHysicsは、これまで「大規模すぎて計算が終わらない」、あるいは「メッシュの歪みで途中で止まる」とされていた、「豪雨・浸透 → 斜面崩壊 → 構造物への衝突・破壊」という一連の連鎖的災害（Cascading disasters）を一つのシステム内で一気通貫かつ超高速にシミュレートできる点が極めて強力であると評価されています。MBS（Project Chrono）との結合は、防護壁や杭などの防災構造物の耐力評価において非常に実用的です。
• 今後の課題と限界:一方で、2相流シミュレーションにおける「初期応力状態（Prestress）の設定」や、非常に透水係数が低いケースにおける「間隙水圧の収束性と計算安定性の担保」にはまだ改善の余地が残されています。また、人工的な数値粘性（Artificial Viscosity）や応力拡散（Stress Diffusion）などのパラメータ設定に一部依存する箇所があり、これらを物理的な土質試験データとどのように厳密にキャリブレーション（校正）していくかが、今後の予測精度向上の鍵として議論されています。

土水2層のSPHです。DualSPHysicsベースです。
私が作成したものに比べ、圧密、不飽和透水係数の実装で、こちらの方が進んでいます。良いですね。弾塑性とビンガム流体の切り替えは、実務で使うには調整が必要そう。私も個々は割り切りましたが、ひと段落付けば他の研究成果も取り入れるべきでしょうね。
バグを見つけて修正する速度がAIの進化で劇的に速くなりました。また、機能の追加も同様です。このようにAIを利用してコードを作っていくと、結局同じようなコードができてしまいそうです。不飽和や圧密も容易に実装できると思います。ま、過去の手法の実装や論文に出ている内容は容易に再現できるようになったので、その先を考えられるために時間をとれるようになったと考えるべきなのでしょう。