With the recent rapid progress of high-end computers, i.e., one million times faster than it was twenty years ago, it is expected that applications of fully-resolved large eddy simulation (LES) which directly computes the streamwise vortices in a turbulent boundary layer will become feasible in many engineering fields within a few years. The objective of this study is to investigate feasibility of applying fully-resolved LES to the prediction of model ship hydrodynamics. For this purpose, bare hull double model computations of KVLCC2 are carried out with one billion computational grids by using 1,536 computing cores. Feasible results are obtained at two different Reynolds numbers. Details of the turbulent boundary layer and the wake stemming from them are simulated well. Since the boundary layer is not tripped in the computations at Reynolds number of one million, laminar-turbulent transition and laminar effects are observed in detail by instantaneous vorticity distributions and the local skin friction coefficients.