[論文レビュー] A rain induced landslide 3D model based on molecular dynamics with fractal and fractional water diffusion

Gianluca Martelloni, Franco Bagnoli+1 more|arXiv (Cornell University)|Nov 18, 2016

Landslides and related hazards参考文献 51被引用数 3

ひとこと要約

本稿では、不連続拡散を用いた3次元分子動力学モデルを提案し、多孔質土壌における水の浸透を記述するため、フラクタルおよび分数階微分方程式を統合する。粒子レベルでのモール-クーロン破壊基準を用いた粒子ベースのMDシミュレーションと、フラクタル拡散偏微分方程式の解析的・数値的解を結合することで、本モデルは滑落の発火、伝播、分離、およびべき乗則スケーリングを的確に再現し、実世界の滑落動態への応用可能性を検証した。

ABSTRACT

We present a three-dimensional model, based on cohesive spherical particles, of rain-induced landslides. The rainwater infiltration into the soil follow the either the fractional or the fractal diffusion equations. We solve analytically the fractal diffusion partial differential equation (PDE) with particular boundary conditions to simulate a rainfall event. Then, for the PDE, we developed a numerical integration scheme that we integrate with MD (Molecular Dynamics) algorithm for the triggering and propagation of the simulated landslide. Therefore we test the numerical integration scheme of fractal diffusion equation with the analytical solution. We adopt the fractal diffusion equation in term of gravimetric water content that we use as input of triggering scheme based on Mohr-Coulomb limit-equilibrium criterion, adapted to particle level. Moreover, taking into account an interacting force Lennard-Jones inspired, we use a standard MD algorithm to update particle positions and velocities. Then we present results for homogeneous and heterogeneous systems (i.e. composed by particles with same or different radius respectively). Interestingly, in the heterogeneous case, we observe segregation effects due to the different volume of the particles. Finally we show the parameter sensibility analysis both for triggering and propagation phase. Our simulations confirm the results of our previous two-dimensional model and therefore the feasible applicability to real cases.

研究の動機と目的

• 本研究の目的は、多孔質土壌における異常拡散をモデル化することで、滑落シミュレーションの物理的リアリズムを向上させることにある。
• 本研究は、従来の2次元モデルを3次元に拡張し、滑落伝播のダイナミクスをよりよく捉えることを目的としている。
• 本研究の目的には、不均一な土壌における非ガウス的かつ記憶依存的水の拡散をモデル化するためのフラクタルおよび分数階微積分の使用を検証することを含む。
• 本モデルは、粒子レベルの相互作用を通じて、亀裂、アーチング、圧縮領域といった重要な滑落特徴を再現することを目的としている。
• 本研究は、モデルのロバストネスと粒状物質の挙動との整合性を保証するため、パラメータ感度を調査する。

提案手法

• モデルは、傾斜面に配置された3次元格子上に初期静止状態の粘着性球状粒子を用いる。
• 水の浸透は、時間および空間微分の次数がそれぞれαおよびβであるフラクタル拡散偏微分方程式（PDE）を用いてシミュレートされ、拡散係数Dは定数または深さ依存とする。
• 特定の境界条件下でフラクタル拡散PDEの解析解が導出され、土壌の不均一性をモデル化するための確率的駆動項を伴うPDEに対して、数値積分スキーム（陰的アダムス・バッシュフォースト＝ムルトン法）が開発された。
• モール-クーロン破壊基準が粒子レベルに適応され、拡散モデルからの水分含有量に基づいて滑落の発火を決定する。
• 粒子の運動は、Lennard-Jones型の相互作用ポテンシャルを用いた標準的な分子動力学アルゴリズムにより更新される。
• 自然な土壌の不均一性を反映させるために、パラメータに確率的変動を組み込み、粒子半径を変化させた均一および不均一な粒子系をテストした。

実験結果

リサーチクエスチョン

• RQ1 フラクタルおよび分数階微分方程式を用いることで、古典的拡散と比較して、多孔質土壌における水の浸透シミュレーションはどの程度向上するか？
• RQ2 フラクタル拡散PDEにおける時間および空間微分の次数（α, β）を変化させた場合、滑落発火時間および伝播ダイナミクスにどのような影響を与えるか？
• RQ3 粒子サイズの違いが、シミュレートされた滑落における分離および流れのパターンに与える影響は何か？
• RQ4 モデルの統計的出力（例：速度、発火時間）がべき乗分布を示す程度はどの程度か。これは、複雑かつ極端な出来事の行動を示唆する。
• RQ5 摩擦係数、粘性係数、相互作用力パラメータの変動に対して、最終的な滑落速度はどの程度感度を示すか？

主な発見

• 本モデルは、3次元における滑落発火および伝播を的確に再現し、従来の2次元シミュレーションおよび実際の滑落動態と整合性がある。
• フラクタル拡散PDEにおける時間微分次数（α）が低下するほど、発火時間が延長され、これは準拡散的挙動が飽和および破壊を遅らせる。
• 発火時間（R² = 0.9998）および空間微分次数への感受性（R² = 0.9996）に対するべき乗則フィットは、モデルが複雑系と統計的に整合していることを確認した。
• 不均一系では、粒子サイズの違いが分離効果を引き起こし、より大きな粒子が上昇する。
• 最終的な滑落平均速度は、低摩擦係数（R² = 0.9929）および粘性係数（R² = 0.983）に対してべき乗則的依存を示し、強い非線形感受性を示した。
• 本モデルの伝播軸に沿った速度プロファイルは、1次元有効モデルの予測と一致し、物理的整合性を裏付けた。