This repository contains the essential files solving numerically a 2D flow model (downscaled to 1D for simplicity). The model solved here was proposed originally by Iverson and Denlinger in 2001 (DOI: 10.1029/2000JB900330). There are 10 files:
(1) The main code "FVM_ID_Final_base.m". Here the solver is implemented. (2) "fCoef_v1.m" an auxiliar function for the coefficients of the system of equations. (3) "fgx1D.m" an auxiliary function for rotating the axis of acting forces parallel to the direction of movement. (4) "fgz1D.m" an auxiliary function for rotating the axis of acting forces perpendicular to the direction of movement. (5) "fh0.m" auxiliar funtion for the definition of the initian conditions. (6) "fjacobi_v3_base.m" auxiliar function with the Jacobi method solving the equivalent system of linear equations. (7) "flambda.m" auxiliar function declaring the pore pressure ratio. (8) "fThetta1D" auxiliar function declaring the inclination angle along the domain.
There are two non-essential additional files for model validation purposes:
(9) A video showcasting the flow depth and instant velocity development "Model_h_u.mp4". (10) A video comparing the experimental granular flow with the calibrated model "DataS3P1_S3P15_Cal.mp4".