For the fluid domain \(\Omega \) , the lattice Boltzmann (LB) equations are employed, where the synchronous motions of fluid particles on a regular lattice are enforced through a particle distribution function [1]}. This distribution function enforces mass and momentum conservation as well as ensuring that the fluid is Galilean invariant and isotropic [2]}. In the present work, a single-relaxation-time (SRT) incompressible LB method is used to solve the incompressible flow [3]}. The evolution of the distribution functions on the lattice is governed by the discrete Boltzmann equation with the Bhatnagar-Gross-Krook (BGK) collision model, given by
\(f_i(x+e_i\Delta t,t+\Delta t)-f_i(x,t)=-\frac{1}{\tau }[f_i(x,t)-f_i^{eq}(x,t)]+\Delta t F_i(x,t), \quad i=0,...,N_0-1,\) 
