The LBM operates on a regular lattice and decomposes the fluid domain into a set of lattice nodes. The fluid is modeled as a group of fluid particles that are only allowed to move between lattice nodes or stay at rest. The composition of the lattice nodes depends on the chosen lattice model. The most commonly used lattice model for 2D simulations is the one using a square lattice with nine discrete velocity directions (denoted as model D2Q9, shown in Fig. 1).The motion of fluid particles is governed by the discrete lattice Boltzmann equation (LBE).