[论文解读] 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

一句话总结

本文提出了一种3D分子动力学模型，用于模拟降雨诱发的滑坡，通过整合分形与分数阶扩散方程，模拟多孔土壤中的水分入渗。通过将分形扩散PDE的解析解与数值解，与基于粒子的MD模拟相结合，并在粒子层面采用Mohr-Coulomb破坏准则，该模型成功再现了滑坡的触发、传播、分异及幂律标度——验证了其在真实世界滑坡动力学中的适用性。

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.

研究动机与目标

• 本研究旨在通过在多孔土壤中建模异常扩散的水分入渗，提升滑坡模拟的物理真实性。
• 旨在将先前的2D模型扩展至3D，以更好地捕捉滑坡传播的动力学行为。
• 目标包括验证分形与分数阶微积分在模拟非高斯、具有记忆依赖性的水分扩散于非均质土壤中的适用性。
• 该模型旨在通过粒子层面的相互作用再现滑坡的关键特征，如断裂、拱形区与压缩区。
• 研究参数敏感性，以确保模型的鲁棒性，并与颗粒物质行为保持一致。

提出的方法

• 模型采用粘结的球形粒子在倾斜平面上排列成3D晶格，初始时静止。
• 水分入渗通过分形扩散偏微分方程（PDE）模拟，时间与空间导数阶数分别为α与β，扩散系数D为常数或随深度变化。
• 在特定边界条件下，推导出分形扩散PDE的解析解，并开发了用于含随机强迫项PDE的数值积分方案（隐式Adams-Bashforth-Moulton格式），以模拟土壤非均质性。
• 在粒子层面适配Mohr-Coulomb破坏准则，基于扩散模型提供的含水量判断滑坡触发。
• 通过标准分子动力学算法，利用Lennard-Jones型相互作用势更新粒子运动。
• 模型引入参数的随机变化以反映自然土壤的非均质性，并测试了均质与非均质粒子系统（粒子半径不同）的情况。

实验结果

研究问题

• RQ1 与经典扩散相比，使用分形与分数阶扩散方程在模拟多孔土壤中水分入渗方面有何改进？
• RQ2 在分形扩散PDE中，时间与空间导数阶数（α, β）的变化对滑坡触发时间与传播动力学有何影响？
• RQ3 粒子尺寸差异如何影响模拟滑坡中的分异与流动模式？
• RQ4 模型的统计输出（如速度、触发时间）在多大程度上符合幂律分布，表明其具有复杂极端事件行为？
• RQ5 最终滑坡速度对摩擦系数、粘度及相互作用力参数的变化有多敏感？

主要发现

• 该模型成功在3D中再现了滑坡的触发与传播，结果与先前2D模拟及真实滑坡动力学一致。
• 随着分形扩散PDE中时间导数阶数（α）减小，触发时间延长，表明亚扩散行为会延迟饱和与破坏。
• 触发时间（R² = 0.9998）与空间导数敏感性（R² = 0.9996）的幂律拟合结果，证实了模型在统计上与复杂系统的一致性。
• 在非均质系统中，粒子尺寸差异导致分异效应，较大粒子向上传播。
• 最终滑坡平均速度对较低摩擦系数（R² = 0.9929）与粘度（R² = 0.983）呈幂律依赖，表明其具有强烈的非线性敏感性。
• 沿传播轴的速度分布与一维有效模型的预测一致，验证了模型的物理一致性。