Configurations and scripts to reproduce the numerical simulations of "A two-fluid model for immersed granular avalanches with dilatancy effects" article

This directory contains the main data to reproduce the results presented in the article "A two-fluid model for immersed granular avalanches with dilatancy effects" by Eduard Puig Montellà, Julien Chauchat, Bruno Chareyre, Cyrille Bonamy and Tian-Jian Hsu.

The numerical results and experimental data extracted from Pailha et ad (2008) can be found inside "NumericalData" and "ExperimentalData" folders respectively. The script "PressureVelocityPlot.py" displays the evolution of the surface particle velocity and the excess of pore pressure with time for cases ranging from loose to dense granular avalanches. 

The input files needed to reproduce a dense granular avalanche (phi=0.592) in 1D and 2D are found in the following folders: "1D_DenseCase" and "2D_DenseCase". To accelerate the simulation, the files are given after 200 seconds of sedimentation in order to reach an equilibrium state. Please read the corresponding README.txt files to launch a 1D and/or a 2D simulation. Additionally, python scripts in each configuration are provided to evaluate the evolution of the main parameters during the avalanche.