High-fidelity models become more and more used in the wave energy sector. They offer a fully nonlinear simulation tool that in theory should encompass all linear and nonlinear forces acting on a wave energy converter (WEC). Studies using high-fidelity models are usually focusing on validation of the model. However, a validated model does not necessarily provide reliable solutions. Solution verification is the methodology to estimate the numerical uncertainties related to a simulation. In this work we test four different approaches: the classical grid convergence index (GCI); a least-squares version (LS-GCI); a simplified version of the least-square method (SLS-GCI); and the ITTC recommended practice. The LS-GCI requires four or more solutions whereas the other three methods only need three solutions. We apply these methods to four different high-fidelity models for the case of a heaving sphere. We evaluate the numerical uncertainties for two parameters in the time domain and two parameters in the frequency domain. It was found that the GCI and ITTC were hard to use on the frequency domain parameters as they require monotonic convergence which sometimes does not happen due to the differences in the solutions being very small. The SLS-GCI performed almost as well as the SL-GCI method and will be further investigated.