Motivated by petro-chemical geological systems, we consider the natural convection boundary layer flow from a vertical isothermal wavy surface adjacent to a saturated non-Darcian high permeability porous medium. High permeability is considered to represent geologically sparsely packed porous media. Both Darcian drag and Forchheimer inertial drag terms are included in the velocity boundary layer equation. A high permeability medium is considered. We employ a sinusoidal relation for the wavy surface. Using a set of transformations, the momentum and heat conservation equations are converted from an (x, y) coordinate system to an (x,η) dimensionless system. The two-point boundary value problem is then solved numerically with a pseudo-spectral method based on combining the Bellman–Kalaba quasi linearization method with the Chebyschev spectral collocation technique (SQLM). The SQLM computations are demonstrated to achieve excellent correlation with smoothed particle hydrodynamic (SPH) Lagrangian solutions. We study the effect of Darcy number (Da), Forchheimer number (Fs), amplitude wavelength (A) and Prandtl number (Pr) on the velocity and temperature distributions in the regime. Local Nusselt number is also computed for selected cases. The study finds important applications in petroleum engineering and also energy systems exploiting porous media and undulating (wavy) surface geometry. The SQLM algorithm is shown to be exceptionally robust and achieves fast convergence and excellent accuracy in nonlinear heat transfer simulations
