Skip to content

zdxying/FreeLB

Repository files navigation

FreeLB: Freely Coupled Lattice Boltzmann Code

This code(FreeLB) aims to facilitate the free dissemination and advancement of knowledge regarding the lattice Boltzmann (LB) coupled with other methods/models(cellular automata (CA), free surface, non-newtonian fluid...).

FreeLB should be helpful for beginners who know the basic of LBM but have not started to implement it yet, especially, for those who want to build their own LBM program.

The following picture shows the capability of FreeLB to couple with CA method to simulate the solidification process of a dendrite with block structure

Concentration_Field

The following picture shows the developed feature of adaptive mesh refinement(AMR) in FreeLB

Partof_Concentration_Field

Streamline

The author is inspired by OpenLB project, but provides a more concise and clear code structure, also, a relatively direct data access pattern, which is more suitable for beginners to understand and modify the code. The following files are modified from OpenLB project:

  • 'src/data_struct/octree.h'
  • 'src/data_struct/octree.hh'
  • 'src/data_struct/Vector.h'
  • 'src/io/base64.h'
  • 'src/io/stlreader.h'
  • 'src/io/stlreader.hh'
  • 'src/parallel/mpi_manager.h'
  • 'src/parallel/mpi_manager.hh'
  • 'src/lbm/freeSurface.h'

This code is written in C++ for clarity and efficiency, methods of LB and CA are stored in classes in ./src folder, specificlly, the declaration of the methods are in .h files and the implementations are in .hh files. To implement the simulation, one should prepare an .ini file storing simulation parameters and write his/her own .cpp/.cu file which calls the methods of this code, also, defines the initial and boundary conditions. There are some examples in ./examples folder, which can be used as a reference.

OpenMP is enabled through the code, MPI is available for most of methods and Cuda is enabled for a few methods, which is still under development.

If you have any suggestions, please contact the author by email: ymmanyuan@outlook.com

Build

To build the code, you just need to run the following command in the root directory of the code:

cd ./examples/cavblock3d_benchmark
make

The executable file will be generated in the same directory with .exe extension.

Benchmark

The following table shows the performance of FreeLB in simulating the lid-driven cavity flow problem(/examples/cavblock3d_benchmark and /examples/cavblock3d_benchmarkcu) with D3Q19 lattice set on a 100x100x100 lattice grid(in million lattice updates per second, MLUPs).

Device Platform FP32 MLUPs
i9-14900K Ubuntu 22.04 608
E5-2699v4 Ubuntu 20.04 WSL2 553
i7-11800H Ubuntu 22.04 WSL2 173
i7-8750H Ubuntu 22.04 WSL2 124
RTX 4080 Ubuntu 22.04 2903
RTX 3060 Laptop Ubuntu 22.04 WSL2 1581

Dependencies

  • C++ compiler supporting C++17
  • GNU Make
    Optional(If parallel computing is needed, install):
  • OpenMP
  • MPI
  • Cuda toolkit

License

FreeLB is licensed under the GNU GPL v3 license.

Acknowledgement

The author would like to thank OpenLB project for the inspiration of this code and thank Prof. Tao Jing for his help in the development of this code.