We present a residual-based isogeometric variational multiscale method to solve laminar and turbulent channel flow. Residual-based variational multiscale method is a new finite element formulation for solving turbulent flows using a large-eddy simulation type modeling. Isogeometric analysis, a new finite element method using CAD blending functions as its basis functions, is employed for higher order approximation of the solution. First, laminar flow with $Re_{ au}=0.55$ through flat channel is considered and linear, quadratic and cubic basis functions, which are $C^0$, $C^1$, and $C^2$ -continuous across element boundaries, respectively are employed and their accuracy is presented by the comparison with analytical result. Next, same methodology is applied to the turbulent channel flow with $Re_{ au}=180$. Current results are validated by the comparison of turbulence statistics using available DNS data.