The present paper deals with the numerical investigation of a 2D laminar fluid flow and heat transfer in a plane channel with two square blocks located at arbitrary positions. The numerical model is based on a coupling between the multiple relaxation time-lattice Boltzmann equation and the finite difference method for incompressible flow. Both the horizontal and the vertical separation distances between the two blocks are varied. Particular attention was paid to the distribution patterns of the time averaged local Nusselt number on the top and bottom walls. Results obtained from the present study show a complex flow patterns developed in the channel due to the change of the square blocks positions.