A numerical investigation was conducted to analyze free convective movement inside a square shape porous inclusion with a curved left wall characterized by varying degrees of waviness. The left curvy wall was sustained at a substantially low temperature, while the lower boundary was heated sinusoidally, and the right vertical boundary was subjected to a linear cooling. This research is predominantly pertinent to industries such as agro-food processing, glass manufacturing, insulation of thermal systems, and cooling of microprocessors, where precise thermal forecasts are crucial. This study examines circulation patterns, local Nusselt numbers, and isothermal lines, while exploring the impacts of key parameters, comprising Rayleigh number (103 ≤ Ra ≤ 10⁶), number of waviness (2, 4, 8), and Darcy number (10⁻4 ≤ Da ≤ 10⁻1). Governing equations were numerically solved using Ansys Fluent employing SIMPLE algorithm for coupling velocity and pressure terms. Current study indicates that the rate of heat transmission significantly increases with higher Da and Ra. Heat transfer has also been found to augment with the number of waviness in the left vertical wall.