Lattice Boltzmann Method (LBM) is widely accepted for solving multi-phase fluid interaction applications.The LBM – a mesoscale method – doesn’t require to solve governing equations at molecular level (unlikemolecular dynamics-based methods) which makes it computationally lesser costlier. In addition, it is also taking care of the physical properties for both the fluids. LBM with multi-relaxation-time (MRT) model is having higher stability and accuracy, and hence, LBM with MRT-based methodology is implemented in this work. LBM-based multi-phase fluid solver is developed, however in this work only single-phase benchmarking problems are solved in order to validate the proposed solver. 2-D channel flow in backward-facing step and 2-D channel flow over a square cylinder are solved, and results are validated through available numerical and analytical solutions available. The problem of backward-facing step is simulated at \(R_{e} = 800\) and flow over a square cylinder is simulated at two different Reynolds number, i.e., \(R_{e} = 60\,\& \,100\) . Efforts are made to showcase the capabilities and benefits of LB-based model in comparison of well-established Navier–Stokes-based solvers. Numerical results of the proposed work are in nice matching with well-established NS-based numerical and analytical solutions available in the literature.