A numerical study is performed to analyze steady state forced convection fluid flow through a confined porous square cylinder. The Darcy-Brinkman-Forchheimer model is adopted for the porous region. The finite volume method and the iterative SIMPLE algorithm are used to solve the governing equations. The results obtained are presented for the streamlines, variation of Nusselt number and drag coefficient for the range of conditions as 5 ≤ Re ≤ 40 and 10−2 ≤ Da ≤ 10−6.