double-layer two-phase formulation その2

振り返ると、double-layer での実装に取り掛かってから１か月以上かかっています。

以下を参考に始めたのですが、試行錯誤の結果、最終的にはPersianSPHの文献と同じような形に落ち着きました。
A coupled fluid-solid SPH approach to modelling flow through deformable porous media - ScienceDirect

SPH approach for simulating hydro-mechanical processes with large deformations and variable permeabilities | Acta Geotechnica (springer.com)

２次元、３次元共に土、水の相互作用までは動くのですが、地下水が地上に出た場合の密度変化をうまく表現できません。相互作用の検証も終わっていませんので、まだ時間がかかりそうです。

SPHでの２層2相表現は定性的にとどまると伺ったことがあります。それほど新しいものでもないので、国内でも実装されている方はいらっしゃいます。その先生がおっしゃるのですから、そうなのかもしれません。

まだ１か月。もう少し追いかけてみましょう。