Highly nonlinear extreme waves are the major, often the dominant, environmental load on offshore wind turbines. The higher-order ‘ringing’ loads associated with the nonlinear waves can cause unexpected resonance of the monopile. Hydrodynamic analysis of these harmonic loads remains a challenge due to the difficulty in extracting the bound harmonics from the force spectrum in an extreme wave event. A phase manipulation approach (four-phase combination) has been recently demonstrated to be able to separate the higher harmonic components of the wave loads in tank tests. In this work, we employ a fully nonlinear potential flow based Numerical Wave Tank (NWT) to simulate the wave diffraction by a fixed vertical column. We present a detailed study of our checks on the numerical accuracy of our model. Phase control is implemented for the wavemaker to manipulate the phase of each wave component. Focused wave groups are generated to represent the incoming extreme waves. With the four-phase decomposition, the higher harmonics of the wave loads are shown to be clearly separated. Comparisons with the existing test results show fairly good agreement at higher harmonics. The structure of the harmonic forces and moments are analysed and we reconstruct the higher harmonics based on the Stokes expansion assumption using the linear force. In addition, the effects of wave steepness on the harmonic components are discussed.