Liggghts_VBP_VTKReader

A simple tool that parses VTK files generated by LIGGGHTS for Vibrating Packed Bed simulations and calculates particle counts in slices of the packed bed.

How to build

You need to install vtk to build this project. It has been tested using vtk6.3. This also depends on matplotlibcpp to plot graphs. A fork of the original matplotlibcpp has been used. Therefore you also need to install matplotlib and the python2-dev libs.

sudo apt-get install libvtk6.3 libvtk6-dev
sudo apt-get install python-matplotlib python2.7-dev
pip install matplotlib
pip install numpy

From root directory of the project:

git clone --recurse-submodules https://github.com/KarthikRIyer/Liggghts_VBP_VTKReader
cd Liggghts_VBP_VTKReader
mkdir build
cd build
cmake -DCMAKE_BUILD_TYPE=Release ..
make

How to run

From build directory:

This will write the output to a file in the same directory as the vtk file with the extension .postprocessed

./Liggghts_VBP_VTKReader <path to vtk file>

This will write the output to files in a new directory _post_processed in the input directory.

./Liggghts_VBP_VTKReader <path to directory with the vtk files>

You can use the following flags to limit your output:

-notall This will process only the first and last files (according to the index)
-onlylast This will process only the last file (according to the index)
-initframe <frame number> This will calculate the bed height for the frame number specified and will consider it for all calculations

To plot the velocity field of the vertical middle slice of the bed, install matplotlib and run in the root directory of the project:

python vel.py <path to the .vel file or path to the 'post' directory where the LIGGGHTS sim data resides>