As the first step to investigate the nonlinear interactions between turbulence and marine structures inside a viscous NWT, a LES technique was applied to solve the turbulent channel flow for =150. The employed turbulence models included 4 types: the Smagorinsky model, the Dynamic SGS model, the Structure Function model, and the Generalized Normal Stress model. The simulated data in time-series for the LESs were averaged in both time and space, and statistical analyses were performed. The results of the LESs were compared with those of a DNS, developed in the present study and two spectral methods by Yoon et al.(2003) and Kim et a1.(1987). Based on this research, the accuracy of LESs has been found to be still related to the number of grids for fine grid size).