LF_DEM Simulation Analysis Scripts

---

Particle Interactions Visualization of contact networks. The lubrication (green) interactions and frictional (red) interactions in a dense suspension.

Particle Velocities Angular velocities shown by particle color (color bar for magnitute); translational velocities by black arrows in a dense suspension.

This repository contains scripts and tools for analyzing data generated by LF_DEM (Lubrication Flow - Discrete Element Method) simulations. It provides utilities to read, process, and visualize particle and interaction data to help researchers investigate the behavior of particulate systems. Feel free to reach out for any questions: pandare.rahul@gmail.com.

Related Publications: This repository contains scripts for work done in the recent publication by our group - Contact network structures and rigidty development in bidisperse suspenssions.

Data availability: This repository does not include simulation data due to its large size. Data can be provided upon request, and we can work out a suitable way to share it. Altough a sample of 2D data is available here.

Users can also generate their own data by running the LF-DEM code. Detailed installation instructions are available here: LF-DEM installation.

Directory Structure

LF-DEM_analysis/
├── assets/ Sample images and figures
├── notebooks/ Jupyter notebooks for data analysis
├── src/ Source code and functions for LF-DEM data processing
├── tests/ Test scripts for validating functions

Features

• Parse and load rigid cluster data, particle sizes, interaction lists, and particle parameters.
• Perform statistical and structural analysis on simulation outputs.
• Unit tests to ensure correctness of data processing functions.

Installation

Clone the repository:

git clone https://github.com/yourusername/LF_DEM_analysis.git
cd LF_DEM_analysis

Install the required libraries:

pip install -r requirements.txt

To use the functions as a system-wide Python package, run the following from the root (after the previous step):

pip install -e .

This allows you to import functions for reading simulation data and calculating rigid clusters in any Python script or notebook.

Perform unit tests for test functions:

cd test/
bash test_functions.sh

If the tests pass, the functions are ready to use. Badge indicating test status: