Computational studies on entropy generation during laminar mixed convection in porous square enclosures have been carried out based on Darcy-Brinkman-Forchheimer model using the penalty finite element method. Finite element simulations are performed for the isothermally hot bottom wall, adiabatic top wall, and isothermally cold side walls (case 1) or linearly heated side walls (case 2) or linearly heated left wall with isothermally cold right wall (case 3) for a wide range of modified Darcy numbers (10 À5 Da m 10 À2 ), Grashof numbers (Gr ¼ 10 3 À 10 5 ), and modified Prandtl numbers (Pr m ¼ 0.026 and 7.2). Further, the effects of Da m on the total entropy generation (S total ), average Bejan number (Be av ), and average Nusselt number ðNu b Þ are discussed. It is found that Re ¼ 100 is preferred over Re ¼ 1 based on larger heat transfer rate with minimum entropy generation for Pr m ¼ 0.026 and 7.2, 10 À5 Da m 10 À2 at Gr ¼ 10 5 for all the cases.