The present study deals with the finite element based numerical simulations of heat transfer and entropy generation rates during natural convection for fluid saturated porous media in enclosures involving curved walls (case 1: lower curvature and case 2: higher curvature) with various thermal boundary conditions. The differential heating (isothermally hot left wall and cold right wall and adiabatic horizontal walls) and Rayleigh-Benard heating (isothermally hot bottom wall and cold top wall involving adiabatic left and right walls) are considered. The locations and magnitudes of the entropy generation due to heat transfer (S-theta) and fluid friction (S-psi) are presented and discussed based on the spatial distributions of isotherms and streamlines, respectively. The magnitudes of local entropy generation (S-theta, S-psi), total entropy generation (S-total) and average heat transfer rates ((N-ur) over bar and (N-ut) over bar) are significantly lesser for the Rayleigh-Benard heating compared to the differential heating for all the cases involving all Da(m) and Pr-m. The Rayleigh-Benard heating is the optimal strategy for all Da(m), and Pr-m, involving both the concave cases except for 10(-3) <= Da(m) <= 10(-2), Pr-m = 10 and case 1 (concave) domain. The Rayleigh-Benard heating is also the optimal strategy compared to the differential heating involving the convex cases at 10(-5) <= Da(m) <= 10(-4) whereas the differential heating is the optimal heating strategy for Da(m) >= 10(-3) involving both Pr-m for the convex cases. (C) 2018 Elsevier Ltd. All rights reserved.