In this piece of work, a numerical investigation of natural convection is carried out on the buoyancy-driven flow of nanofluids and heat transfer through porous media packed inside a wavy cavity. The cavity is placed horizontal, and its right vertical wall is of wavy nature, the bottom and top walls of the cavity are adiabatic, and there is a temperature difference between the left and right vertical wall. The dimensionless governing equations for the flow of nanofluids through the Darcian porous media are solved iteratively by using finite difference method. The study is conducted for wide range of governing parameters, such as Rayleigh-Darcy number [Formula: see text], nanoparticle volume fraction [Formula: see text] for three types of nanofluids [Formula: see text]-[Formula: see text], Cu-[Formula: see text], TiO2-[Formula: see text], the waviness of the vertical wall controlled by dimensionless length of amplitude of the wave [Formula: see text] and number of undulations per unit length ( N = 1, 3, 5). The simulated results reveals that the presence of nanoparticles enhances the convective heat transfer process at low Ra, and the wall affects the local convection rate and it also controls the overall heat transfer rate. For a cavity with N = 3, [Formula: see text] is increased by 33% at Ra = 10, and at [Formula: see text] has a drop by 10% as the a is increased from 0.05 to 0.25 having 20% of nanoparticles.