This study analyzed the difference in wave diffraction force of a surface piercing body according to the grid update scheme of free surface nodes in the time integration process. A two-dimensional fully nonlinear numerical wave tank based on potential flow theory and boundary element method was used for numerical analysis. The MEL (mixed Eulerian–Lagrangian) method was used for free surface node treatment. This method was divided into a semi-Lagrangian approach and a full Lagrangian approach (material node approach).The Runge–Kutta 4th-order integration method was used to integrate the free surface boundary condition temporally. In the internal sub-time step of the time integration process, the frozen coefficient scheme and the full update scheme were used according to the node update scheme, respectively. The vertical diffraction forces of each frequency component for node update schemes were compared. The difference in the 2nd-order vertical forces according to wave frequency, wave steepness, time step interval, and node spacing was investigated. A more precise time-domain numerical analysis method was proposed for high-order external forces in the wave-floating body interaction.