This study explores the numerical investigation of two-dimensional natural convection in a porous corrugated enclosure with heat sources of different heights and electrically conducting fluids. The governing equations are solved using a higher-order compact finite difference scheme. Boundaries on the top and sides are subject to thermal cooling and adiabatic conditions, respectively. Numerical results are presented for wide range of Rayleigh (103 to 106), Darcy (10−5 to 10−1), Hartmann (25 to 150), and Prandtl (0.015 to 10) numbers. For fixed high Rayleigh-Darcy values, the porous medium has higher permeability. Most significant possible multiple steady solutions can be found in the case 1 at Ra=106,Da=10−2,Ha=50, and Pr=10.0. In the low to moderate critical values of Prandtl numbers, it encompasses a wide range of fluid flow phenomena, from stable to unstable. The maximum enhancement of the averaged Nusselt numbers on the boundaries in increasing and decreasing trends are viz. Ra(139.13% and 54.37%), Da(46.14% and 74.25%), Ha(633.04% and 74.73%), and Pr(142.52% and 45.92%) for case 1, Ra(129.64% and 0.32%), Da(40.68% and 69.66%), Ha(46.65% and 80.85%), and Pr(18.41% and 18.23%) for case 2, and Ra(163.98% and 77.20%), Da(8.25% and 58.68%), Ha(81.38% and 59.17%), and Pr(30.60% and 22.46%) for case 3.