PyLaBolt
PyLaBolt is a single phase, 2D, parallel lattice Boltzmann solver for fluid flow. It uses Numba accelerated Python code to run lattice Boltzmann simulations on 2D lattices. Simulations can be run on CPU in parallel via Numba's own OpenMP parallelization and the mpi4py library. For running on NVIDIA GPUs, PyLaBolt uses Numba's CUDA bindings.
Installation and Usage
PyLaBolt can be installed via pip
via the command:
$ pip install pylabolt
More details on the dependencies and their configuration can be found in the PyLaBolt documentation here.
Tutorial cases are provided in the tutorials
folder. For example, consider the lid driven cavity
problem in tutorials/cavity/Re_100/
. The configuration file that defines the simulation is
called simulation.py
. After installation just run the following command from the tutorials/cavity/Re_100/
folder:
$ pylabolt --solver fluidLB
The output data is written in the output
folder. By default the data is written into binary files with
.dat
extension. To visualize the data in Paraview /
Mayavi, the VTK library is used.
For example, to convert the last time-step data to a .vtk
file, run the following command from
the working directory:
$ pylabolt --toVTK last
The output_<time-step>.vtk
files are stored in output/VTK
directory which can be opened in
Paraview/Mayavi. Sample results for lid driven cavity, plane Poiseuille flow, and flow past a cylinder are shown below.
More details on setting up and running simulations can be found in the documentation.
Features
PyLaBolt currently supports the following collision schemes
- Bhatnagar-Gross-Krook (BGK) scheme - Physical Review, vol. 94, Issue 3, pp. 511-525
The MRT and TRT collision schemes will be added in future releases.
The boundary conditions available are
- No slip boundary via the halfway bounce back method - Journal of Fluid Mechanics , Volume 271 , 25 July 1994 , pp. 285 - 309
- Moving wall boundary condition via the halfway bounce back method - Journal of Fluid Mechanics , Volume 271 , 25 July 1994 , pp. 285 - 309 , Journal of Statistical Physics volume 104, pages 1191–1251 (2001)
- The fixed pressure boundary condition via the anti-bounce back method - Commun. Comput. Phys. 3, 427 (2008)
- Zero gradient boundary condition
- Periodic boundary condition - The Lattice Boltzmann Method
For more information on the schemes and boundary conditions, we urge the reader the go through the following books
- The Lattice Boltzmann Method - Timm Krüger, Halim Kusumaatmaja, Alexandr Kuzmin, Orest Shardt, Goncalo Silva, Erlend Magnus Viggen
- The Lattice Boltzmann Equation: For Complex States of Flowing Matter - Sauro Succi
PyLaBolt leverages the performance advantages on multi-core CPUs, High Perfomance computing clusters, and GPUs to run large simulations. Currently the parallel computing features supported by PyLaBolt are:
- Numba accelerated Python code can run in parallel on multi-core CPUs through OpenMP.
- To run on multiple machines/clusters, PyLaBolt uses OpenMPI via the mpi4py library.
- PyLaBolt can also run simulations on NVIDIA GPUs through Numba's CUDA bindings.
- PyLaBolt provides support to convert output to VTK format, which can post-processed in Paraview/Mayavi.
For more details, refer to the documentation.
Acknowledgements
PyLaBolt grew from the course project of AE6102 course, of the Department of Aerospace Engineering at Indian Institute of Technology, Bombay. We are grateful to Prof. Prabhu Ramachandran of Department of Aerospace Engineering at Indian Institute of Technology, Bombay for exposing us to the necessary skills to write high performance scientific codes.
We are also grateful to Dr. Amol Subhedar at Department of Chemical Engineering, Indian Institute of Technology, Bombay for providing us critical suggestions on the theoretical foundations of the algorithms.
For contributors see the Github contributors page.