It is known that, when the immersed boundary method (IBM) is implemented within spectral-like methods, the Gibbs oscillation seriously deteriorates the calculation of derivatives near the body surface. In this paper, a radial basis function (RBF) based smoothing technique is proposed with the intention of eliminating or efficiently reducing the Gibbs oscillation without affecting the flow field outside the body. Based on this technique, a combined IBM/spectral scheme is developed to solve the incompressible Navier-Stokes equations. Numerical simulations of flow through a periodic lattice of cylinders of various cross sections are performed. The results demonstrate that the proposed methodology is able to give accurate and nearly oscillation-free numerical solutions of incompressible viscous flows. (C) 2011 Elsevier Inc. All rights reserved.