Here some simple codes for the Lattice Boltzmann Method are presented. Along with the codes I as well present some simple Python scripts to analyze results produced by those codes.

1. Shan-Chen gas-liquid model with the droplet in the center of a domain.

2. Binary-liquid free energy model with the droplet on the surface.

3. Shan-Chen gas-liquid model programmed on GPU with OpenCL technology.

4. Binary liquid model of the long gas bubbles in 3D microchannels (Bretherton flow). The implementation is parallelized in MPI. Two different implementations are covering full scale simulation and quarter channel simulations with symmetry conditions.

5. The binary liquid droplets on the curved surface. The boundary condition for a phase field is established in a such way as to control the wetting angle.

6. Deposition of the droplet to the curved solid.

7. Propagation of the Gaussian hill via advection-diffusion equation modeled by two-relaxation-times LB model.

8. 2D multiple-relaxation-times codes with two different approaches: through Gram-Schmidt procedure and through Hermite polynomials approach.

9. The free-surface code that reads velocities and geometry of the bubble from files, and then imposes the velocity of the walls that bubble stays in the center (reference frame of the bubble). The bubble surface is the subject to the free surface boundary condition, that no normal velocity components do exist. That allows to cleanup the LB multiphase simulation to get rid of non-zero divergence at the bubble surface to be able to perform proper mass transfer calculations. Otherwise non-zero divergence at the surface of the bubble makes mass transfer simulations unstable leading to an infinite accumalation of mass.

10. The diffusion from the boundary of 2D cylinder into the inside cylinder domain. Two boundary conditions of the constant concentration are presented: Inamuro and anti-bounce-back conditions.

11. The thermal layer development from the wall with the uniform velocity profile.

12. The thermal layer development from the wall with the Poiseuille velocity profile.

13. Steady cylinder inside the channel simulated by immersed boundary method.

14. Moving cylinder simulated with the immersed boundary method. The cylinder is located off center-line inside the channel with the moving walls into the opposite directions. The cylinder is moving to the center of the channel due to lift forces.

15. Falling ellipse by gravity simulated with the immersed boundary method.