In the project a literature study [1] was first carried out to collect the state-of-the-art numerical modelling techniques using CFD methods. This was followed by a validation study where CFD results for an LNG vessel were compared with SSPA’s model test data. As a result, a modelling strategy and numerical schemes were proposed for the present work. Systematic simulations were then carried out for targeted ships at full scale, with variation of a few key parameters defined in a simulation matrix. The report presents the main findings from the computational results carried out so far for the first targeted vessel, Tanker-1. The study clearly indicates that the seabed is subjected to a higher bottom shear stress (BBS) when the following parameter changes: at a higher ship speed, a fuller ship with higher block coefficient, or at a lower water depth h/T. The maximum BBS occurs underneath the fore shoulder of ship. The current study demonstrates that it is feasible to predict the bottom shear stress and shallow water ship resistance with the proposed CFD methodology. The increase of resistance and the characteristic distribution of bottom shear stress are fully in line with the theoretical expectations. Though the absolute accuracy level of the method is subjected to questions, for engineering applications this level of accuracy may still be enough. The full scale computation with VOF module takes substantially long time. Nevertheless, considering the relatively low cost of using CFD tool with the costly on-site measurement of BSS, the CFD method offers a relatively efficient tool to predict these quantities. Future work should include the propeller effect and the effect of bottom roughness, for each studied ship. Furthermore, the CFD calculation for the rest of 7 vessels defined in the simulation matrix should be continued and completed before a sufficient database is acquired for regression analysis.