A rain induced landslide 3D model based on molecular dynamics with fractal and fractional water diffusion

TL;DR

This paper introduces a 3D molecular dynamics model for rain-induced landslides, incorporating fractal and fractional water diffusion equations to simulate infiltration and triggering, validated through analytical and numerical methods.

Contribution

The study develops a novel 3D MD-based landslide model integrating fractal water diffusion and particle interactions, extending previous 2D models to more realistic scenarios.

Findings

Model successfully simulates landslide triggering and propagation.

Heterogeneous systems show segregation effects.

Parameter sensitivity analysis highlights key influences.

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.