The Darcy flow model with the Boussinesq approximation is used to investigate numerically the natural convection inside an inclined wavy cavity filled with a porous medium. Finite Element Method is used to discretize the governing differential equations with non-staggered variable arrangement. Results are presented for 0 degrees <= phi <= 90 degrees, 10 <= Ra <= 10(3), 1 <= A <= 3 and 0 <= lambda <= 0.3, where phi, Ra, A and lambda correspond to the cavity inclination angle, Rayleigh number, aspect ratio and surface waviness parameter, respectively. Stream and isotherm lines representing the corresponding flow and thermal fields, and local and average Nusselt numbers distribution expressing the rate of heat transfer are determined and shown on graphs and tables. A good agreement is observed between the present results and those known from the open literature. The flow and thermal structures found to be highly dependent on surface waviness for inclination angles less than 45 degrees, especially for high Rayleigh numbers.