Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${alpha}/{lambda}$) of wavy wall is 0.05. The computational domain is $2{lambda}{ imes}1.05{lambda}{ imes}{lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.