Discrete Element Model Simulations for a Granular Bed Undergoing Shear Deformation

Discrete Element Model (DEM) Archive for, "Shear Variation at the Ice-Till Interface Changes the Spatial Distribution of Till Porosity and Meltwater Drainage" by, Indraneel Kasmalkar, Anders Damsgaard, Liran Goren, Jenny Suckale This dataset contains the output files of the DEM simulations for understanding how the porosity of a granular bed evolves with an imposed laterally varying shear. The DEM simulations are performed using Sphere (https://src.adamsgaard.dk/sphere/), developed by Anders Damsgaard. The output files are in the 'output' folder. The porosity values are curated and stored in 'porosities' folder as 4D numpy arrays. The results are generated as .pdf files in the 'Images' folder. The scripts are stored in the 'python' folder. The main file in the 'python' folder is MyVisualization.py. To run the script you will need Python 3.6 or newer, with the following packages: numpy, scipy, matplotlib, colorbrewer, math, pickle In MyVisualization.py, go to line if __name__ == '__main__' Before the line, you can change the variable plot_type to get the plot you need. The options for plot_type are listed in a comment right before the variable. You can limit the plotting to either simple shear ('unif') or laterally varying shear ('hl') within each if condition thereafter