In this paper, a numerical simulation technique is developed to investigate the qualitative and quantitative behaviour of Cu-nanoparticles in a porous medium vis-a-vis the heat transfer enhancements-buoyancy driven flow in a two-dimensional square cavity, with moving walls is presented. The model utilizes the finite volume approach to solve the Brinkman-Darcy equations for Cu-nanoparticles in a porous media. Discretization is carried out for convective and diffusive fluxes using Quadratic Upwind Interpolation for Convective Kinematics (QUICK) and central difference schemes, respectively. Tri-Diagonal Matrix Algorithm is invoked to solve the set of algebraic equations. The Darcy number (Da), Prandtl number (Pr), and volume fraction (χ) are varied from 10-3 to 10-1, 3 to 7, and 0% to 20%, respectively. Insight into the cause of variations in isotherms, streamlines, Nusselt number (Nu), and mid-plane velocities is explicated. The present numerical results are compared with the existing literature and found to be in good agreement. Even though nanoparticles slightly hinder the activity of the fluid, they can augment the average Nu by 90% for Pr = 7, Da = 0.1, and χ = 20% as compared to the absence of nanoparticles. Their efficacy is more prominent for flows with higher Da and Pr. Quantitative values for Nu were obtained for various combinations of Pr, Da, and χ. © 2014 Wiley Periodicals, Inc.